Evolution of prokaryotic and eukaryotic virulence effectors

Ecological dependencies and the illusion of cooperation in microbial communities

Ecological dependencies - where organisms rely on other organisms for survival - are a ubiquitous feature of life on earth. Multicellular hosts rely on symbionts to provide essential vitamins and amino acids. Legume plants similarly rely on nitrogen-fixing rhizobia to convert atmospheric nitrogen to ammonia. In some cases, dependencies can arise via loss-of-function mutations that allow one partner to benefit from the actions of another. It is common in microbiology to label ecological dependencies between species as cooperation - making it necessary to invoke cooperation-specific frameworks to explain the phenomenon. However, in many cases, such traits are not (at least initially) cooperative, because they are not selected for because of the benefits they confer on a partner species. In contrast, dependencies in microbial communities may originate from fitness benefits gained from genomic-streamlining (i.e. Black Queen Dynamics). Here, we outline how the Black Queen Hypothesis predicts the formation of metabolic dependencies via loss-of-function mutations in microbial communities, without needing to invoke any cooperation-specific explanations. Furthermore we outline how the Black Queen Hypothesis can act as a blueprint for true cooperation as well as discuss key outstanding questions in the field. The nature of interactions in microbial communities can predict the ability of natural communities to withstand and recover from disturbances. Hence, it is vital to gain a deeper understanding of the factors driving these dynamic interactions over evolutionary time.

Pathogen protein modularity enables elaborate mimicry of a host phosphatase

Pathogens produce diverse effector proteins to manipulate host cellular processes. However, how functional diversity is generated in an effector repertoire is poorly understood. Many effectors in the devastating plant pathogen Phytophthora contain tandem repeats of the "(L)WY" motif, which are structurally conserved but variable in sequences. Here, we discovered a functional module formed by a specific (L)WY-LWY combination in multiple Phytophthora effectors, which efficiently recruits the serine/threonine protein phosphatase 2A (PP2A) core enzyme in plant hosts. Crystal structure of an effector-PP2A complex shows that the (L)WY-LWY module enables hijacking of the host PP2A core enzyme to form functional holoenzymes. While sharing the PP2A-interacting module at the amino terminus, these effectors possess divergent C-terminal LWY units and regulate distinct sets of phosphoproteins in the host. Our results highlight the appropriation of an essential host phosphatase through molecular mimicry by pathogens and diversification promoted by protein modularity in an effector repertoire.

The absence of PNPase activity in Enterococcus faecalis results in alterations of the bacterial cell-wall but induces high proteolytic and adhesion activities

Enterococci are lactic acid bacteria (LAB) used as starters and probiotics, delineating their positive attributes. Nevertheless, enterococci can be culprit for thousands of infectious diseases, including urinary tract infections, bacteremia and endocarditis. Here, we aim to determine the impact of polynucleotide phosphorylase (PNPase) in the biology of Enterococcus faecalis 14; a human isolate from meconium. Thus, a mutant strain deficient in PNPase synthesis, named ΔpnpA mutant, was genetically obtained. After that, a transcriptomic study revealed a set of 244 genes differentially expressed in the ΔpnpA mutant compared with the wild-type strain, when exploiting RNAs extracted from these strains after 3 and 6 h of growth. Differentially expressed genes include those involved in cell wall synthesis, adhesion, biofilm formation, bacterial competence and conjugation, stress response, transport, DNA repair and many other functions related to the primary and secondary metabolism of the bacteria. Moreover, the ΔpnpA mutant showed an altered cell envelope ultrastructure compared with the WT strain, and is also distinguished by a strong adhesion capacity on eukaryotic cell as well as a high proteolytic activity. This study, which combines genetics, physiology and transcriptomics enabled us to show further biological functions that could be directly or indirectly controlled by the PNPase in E. faecalis 14.

Transcriptional, hormonal, and metabolic changes in susceptible grape berries under powdery mildew infection

Grapevine (Vitis vinifera) berries are extremely sensitive to infection by the biotrophic pathogen Erysiphe necator, causing powdery mildew disease with deleterious effects on grape and wine quality. The combined analysis of the transcriptome and metabolome associated with this common fungal infection has not been previously carried out in any fruit. In order to identify the molecular, hormonal, and metabolic mechanisms associated with infection, healthy and naturally infected V. vinifera cv. Carignan berries were collected at two developmental stages: late green (EL33) and early véraison (EL35). RNA sequencing combined with GC-electron impact ionization time-of-flight MS, GC-electron impact ionization/quadrupole MS, and LC-tandem MS analyses revealed that powdery mildew-susceptible grape berries were able to activate defensive mechanisms with the involvement of salicylic acid and jasmonates and to accumulate defense-associated metabolites (e.g. phenylpropanoids, fatty acids). The defensive strategies also indicated organ-specific responses, namely the activation of fatty acid biosynthesis. However, defense responses were not enough to restrict fungal growth. The fungal metabolic program during infection involves secretion of effectors related to effector-triggered susceptibility, carbohydrate-active enzymes and activation of sugar, fatty acid, and nitrogen uptake, and could be under epigenetic regulation. This study also identified potential metabolic biomarkers such as gallic, eicosanoic, and docosanoic acids and resveratrol, which can be used to monitor early stages of infection.

Sec-Delivered Effector 1 (SDE1) of 'Candidatus Liberibacter asiaticus' Promotes Citrus Huanglongbing

Sec-delivered effector 1 (SDE1) from the huanglongbing (HLB)-associated bacterium 'Candidatus Liberibacter asiaticus' was previously characterized as an inhibitor of defense-related, papain-like cysteine proteases in vitro and in planta. Here, we investigated the contributions of SDE1 to HLB progression. We found that SDE1 expression in the model plant Arabidopsis thaliana caused severe yellowing in mature leaves, reminiscent of both 'Ca. L. asiaticus' infection symptoms and accelerated leaf senescence. Induction of senescence signatures was also observed in the SDE1-expressing A. thaliana lines. These signatures were apparent in older leaves but not in seedlings, suggesting an age-associated effect. Furthermore, independent lines of transgenic Citrus paradisi (L.) Macfadyen (Duncan grapefruit) that express SDE1 exhibited hypersusceptibility to 'Ca. L. asiaticus'. Similar to A. thaliana, transgenic citrus expressing SDE1 showed altered expression of senescence-associated genes, but only after infection with 'Ca. L. asiaticus'. These findings suggest that SDE1 is a virulence factor that contributes to HLB progression, likely by inducing premature or accelerated senescence in citrus. This work provides new insight into HLB pathogenesis.[Formula: see text] Copyright © 2020 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

autoBioSeqpy: A Deep Learning Tool for the Classification of Biological Sequences

Deep learning has proven to be a powerful method with applications in various fields including image, language, and biomedical data. Thanks to the libraries and toolkits such as TensorFlow, PyTorch, and Keras, researchers can use different deep learning architectures and data sets for rapid modeling. However, the available implementations of neural networks using these toolkits are usually designed for a specific research and are difficult to transfer to other work. Here, we present autoBioSeqpy, a tool that uses deep learning for biological sequence classification. The advantage of this tool is its simplicity. Users only need to prepare the input data set and then use a command line interface. Then, autoBioSeqpy automatically executes a series of customizable steps including text reading, parameter initialization, sequence encoding, model loading, training, and evaluation. In addition, the tool provides various ready-to-apply and adapt model templates to improve the usability of these networks. We introduce the application of autoBioSeqpy on three biological sequence problems: the prediction of type III secreted proteins, protein subcellular localization, and CRISPR/Cas9 sgRNA activity. autoBioSeqpy is freely available with examples at https://github.com/jingry/autoBioSeqpy.

Type III Secretion Effectors with Arginine N-Glycosyltransferase Activity

Type III secretion systems are used by many Gram-negative bacterial pathogens to inject proteins, known as effectors, into the cytosol of host cells. These virulence factors interfere with a diverse array of host signal transduction pathways and cellular processes. Many effectors have catalytic activities to promote post-translational modifications of host proteins. This review focuses on a family of effectors with glycosyltransferase activity that catalyze addition of N-acetyl-d-glucosamine to specific arginine residues in target proteins, leading to reduced NF-κB pathway activation and impaired host cell death. This family includes NleB from Citrobacter rodentium, NleB1 and NleB2 from enteropathogenic and enterohemorrhagic Escherichia coli, and SseK1, SseK2, and SseK3 from Salmonella enterica. First, we place these effectors in the general framework of the glycosyltransferase superfamily and in the particular context of the role of glycosylation in bacterial pathogenesis. Then, we provide detailed information about currently known members of this family, their role in virulence, and their targets.

Community diversity and habitat structure shape the repertoire of extracellular proteins in bacteria

We test the hypothesis that the frequency and cost of extracellular proteins produced by bacteria, which often depend on cooperative processes, vary with habitat structure and community diversity. The integration of the environmental distribution of bacteria (using 16S datasets) and their genomes shows that bacteria living in more structured habitats encode more extracellular proteins. In contrast, the effect of community diversity depends on protein function: it’s positive for proteins implicated in antagonistic interactions and negative for those involved in nutrient acquisition. Extracellular proteins are costly and endure stronger selective pressure for low cost and for low diffusivity in less structured habitats and in more diverse communities. Finally, Bacteria found in multiple types of habitats, including host-associated generalists, encode more extracellular proteins than niche-restricted bacteria. These results show that ecological variables, notably habitat structure and community diversity, shape the evolution of the repertoires of genes encoding extracellular proteins and thus affect the ability of bacteria to manipulate their environment.

Microbes secrete a repertoire of extracellular proteins to serve various functions depending on the ecological context. Here the authors examine how bacterial community composition and habitat structure affect the extracellular proteins, showing that generalist species and those living in more structured environments produce more extracellular proteins, and that costs of production are lower in more diverse communities.

The RXLR Effector PcAvh1 Is Required for Full Virulence of Phytophthora capsici

Plant pathogens employ diverse secreted effector proteins to manipulate host physiology and defense in order to foster diseases. The destructive Phytophthora pathogens encode hundreds of cytoplasmic effectors, which are believed to function inside the plant cells. Many of these cytoplasmic effectors contain the conserved N-terminal RXLR motif. Understanding the virulence function of RXLR effectors will provide important knowledge of Phytophthora pathogenesis. Here, we report the characterization of RXLR effector PcAvh1 from the broad-host range pathogen Phytophthora capsici. Only expressed during infection, PcAvh1 is quickly induced at the early infection stages. CRISPR/Cas9-knockout of PcAvh1 in P. capsici severely impairs virulence while overexpression enhances disease development in Nicotiana benthamiana and bell pepper, demonstrating that PcAvh1 is an essential virulence factor. Ectopic expression of PcAvh1 induces cell death in N. benthamiana, tomato, and bell pepper. Using yeast two-hybrid screening, we found that PcAvh1 interacts with the scaffolding subunit of the protein phosphatase 2A (PP2Aa) in plant cells. Virus-induced gene silencing of PP2Aa in N. benthamiana attenuates resistance to P. capsici and results in dwarfism, suggesting that PP2Aa regulates plant immunity and growth. Collectively, these results suggest that PcAvh1 contributes to P. capsici infection, probably through its interaction with host PP2Aa.

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.

Where are we going with genomics in plant pathogenic bacteria?

Genome sequencing is commonly used in research laboratories right now thanks to the rise of high-throughput sequencing with higher speed and output-to-cost ratios. Here, we summarized the application of genomics in different aspects of plant bacterial pathosystems. Genomics has been used in studying the mechanisms of plant-bacteria interactions, and host specificity. It also helps with taxonomy, study of non-cultured bacteria, identification of causal agent, single cell sequencing, population genetics, and meta-transcriptomic. Overall, genomics has significantly improved our understanding of plant-microbe interaction.

Post-translational Mechanisms of Host Subversion by Bacterial Effectors

Bacterial effector proteins are a specialized class of secreted proteins that are translocated directly into the host cytoplasm by bacterial pathogens. Effector proteins have diverse activities and targets, and many mediate post-translational modifications of host proteins. Effector proteins offer potential in novel biotechnological and medical applications as enzymes that may modify human proteins. Here, we discuss the mechanisms used by effectors to subvert the human host through blocking, blunting, or subverting immune mechanisms. This capacity allows bacteria to control host cell function to support pathogen survival, replication and dissemination to other hosts. In addition, we highlight that knowledge of effector protein activity may be used to develop chemical inhibitors as a new approach to treat bacterial infections.

Rapid turnover of effectors in grass powdery mildew (Blumeria graminis)

Background

Grass powdery mildew (Blumeria graminis, Ascomycota) is a major pathogen of cereal crops and has become a model organism for obligate biotrophic fungal pathogens of plants. The sequenced genomes of two formae speciales (ff.spp.), B.g. hordei and B.g. tritici (pathogens of barley and wheat), were found to be enriched in candidate effector genes (CEGs). Similar to other filamentous pathogens, CEGs in B. graminis are under positive selection. Additionally, effectors are more likely to have presence-absence polymorphisms than other genes among different strains.

Results

Here we identified effectors in the genomes of three additional host-specific lineages of B. graminis (B.g. poae, B.g. avenae and B.g. infecting Lolium) which diverged between 24 and 5 million years ago (Mya). We found that most CEGs in B. graminis are clustered in families and that most families are present in both reference genomes (B.g. hordei and B.g. tritici) and in the genomes of all three newly annotated lineages. We identified conserved protein domains including a novel lipid binding domain. The phylogenetic analysis showed that frequent gene duplications and losses shaped the diversity of the effector repertoires of the different lineages through their evolutionary history. We observed several lineage-specific expansions where large clades of CEGs originated in only one lineage from a single gene through repeated gene duplications. When we applied a birth-death model we found that the turnover rate (the rate at which genes are deleted and duplicated) of CEG families is much higher than for non-CEG families. The analysis of genomic context revealed that the immediate surroundings of CEGs are enriched in transposable elements (TE) which could play a role in the duplication and deletion of CEGs.

Conclusions

The CEG repertoires of related pathogens diverged dramatically in short evolutionary times because of rapid turnover and of positive selection fixing non-synonymous mutations. While signatures of positive selection on effector sequences are the expected outcome of the evolutionary “arms race” between pathogen and plant immune system, it is more difficult to infer the mechanisms and evolutionary forces that maintained an extreme turnover rate in CEG families of B. graminis for several millions of years.

Electronic supplementary material

The online version of this article (10.1186/s12862-017-1064-2) contains supplementary material, which is available to authorized users.

Dissection of genomic features and variations of three pathotypes of Puccinia striiformis through whole genome sequencing

Stripe rust of wheat, caused by Puccinia striiformis f. sp. tritici, is one of the important diseases of wheat. We used NGS technologies to generate a draft genome sequence of two highly virulent (46S 119 and 31) and a least virulent (K) pathotypes of P. striiformis from the Indian subcontinent. We generated ~24,000-32,000 sequence contigs (N50;7.4-9.2 kb), which accounted for ~86X-105X sequence depth coverage with an estimated genome size of these pathotypes ranging from 66.2-70.2 Mb. A genome-wide analysis revealed that pathotype 46S 119 might be highly evolved among the three pathotypes in terms of year of detection and prevalence. SNP analysis revealed that ~47% of the gene sets are affected by nonsynonymous mutations. The extracellular secreted (ES) proteins presumably are well conserved among the three pathotypes, and perhaps purifying selection has an important role in differentiating pathotype 46S 119 from pathotypes K and 31. In the present study, we decoded the genomes of three pathotypes, with 81% of the total annotated genes being successfully assigned functional roles. Besides the identification of secretory genes, genes essential for pathogen-host interactions shall prove this study as a huge genomic resource for the management of this disease using host resistance.

YopJ Family Effectors Promote Bacterial Infection through a Unique Acetyltransferase Activity

Gram-negative bacterial pathogens rely on the type III secretion system to inject virulence proteins into host cells. These type III secreted "effector" proteins directly manipulate cellular processes to cause disease. Although the effector repertoires in different bacterial species are highly variable, the Yersinia outer protein J (YopJ) effector family is unique in that its members are produced by diverse animal and plant pathogens as well as a nonpathogenic microsymbiont. All YopJ family effectors share a conserved catalytic triad that is identical to that of the C55 family of cysteine proteases. However, an accumulating body of evidence demonstrates that many YopJ effectors modify their target proteins in hosts by acetylating specific serine, threonine, and/or lysine residues. This unique acetyltransferase activity allows the YopJ family effectors to affect the function and/or stability of their targets, thereby dampening innate immunity. Here, we summarize the current understanding of this prevalent and evolutionarily conserved type III effector family by describing their enzymatic activities and virulence functions in animals and plants. In particular, the molecular mechanisms by which representative YopJ family effectors subvert host immunity through posttranslational modification of their target proteins are discussed.

Comparative Genomic Analysis among Four Representative Isolates of Phytophthora sojae Reveals Genes under Evolutionary Selection

Comparative genomic analysis is useful for identifying genes affected by evolutionary selection and for studying adaptive variation in gene functions. In Phytophthora sojae, a model oomycete plant pathogen, the related study is lacking. We compared sequence data among four isolates of P. sojae, which represent its four major genotypes. These isolates exhibited >99.688%, >99.864%, and >98.981% sequence identities at genome, gene, and non-gene regions, respectively. One hundred and fifty-three positive selection and 139 negative selection candidate genes were identified. Between the two categories of genes, the positive selection genes were flanked by larger intergenic regions, poorly annotated in function, and less conserved; they had relatively lower transcription levels but many genes had increased transcripts during infection. Genes coding for predicted secreted proteins, particularly effectors, were overrepresented in positive selection. Several RxLR effector genes were identified as positive selection genes, exhibiting much stronger positive selection levels. In addition, candidate genes with presence/absence polymorphism were analyzed. This study provides a landscape of genomic variation among four representative P. sojae isolates and characterized several evolutionary selection-affected gene candidates. The results suggest a relatively covert two-speed genome evolution pattern in P. sojae and will provide clues for identification of new virulence factors in the oomycete plant pathogens.

Molecular Soybean-Pathogen Interactions

Soybean hosts a wide variety of pathogens that cause significant yield losses. The importance of soybean as a major oilseed crop has led to research focused on its interactions with pathogens, such as Soybean mosaic virus, Pseudomonas syringae, Phytophthora sojae, Phakopsora pachyrhizi, and Heterodera glycines. Pioneering work on soybean's interactions with these organisms, which represent the five major pathogen groups (viruses, bacteria, oomycetes, fungi, and nematodes), has contributed to our understanding of the molecular mechanisms underlying virulence and immunity. These mechanisms involve conserved and unique features that validate the need for research in both soybean and homologous model systems. In this review, we discuss identification of effectors and their functions as well as resistance gene-mediated recognition and signaling. We also point out areas in which model systems and recent advances in resources and tools have provided opportunities to gain deeper insights into soybean-pathogen interactions.

Weeds, as ancillary hosts, pose disproportionate risk for virulent pathogen transfer to crops

Background

The outcome of the arms race between hosts and pathogens depends heavily on the interactions between their genetic diversity, population size and transmission ability. Theory predicts that genetically diverse hosts will select for higher virulence and more diverse pathogens than hosts with low genetic diversity. Cultivated hosts typically have lower genetic diversity and thus small effective population sizes, but can potentially harbour large pathogen population sizes. On the other hand, hosts, such as weeds, which are genetically more diverse and thus have larger effective population sizes, usually harbour smaller pathogen population sizes. Large pathogen population sizes may lead to more opportunities for mutation and hence more diverse pathogens. Here we test the predictions that pathogen neutral genetic diversity will increase with large pathogen population sizes and host diversity, whereas diversity under selection will increase with host diversity. We assessed and compared the diversity of a fungal pathogen, Rhynchosporium commune, on weedy barley grass (which have a large effective population size) and cultivated barley (low genetic diversity) using microsatellites, effector locus nip1 diversity and pathogen aggressiveness in order to assess the importance of weeds in the evolution of the neutral and selected diversity of pathogens.

Results

The findings indicated that the large barley acreage and low host diversity maintains higher pathogen neutral genetic diversity and lower linkage disequilibrium, while the weed maintains more pathotypes and higher virulence diversity at nip1. Strong evidence for more pathogen migration from barley grass to barley suggests transmission of virulence from barley grass to barley is common.

Conclusions

Pathogen census population size is a better predictor for neutral genetic diversity than host diversity. Despite maintaining a smaller pathogen census population size, barley grass acts as an important ancillary host to R. commune, harbouring highly virulent pathogen types capable of transmission to barley. Management of disease on crops must therefore include management of weedy ancillary hosts, which may harbour disproportionate supplies of virulent pathogen strains.

Electronic supplementary material

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

Filamentous pathogen effectors interfering with small RNA silencing in plant hosts

Filamentous eukaryotic pathogens including fungi and oomycetes are major threats of plant health. During the co-evolutionary arms race with the hosts, these pathogens have evolved a large repertoire of secreted virulence proteins, called effectors, to facilitate colonization and infection. Many effectors are believed to directly manipulate targeted processes inside the host cells; and a fundamental function of the effectors is to dampen immunity. Recent evidence suggests that the destructive oomycete pathogens in the genus Phytophthora encode RNA silencing suppressors. These effectors play an important virulence role during infection, likely through their inhibitory effect on host small RNA-mediated defense.

Identification of protein secretion systems in bacterial genomes

Bacteria with two cell membranes (diderms) have evolved complex systems for protein secretion. These systems were extensively studied in some model bacteria, but the characterisation of their diversity has lagged behind due to lack of standard annotation tools. We built online and standalone computational tools to accurately predict protein secretion systems and related appendages in bacteria with LPS-containing outer membranes. They consist of models describing the systems’ components and genetic organization to be used with MacSyFinder to search for T1SS-T6SS, T9SS, flagella, Type IV pili and Tad pili. We identified ~10,000 candidate systems in bacterial genomes, where T1SS and T5SS were by far the most abundant and widespread. All these data are made available in a public database. The recently described T6SSⁱⁱⁱ and T9SS were restricted to Bacteroidetes, and T6SSⁱⁱ to Francisella. The T2SS, T3SS, and T4SS were frequently encoded in single-copy in one locus, whereas most T1SS were encoded in two loci. The secretion systems of diderm Firmicutes were similar to those found in other diderms. Novel systems may remain to be discovered, since some clades of environmental bacteria lacked all known protein secretion systems. Our models can be fully customized, which should facilitate the identification of novel systems.

Evolution of the EKA family of powdery mildew avirulence-effector genes from the ORF 1 of a LINE retrotransposon

Background

The Avrk1 and Avra10 avirulence (AVR) genes encode effectors that increase the pathogenicity of the fungus Blumeria graminis f.sp. hordei (Bgh), the powdery mildew pathogen, in susceptible barley plants. In resistant barley, MLK1 and MLA10 resistance proteins recognize the presence of AVRK1 and AVRA10, eliciting the hypersensitive response typical of gene for gene interactions. Avrk1 and Avra10 have more than 1350 homologues in Bgh genome, forming the EKA (Effectors homologous to Avrk1 and Avra10) gene family.

Results

We tested the hypothesis that the EKA family originated from degenerate copies of Class I LINE retrotransposons by analysing the EKA family in the genome of Bgh isolate DH14 with bioinformatic tools specially developed for the analysis of Transposable Elements (TE) in genomes. The Class I LINE retrotransposon copies homologous to Avrk1 and Avra10 represent 6.5 % of the Bgh annotated genome and, among them, we identified 293 AVR/effector candidate genes. We also experimentally identified peptides that indicated the translation of several predicted proteins from EKA family members, which had higher relative abundance in haustoria than in hyphae.

Conclusions

Our analyses indicate that Avrk1 and Avra10 have evolved from part of the ORF1 gene of Class I LINE retrotransposons. The co-option of Avra10 and Avrk1 as effectors from truncated copies of retrotransposons explains the huge number of homologues in Bgh genome that could act as dynamic reservoirs from which new effector genes may evolve. These data provide further evidence for recruitment of retrotransposons in the evolution of new biological functions.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2185-x) contains supplementary material, which is available to authorized users.

The fungal core effector Pep1 is conserved across smuts of dicots and monocots

The secreted fungal effector Pep1 is essential for penetration of the host epidermis and establishment of biotrophy in the Ustilago maydis-maize pathosystem. Previously, Pep1 was found to be an inhibitor of apoplastic plant peroxidases, which suppresses the oxidative burst, a primary immune response of the host plant and enables fungal colonization. To investigate the conservation of Pep1 in other pathogens, genomes of related smut species were screened for pep1 orthologues. Pep1 proteins were produced in Escherichia coli for functional assays. The biological function of Pep1 was tested by heterologous expression in U. maydis and Hordeum vulgare. Pep1 orthologues revealed a remarkable degree of sequence conservation, indicating that this effector might play a fundamental role in virulence of biotrophic smut fungi. Pep1 function and its role in virulence are conserved in different pathogenic fungi, even across the monocot-dicot border of host plants. The findings described in this study classify Pep1 as a phylogenetically conserved fungal core effector. Furthermore, we documented the influence of Pep1 on the disease caused by Blumeria graminis f. sp. hordei which is a non-smut-related pathosystem.

Adaptive genomic structural variation in the grape powdery mildew pathogen, Erysiphe necator

BACKGROUND

Powdery mildew, caused by the obligate biotrophic fungus Erysiphe necator, is an economically important disease of grapevines worldwide. Large quantities of fungicides are used for its control, accelerating the incidence of fungicide-resistance. Copy number variations (CNVs) are unbalanced changes in the structure of the genome that have been associated with complex traits. In addition to providing the first description of the large and highly repetitive genome of E. necator, this study describes the impact of genomic structural variation on fungicide resistance in Erysiphe necator.

RESULTS

A shotgun approach was applied to sequence and assemble the genome of five E. necator isolates, and RNA-seq and comparative genomics were used to predict and annotate protein-coding genes. Our results show that the E. necator genome is exceptionally large and repetitive and suggest that transposable elements are responsible for genome expansion. Frequent structural variations were found between isolates and included copy number variation in EnCYP51, the target of the commonly used sterol demethylase inhibitor (DMI) fungicides. A panel of 89 additional E. necator isolates collected from diverse vineyard sites was screened for copy number variation in the EnCYP51 gene and for presence/absence of a point mutation (Y136F) known to result in higher fungicide tolerance. We show that an increase in EnCYP51 copy number is significantly more likely to be detected in isolates collected from fungicide-treated vineyards. Increased EnCYP51 copy numbers were detected with the Y136F allele, suggesting that an increase in copy number becomes advantageous only after the fungicide-tolerant allele is acquired. We also show that EnCYP51 copy number influences expression in a gene-dose dependent manner and correlates with fungal growth in the presence of a DMI fungicide.

CONCLUSIONS

Taken together our results show that CNV can be adaptive in the development of resistance to fungicides by providing increasing quantitative protection in a gene-dosage dependent manner. The results of this work not only demonstrate the effectiveness of using genomics to dissect complex traits in organisms with very limited molecular information, but also may have broader implications for understanding genomic dynamics in response to strong selective pressure in other pathogens with similar genome architectures.

Phytophthora suppressor of RNA silencing 2 is a conserved RxLR effector that promotes infection in soybean and Arabidopsis thaliana

The genus Phytophthora consists of notorious and emerging pathogens of economically important crops. Each Phytophthora genome encodes several hundreds of cytoplasmic effectors, which are believed to manipulate plant immune response inside the host cells. However, the majority of Phytophthora effectors remain functionally uncharacterized. We recently discovered two effectors from the soybean stem and root rot pathogen Phytophthora sojae with the activity to suppress RNA silencing in plants. These effectors are designated Phytophthora suppressor of RNA silencing (PSRs). Here, we report that the P. sojae PSR2 (PsPSR2) belongs to a conserved and widespread effector family in Phytophthora. A PsPSR2-like effector produced by P. infestans (PiPSR2) can also suppress RNA silencing in plants and promote Phytophthora infection, suggesting that the PSR2 family effectors have conserved functions in plant hosts. Using Agrobacterium rhizogenes-mediated hairy roots induction, we demonstrated that the expression of PsPSR2 rendered hypersusceptibility of soybean to P. sojae. Enhanced susceptibility was also observed in PsPSR2-expressing Arabidopsis thaliana plants during Phytophthora but not bacterial infection. These experiments provide strong evidence that PSR2 is a conserved Phytophthora effector family that performs important virulence functions specifically during Phytophthora infection of various plant hosts.

Positive selection and intragenic recombination contribute to high allelic diversity in effector genes of Mycosphaerella fijiensis, causal agent of the black leaf streak disease of banana

Previously, we have determined the nonhost-mediated recognition of the MfAvr4 and MfEcp2 effector proteins from the banana pathogen Mycosphaerella fijiensis in tomato, by the cognate Cf-4 and Cf-Ecp2 resistance proteins, respectively. These two resistance proteins could thus mediate resistance against M. fijiensis if genetically transformed into banana (Musa spp.). However, disease resistance controlled by single dominant genes can be overcome by mutated effector alleles, whose products are not recognized by the cognate resistance proteins. Here, we surveyed the allelic variation within the MfAvr4, MfEcp2, MfEcp2-2 and MfEcp2-3 effector genes of M. fijiensis in a global population of the pathogen, and assayed its impact on recognition by the tomato Cf-4 and Cf-Ecp2 resistance proteins, respectively. We identified a large number of polymorphisms that could reflect a co-evolutionary arms race between host and pathogen. The analysis of nucleotide substitution patterns suggests that both positive selection and intragenic recombination have shaped the evolution of M. fijiensis effectors. Clear differences in allelic diversity were observed between strains originating from South-East Asia relative to strains from other banana-producing continents, consistent with the hypothesis that M. fijiensis originated in the Asian-Pacific region. Furthermore, transient co-expression of the MfAvr4 effector alleles and the tomato Cf-4 resistance gene, as well as of MfEcp2, MfEcp2-2 and MfEcp2-3 and the putative Cf-Ecp2 resistance gene, indicated that effector alleles able to overcome these resistance genes are already present in natural populations of the pathogen, thus questioning the durability of resistance that can be provided by these genes in the field.

Diversifying selection in the wheat stem rust fungus acts predominantly on pathogen-associated gene families and reveals candidate effectors

Plant pathogens cause severe losses to crop plants and threaten global food production. One striking example is the wheat stem rust fungus, Puccinia graminis f. sp. tritici, which can rapidly evolve new virulent pathotypes in response to resistant host lines. Like several other filamentous fungal and oomycete plant pathogens, its genome features expanded gene families that have been implicated in host-pathogen interactions, possibly encoding effector proteins that interact directly with target host defense proteins. Previous efforts to understand virulence largely relied on the prediction of secreted, small and cysteine-rich proteins as candidate effectors and thus delivered an overwhelming number of candidates. Here, we implement an alternative analysis strategy that uses the signal of adaptive evolution as a line of evidence for effector function, combined with comparative information and expression data. We demonstrate that in planta up-regulated genes that are rapidly evolving are found almost exclusively in pathogen-associated gene families, affirming the impact of host-pathogen co-evolution on genome structure and the adaptive diversification of specialized gene families. In particular, we predict 42 effector candidates that are conserved only across pathogens, induced during infection and rapidly evolving. One of our top candidates has recently been shown to induce genotype-specific hypersensitive cell death in wheat. This shows that comparative genomics incorporating the evolutionary signal of adaptation is powerful for predicting effector candidates for laboratory verification. Our system can be applied to a wide range of pathogens and will give insight into host-pathogen dynamics, ultimately leading to progress in strategies for disease control.

Effective identification of Gram-negative bacterial type III secreted effectors using position-specific residue conservation profiles

BACKGROUND

Type III secretion systems (T3SSs) are central to the pathogenesis and specifically deliver their secreted substrates (type III secreted proteins, T3SPs) into host cells. Since T3SPs play a crucial role in pathogen-host interactions, identifying them is crucial to our understanding of the pathogenic mechanisms of T3SSs. This study reports a novel and effective method for identifying the distinctive residues which are conserved different from other SPs for T3SPs prediction. Moreover, the importance of several sequence features was evaluated and further, a promising prediction model was constructed.

RESULTS

Based on the conservation profiles constructed by a position-specific scoring matrix (PSSM), 52 distinctive residues were identified. To our knowledge, this is the first attempt to identify the distinct residues of T3SPs. Of the 52 distinct residues, the first 30 amino acid residues are all included, which is consistent with previous studies reporting that the secretion signal generally occurs within the first 30 residue positions. However, the remaining 22 positions span residues 30-100 were also proven by our method to contain important signal information for T3SP secretion because the translocation of many effectors also depends on the chaperone-binding residues that follow the secretion signal. For further feature optimisation and compression, permutation importance analysis was conducted to select 62 optimal sequence features. A prediction model across 16 species was developed using random forest to classify T3SPs and non-T3 SPs, with high receiver operating curve of 0.93 in the 10-fold cross validation and an accuracy of 94.29% for the test set. Moreover, when performing on a common independent dataset, the results demonstrate that our method outperforms all the others published to date. Finally, the novel, experimentally confirmed T3 effectors were used to further demonstrate the model's correct application. The model and all data used in this paper are freely available at http://cic.scu.edu.cn/bioinformatics/T3SPs.zip.

Global diversity and distribution of three necrotrophic effectors in Phaeosphaeria nodorum and related species

Population genetic and phylogenetic studies have shown that Phaeosphaeria nodorum is a member of a species complex that probably shares its center of origin with wheat (Triticum aestivum and Triticum durum). We examined the evolutionary histories of three known necrotrophic effectors (NEs) produced by P. nodorum and compared them with neutral loci. We screened over 1000 individuals for the presence/absence of each effector and assigned each individual to a multi-effector genotype. Diversity at each NE locus was assessed by sequencing c. 200 individuals for each locus. We found significant differences in effector frequency among populations. We propose that these differences reflect the presence/absence of the corresponding susceptibility gene in wheat cultivars. The population harboring the highest sequence diversity was different for each effector locus and never coincided with populations harboring the highest diversity at neutral loci. Coalescent and phylogenetic analyses showed a discontinuous presence of all three NEs among nine closely related Phaeosphaeria species. Only two of the nine species were found to harbor NEs. We present evidence that the three described NEs of P. nodorum were transmitted to its sister species, Phaeosphaeria avenaria tritici 1, via interspecific hybridization.

Genomic tillage and the harvest of fungal phytopathogens

Genome sequencing has been carried out on a small selection of major fungal ascomycete pathogens. These studies show that simple models whereby pathogens evolved from phylogenetically related saprobes by the acquisition or modification of a small number of key genes cannot be sustained.The genomes show that pathogens cannot be divided into three clearly delineated classes (biotrophs, hemibiotrophs and necrotrophs) but rather into a complex matrix of categories each with subtly different properties. It is clear that the evolution of pathogenicity is ancient, rapid and ongoing. Fungal pathogens have undergone substantial genomic rearrangements that can be appropriately described as 'genomic tillage'. Genomic tillage underpins the evolution and expression of large families of genes - known as effectors - that manipulate and exploit metabolic and defence processes of plants so as to allow the proliferation of pathogens.

Rhynchosporium commune: a persistent threat to barley cultivation

Rhynchosporium commune is a haploid fungus causing scald or leaf blotch on barley, other Hordeum spp. and Bromus diandrus.

TAXONOMY

Rhynchosporium commune is an anamorphic Ascomycete closely related to the teleomorph Helotiales genera Oculimacula and Pyrenopeziza.

DISEASE SYMPTOMS

Rhynchosporium commune causes scald-like lesions on leaves, leaf sheaths and ears. Early symptoms are generally pale grey oval lesions. With time, the lesions acquire a dark brown margin with the centre of the lesion remaining pale green or pale brown. Lesions often merge to form large areas around which leaf yellowing is common. Infection frequently occurs in the leaf axil, which can lead to chlorosis and eventual death of the leaf.

LIFE CYCLE

Rhynchosporium commune is seed borne, but the importance of this phase of the disease is not fully understood. Debris from previous crops and volunteers, infected from the stubble from previous crops, are considered to be the most important sources of the disease. Autumn-sown crops can become infected very soon after sowing. Secondary spread of disease occurs mainly through splash dispersal of conidia from infected leaves. Rainfall at the stem extension growth stage is the major environmental factor in epidemic development. DETECTION AND QUANTIFICATION: Rhynchosporium commune produces unique beak-shaped, one-septate spores both on leaves and in culture. The development of a specific polymerase chain reaction (PCR) and, more recently, quantitative PCR (qPCR) has allowed the identification of asymptomatic infection in seeds and during the growing season.

DISEASE CONTROL

The main measure for the control of R. commune is the use of fungicides with different modes of action, in combination with the use of resistant cultivars. However, this is constantly under review because of the ability of the pathogen to adapt to host plant resistance and to develop fungicide resistance.

Rapid evolution of the sequences and gene repertoires of secreted proteins in bacteria

Proteins secreted to the extracellular environment or to the periphery of the cell envelope, the secretome, play essential roles in foraging, antagonistic and mutualistic interactions. We hypothesize that arms races, genetic conflicts and varying selective pressures should lead to the rapid change of sequences and gene repertoires of the secretome. The analysis of 42 bacterial pan-genomes shows that secreted, and especially extracellular proteins, are predominantly encoded in the accessory genome, i.e. among genes not ubiquitous within the clade. Genes encoding outer membrane proteins might engage more frequently in intra-chromosomal gene conversion because they are more often in multi-genic families. The gene sequences encoding the secretome evolve faster than the rest of the genome and in particular at non-synonymous positions. Cell wall proteins in Firmicutes evolve particularly fast when compared with outer membrane proteins of Proteobacteria. Virulence factors are over-represented in the secretome, notably in outer membrane proteins, but cell localization explains more of the variance in substitution rates and gene repertoires than sequence homology to known virulence factors. Accordingly, the repertoires and sequences of the genes encoding the secretome change fast in the clades of obligatory and facultative pathogens and also in the clades of mutualists and free-living bacteria. Our study shows that cell localization shapes genome evolution. In agreement with our hypothesis, the repertoires and the sequences of genes encoding secreted proteins evolve fast. The particularly rapid change of extracellular proteins suggests that these public goods are key players in bacterial adaptation.

Distinct effects on diversifying selection by two mechanisms of immunity against Streptococcus pneumoniae

Antigenic variation to evade host immunity has long been assumed to be a driving force of diversifying selection in pathogens. Colonization by Streptococcus pneumoniae, which is central to the organism's transmission and therefore evolution, is limited by two arms of the immune system: antibody- and T cell- mediated immunity. In particular, the effector activity of CD4⁺ T_H17 cell mediated immunity has been shown to act in trans, clearing co-colonizing pneumococci that do not bear the relevant antigen. It is thus unclear whether T_H17 cell immunity allows benefit of antigenic variation and contributes to diversifying selection. Here we show that antigen-specific CD4⁺ T_H17 cell immunity almost equally reduces colonization by both an antigen-positive strain and a co-colonized, antigen-negative strain in a mouse model of pneumococcal carriage, thus potentially minimizing the advantage of escape from this type of immunity. Using a proteomic screening approach, we identified a list of candidate human CD4⁺ T_H17 cell antigens. Using this list and a previously published list of pneumococcal Antibody antigens, we bioinformatically assessed the signals of diversifying selection among the identified antigens compared to non-antigens. We found that Antibody antigen genes were significantly more likely to be under diversifying selection than the T_H17 cell antigen genes, which were indistinguishable from non-antigens. Within the Antibody antigens, epitopes recognized by human antibodies showed stronger evidence of diversifying selection. Taken together, the data suggest that T_H17 cell-mediated immunity, one form of T cell immunity that is important to limit carriage of antigen-positive pneumococcus, favors little diversifying selection in the targeted antigen. The results could provide new insight into pneumococcal vaccine design.

Streptococcus pneumoniae, or pneumococcus, is a leading cause of morbidity and mortality in young children and elderly persons worldwide. Current pneumococcus vaccines target a limited number of clinically important serotypes, while strains with serotypes not targeted by current vaccines are increasing in importance in both carriage and invasive disease. As a result, there has been a substantial interest to develop novel, cost-effective vaccines based on protein antigens from pneumococcus. To this end, it is critical to understand how the human immune system exerts selection pressures on the targeted antigens. Two immune mechanisms targeting pneumococcal protein antigens have been documented, mediated by antibody and T cells, respectively. In this study, we screened for pneumococcal antigens that are commonly recognized by human CD4⁺ T_H17 cells. Using a mouse model of pneumococcal colonization, we demonstrate that T_H17 cell-based immunity almost equally reduces colonization by both an antigen-positive strain and a co-colonizing, antigen-negative strain. Furthermore, we demonstrate that the DNA sequences of T_H17 cell antigens demonstrate no detectable signs of being under selective pressure, unlike pneumococcal antigens known to be strong antibody targets. Thus, one form of the T cell-mediated immunity that is important to limit carriage of antigen-positive pneumococcus favors little diversifying selection in the targeted antigen. These results suggest evolution of escape from T_H17 -based vaccines may be slower than from antibody-based vaccines.

Quantitative variation in effector activity of ToxA isoforms from Stagonospora nodorum and Pyrenophora tritici-repentis

ToxA is a proteinaceous necrotrophic effector produced by Stagonospora nodorum and Pyrenophora tritici-repentis. In this study, all eight mature isoforms of the ToxA protein were purified and compared. Circular dichroism spectra indicated that all isoforms were structurally intact and had indistinguishable secondary structural features. ToxA isoforms were infiltrated into wheat lines that carry the sensitivity gene Tsn1. It was observed that different wheat lines carrying identical Tsn1 alleles varied in sensitivity to ToxA. All ToxA isoforms induced necrosis when introduced into any Tsn1 wheat line but we observed quantitative variation in effector activity, with the least-active version found in isolates of P. tritici-repentis. Pathogen sporulation increased with higher doses of ToxA. The isoforms that induced the most rapid necrosis also induced the most sporulation, indicating that pathogen fitness is affected by differences in ToxA activity. We show that differences in toxin activity encoded by a single gene can contribute to the quantitative inheritance of necrotrophic virulence. Our findings support the hypothesis that the variation at ToxA results from selection that favors increased toxin activity.

The cysteine rich necrotrophic effector SnTox1 produced by Stagonospora nodorum triggers susceptibility of wheat lines harboring Snn1

The wheat pathogen Stagonospora nodorum produces multiple necrotrophic effectors (also called host-selective toxins) that promote disease by interacting with corresponding host sensitivity gene products. SnTox1 was the first necrotrophic effector identified in S. nodorum, and was shown to induce necrosis on wheat lines carrying Snn1. Here, we report the molecular cloning and validation of SnTox1 as well as the preliminary characterization of the mechanism underlying the SnTox1-Snn1 interaction which leads to susceptibility. SnTox1 was identified using bioinformatics tools and verified by heterologous expression in Pichia pastoris. SnTox1 encodes a 117 amino acid protein with the first 17 amino acids predicted as a signal peptide, and strikingly, the mature protein contains 16 cysteine residues, a common feature for some avirulence effectors. The transformation of SnTox1 into an avirulent S. nodorum isolate was sufficient to make the strain pathogenic. Additionally, the deletion of SnTox1 in virulent isolates rendered the SnTox1 mutated strains avirulent on the Snn1 differential wheat line. SnTox1 was present in 85% of a global collection of S. nodorum isolates. We identified a total of 11 protein isoforms and found evidence for strong diversifying selection operating on SnTox1. The SnTox1-Snn1 interaction results in an oxidative burst, DNA laddering, and pathogenesis related (PR) gene expression, all hallmarks of a defense response. In the absence of light, the development of SnTox1-induced necrosis and disease symptoms were completely blocked. By comparing the infection processes of a GFP-tagged avirulent isolate and the same isolate transformed with SnTox1, we conclude that SnTox1 may play a critical role during fungal penetration. This research further demonstrates that necrotrophic fungal pathogens utilize small effector proteins to exploit plant resistance pathways for their colonization, which provides important insights into the molecular basis of the wheat-S. nodorum interaction, an emerging model for necrotrophic pathosystems.

In this manuscript we describe the cloning of SnTox1 from Stagonospora nodorum, the gene encoding the first host selective toxin (SnTox1) identified in this fungus. SnTox1 induces necrosis and promotes disease on wheat lines harboring the Snn1 gene. We verified the function of the SnTox1 gene by expressing it in a yeast culture where the resulting culture filtrate induced necrosis but only on wheat lines that carried a functional Snn1. The SnTox1 gene was also transformed into an avirulent S. nodorum isolate, resulting in an isolate that was virulent on wheat lines harboring Snn1. SnTox1 was also disrupted in virulent S. nodorum isolates resulting in the elimination of disease on Snn1 differential wheat lines. Additionally, we investigated the host response to SnTox1 and S. nodorum strains producing SnTox1 and discovered that several hallmarks of a resistance response were present during the susceptible reaction, showing that the necrotrophic pathogen S. nodorum is likely using SnTox1 to stimulate a host resistance pathway involving Snn1 to induce disease.

Functional diversification of the GALA type III effector family contributes to Ralstonia solanacearum adaptation on different plant hosts

Type III effectors from phytopathogenic bacteria exhibit a high degree of functional redundancy, hampering the evaluation of their precise contribution to pathogenicity. This is illustrated by the GALA type III effectors from Ralstonia solanacearum, which have been shown to be collectively, but not individually, required for disease on Arabidopsis thaliana and tomato. We investigated evolution, redundancy and diversification of this family in order to understand the individual contribution of the GALA effectors to pathogenicity. From sequences available, we reconstructed GALA phylogeny and performed selection studies. We then focused on the GALAs from the reference strain GMI1000 to examine their ability to suppress plant defense responses and contribution to pathogenicity on three different host plants: A. thaliana, tomato (Lycopersicum esculentum) and eggplant (Solanum melongena). The GALA family is well conserved within R. solanacearum species. Patterns of selection detected on some GALA family members, together with experimental results, show that GALAs underwent functional diversification. We conclude that functional divergence of the GALA family likely accounts for its remarkable conservation during R. solanacearum evolution and could contribute to R. solanacearum's adaptation on several host plants.

Bacterial genomes: evolution of pathogenicity

Bacterial pathogens continue to pose a major threat to economically important plant resources. Disease outbreaks can occur through rapid evolution of a pathogen to overcome host defences. The advent of genome sequencing, especially next-generation technologies, has seen a revolution in the study of plant pathogen evolution over the past five years. This review highlights recent developments in understanding bacterial plant pathogen evolution, enabled by genomics and specifically focusing on type III protein effectors. The genotypic changes and mechanisms involved in pathogen evolution are now much better understood. However, there is still much to be learned about the drivers of pathogen evolution, both in terms of plant resistance and bacterial lifestyle.

A systems biology perspective on plant-microbe interactions: biochemical and structural targets of pathogen effectors

Plants have biochemical defences against stresses from predators, parasites and pathogens. In this review we discuss the interaction of plant defences with microbial pathogens such as bacteria, fungi and oomycetes, and viruses. We examine principles of complex dynamic networks that allow identification of network components that are differentially and predictably sensitive to perturbation, thus making them likely effector targets. We relate these principles to recent developments in our understanding of known effector targets in plant-pathogen systems, and propose a systems-level framework for the interpretation and modelling of host-microbe interactions mediated by effectors. We describe this framework briefly, and conclude by discussing useful experimental approaches for populating this framework.

Comprehensive appraisal of the extracellular proteins from a monoderm bacterium: theoretical and empirical exoproteomes of Listeria monocytogenes EGD-e by secretomics

Defined as proteins actively transported via secretion systems, secreted proteins can have radically different subcellular destinations in monoderm (Gram-positive) bacteria. From degradative enzymes in saprophytes to virulence factors in pathogens, secreted proteins are the main tools used by bacteria to interact with their surroundings. The etiological agent of listeriosis, Listeria monocytogenes, is a Gram-positive facultative intracellular foodborne pathogen, whose ecological niche is the soil and as such should be primarily considered as a ubiquitous saprophyte. Recent advances on protein secretion systems in this species prompted us to investigate the exoproteome. First, an original and rational bioinformatic strategy was developed to mimic the protein exportation steps leading to the extracellular localization of secreted proteins; 79 exoproteins were predicted as secreted via Sec, 1 exoprotein via Tat, 4 bacteriocins via ABC exporters, 3 exoproteins via holins, and 3 exoproteins via the WXG100 system. This bioinformatic analysis allowed for defining a databank of the mature protein set in L. monocytogenes, which was used for generating the theoretical exoproteome and for subsequent protein identification by proteomics. 2-DE proteomic analyses were performed over a wide pI range to experimentally cover the largest protein spectrum possible. A total of 120 spots could be resolved and identified, which corresponded to 50 distinct proteins. These exoproteins were essentially virulence factors, degradative enzymes, and proteins of unknown functions, which exportation would essentially rely on the Sec pathway or nonclassical secretion. This investigation resulted in the first comprehensive appraisal of the exoproteome of L. monocytogenes EGD-e based on theoretical and experimental secretomic analyses, which further provided indications on listerial physiology in relation with its habitat and lifestyle. The novel and rational strategy described here is generic and has been purposely designed for the prediction of proteins localized extracellularly in monoderm bacteria.

T3SEdb: data warehousing of virulence effectors secreted by the bacterial Type III Secretion System

Background

Effectors of Type III Secretion System (T3SS) play a pivotal role in establishing and maintaining pathogenicity in the host and therefore the identification of these effectors is important in understanding virulence. However, the effectors display high level of sequence diversity, therefore making the identification a difficult process. There is a need to collate and annotate existing effector sequences in public databases to enable systematic analyses of these sequences for development of models for screening and selection of putative novel effectors from bacterial genomes that can be validated by a smaller number of key experiments.

Results

Herein, we present T3SEdb http://effectors.bic.nus.edu.sg/T3SEdb, a specialized database of annotated T3SS effector (T3SE) sequences containing 1089 records from 46 bacterial species compiled from the literature and public protein databases. Procedures have been defined for i) comprehensive annotation of experimental status of effectors, ii) submission and curation review of records by users of the database, and iii) the regular update of T3SEdb existing and new records. Keyword fielded and sequence searches (BLAST, regular expression) are supported for both experimentally verified and hypothetical T3SEs. More than 171 clusters of T3SEs were detected based on sequence identity comparisons (intra-cluster difference up to ~60%). Owing to this high level of sequence diversity of T3SEs, the T3SEdb provides a large number of experimentally known effector sequences with wide species representation for creation of effector predictors. We created a reliable effector prediction tool, integrated into the database, to demonstrate the application of the database for such endeavours.

Conclusions

T3SEdb is the first specialised database reported for T3SS effectors, enriched with manual annotations that facilitated systematic construction of a reliable prediction model for identification of novel effectors. The T3SEdb represents a platform for inclusion of additional annotations of metadata for future developments of sophisticated effector prediction models for screening and selection of putative novel effectors from bacterial genomes/proteomes that can be validated by a small number of key experiments.

The bacterial pathogen Xanthomonas oryzae overcomes rice defenses by regulating host copper redistribution

Pathogen effectors are virulence factors causing plant diseases. How the host targets of these effectors facilitate pathogen infection is largely unknown. An effector of Xanthomonas oryzae pv oryzae (Xoo) transcriptionally activates rice (Oryza sativa) susceptibility gene Xa13 to cause bacterial blight disease. Xa13 encodes an indispensable plasma membrane protein of the MtN3/saliva family, which is prevalent in eukaryotes with unknown biochemical function. We show that the XA13 protein cooperates with two other proteins, COPT1 and COPT5, to promote removal of copper from xylem vessels, where Xoo multiplies and spreads to cause disease. Copper, an essential micronutrient of plants and an important element for a number of pesticides in agriculture, suppresses Xoo growth. Xoo strain PXO99 is more sensitive to copper than other strains; its infection of rice is associated with activation of XA13, COPT1, and COPT5, which modulate copper redistribution in rice. The involvement of XA13 in copper redistribution has led us to propose a mechanism of bacterial virulence.

S-glycoprotein-like protein regulates defense responses in Nicotiana plants against Ralstonia solanacearum

RsRGA4 (for Ralstonia solanacearum-responsive gene A4) encodes a polypeptide similar to S-locus glycoprotein (SGP) from Brassica rapa and SGP-like proteins from Ipomoea trifida and Medicago truncatula. Therefore, we designated RsRGA4 as NtSGLP (for Nicotiana tabacum SGP-like protein) and NbSGLP (its Nicotiana benthamiana ortholog). NbSGLP is expressed in root, leaf, petal, gynoecium, and stamen. Expression of NbSGLP was strongly induced by inoculation with an avirulent strain of R. solanacearum (Rs8107) and slightly enhanced by inoculation with virulent R. solanacearum (RsOE1-1). Expression of NbSGLP was induced by inoculation with an hrpY-deficient mutant of RsOE1-1 and Rs8107. Expression was also induced by aminocyclopropane carboxylic acid and salicylic acid. Virus-induced gene silencing of NbSGLP enhanced the growth of Rs8107. Growth of RsOE1-1 and appearance of wilt symptoms were also accelerated in silenced plants. Expression of PR-1a and EREBP was reduced, and markers for basal defense, such as callose deposition and reduced vascular flow, were compromised in NbSGLP-silenced plants. Moreover, growth of Pseudomonas cichorii, Pseudomonas syringae pv tabaci, and P. syringae pv mellea was also enhanced in the silenced plants. On the other hand, silencing of NbSGLP did not interfere with the appearance of the hypersensitive response. NbSGLP was secreted in a signal peptide-dependent manner. Agrobacterium tumefaciens-mediated expression of NbSGLP induced PR-1a and EREBP expression, callose deposition, and reduced vascular flow. NbSGLP-induced callose deposition and reduced vascular flow were not observed in salicylic acid-deficient N. benthamiana NahG plants. Taken together, SGLP might have a role in the induction of basal defense in Nicotiana plants.

Catalytic domain of the diversified Pseudomonas syringae type III effector HopZ1 determines the allelic specificity in plant hosts

The type III secretion systems (T3SS) and secreted effectors (T3SEs) are essential virulence factors in Gram-negative bacteria. During the arms race, plants have evolved resistance (R) genes to detect specific T3SEs and activate defence responses. However, this immunity can be efficiently defeated by the pathogens through effector evolution. HopZ1 of the plant pathogen Pseudomonas syringae is a member of the widely distributed YopJ T3SE family. Three alleles are known to be present in P. syringae, with HopZ1a most resembling the ancestral allelic form. In this study, molecular mechanisms underlying the sequence diversification-enabled HopZ1 allelic specificity is investigated. Using domain shuffling experiments, we present evidence showing that a central domain upstream of the conserved catalytic cysteine residue determines HopZ1 recognition specificity. Random and targeted mutagenesis identified three amino acids involved in HopZ1 allelic specificity. Particularly, the exchange of cysteine141 in HopZ1a with lysine137 at the corresponding position in HopZ1b abolished HopZ1a recognition in soybean. This position is under strong positive selection, suggesting that the cysteine/lysine mutation might be a key step driving the evolution of HopZ1. Our data support a model in which sequence diversification imposed by the plant R gene-associated immunity has driven HopZ1 evolution by allowing allele-specific substrate-binding.

The type III effector HopF2Pto targets Arabidopsis RIN4 protein to promote Pseudomonas syringae virulence

Plant immunity can be induced by two major classes of pathogen-associated molecules. Pathogen- or microbe-associated molecular patterns (PAMPs or MAMPs) are conserved molecular components of microbes that serve as "non-self" features to induce PAMP-triggered immunity (PTI). Pathogen effector proteins used to promote virulence can also be recognized as "non-self" features or induce a "modified-self" state that can induce effector-triggered immunity (ETI). The Arabidopsis protein RIN4 plays an important role in both branches of plant immunity. Three unrelated type III secretion effector (TTSE) proteins from the phytopathogen Pseudomonas syringae, AvrRpm1, AvrRpt2, and AvrB, target RIN4, resulting in ETI that effectively restricts pathogen growth. However, no pathogenic advantage has been demonstrated for RIN4 manipulation by these TTSEs. Here, we show that the TTSE HopF2(Pto) also targets Arabidopsis RIN4. Transgenic plants conditionally expressing HopF2(Pto) were compromised for AvrRpt2-induced RIN4 modification and associated ETI. HopF2(Pto) interfered with AvrRpt2-induced RIN4 modification in vitro but not with AvrRpt2 activation, suggestive of RIN4 targeting by HopF2(Pto). In support of this hypothesis, HopF2 (Pto) interacted with RIN4 in vitro and in vivo. Unlike AvrRpm1, AvrRpt2, and AvrB, HopF2(Pto) did not induce ETI and instead promoted P. syringae growth in Arabidopsis. This virulence activity was not observed in plants genetically lacking RIN4. These data provide evidence that RIN4 is a major virulence target of HopF2(Pto) and that a pathogenic advantage can be conveyed by TTSEs that target RIN4.

Computational prediction of type III secreted proteins from gram-negative bacteria

Background

Type III secretion system (T3SS) is a specialized protein delivery system in gram-negative bacteria that injects proteins (called effectors) directly into the eukaryotic host cytosol and facilitates bacterial infection. For many plant and animal pathogens, T3SS is indispensable for disease development. Recently, T3SS has also been found in rhizobia and plays a crucial role in the nodulation process. Although a great deal of efforts have been done to understand type III secretion, the precise mechanism underlying the secretion and translocation process has not been fully understood. In particular, defined secretion and translocation signals enabling the secretion have not been identified from the type III secreted effectors (T3SEs), which makes the identification of these important virulence factors notoriously challenging. The availability of a large number of sequenced genomes for plant and animal-associated bacteria demands the development of efficient and effective prediction methods for the identification of T3SEs using bioinformatics approaches.

Results

We have developed a machine learning method based on the N-terminal amino acid sequences to predict novel type III effectors in the plant pathogen Pseudomonas syringae and the microsymbiont rhizobia. The extracted features used in the learning model (or classifier) include amino acid composition, secondary structure and solvent accessibility information. The method achieved a precision of over 90% on P. syringae in a cross validation study. In combination with a promoter screen for the type III specific promoters, this classifier trained on the P. syringae data was applied to predict novel T3SEs from the genomic sequences of four rhizobial strains. This application resulted in 57 candidate type III secreted proteins, 17 of which are confirmed effectors.

Conclusion

Our experimental results demonstrate that the machine learning method based on N-terminal amino acid sequences combined with a promoter screen could prove to be a very effective computational approach for predicting novel type III effectors in gram-negative bacteria. Our method and data are available to the public upon request.

Positive selection identifies an in vivo role for FimH during urinary tract infection in addition to mannose binding

FimH, the type 1 pilus adhesin of uropathogenic Escherichia coli (UPEC), contains a receptor-binding domain with an acidic binding pocket specific for mannose. The fim operon, and thus type 1 pilus production, is under transcriptional control via phase variation of an invertible promoter element. FimH is critical during urinary tract infection for mediating colonization and invasion of the bladder epithelium and establishment of intracellular bacterial communities (IBCs). In silico analysis of FimH gene sequences from 279 E. coli strains identified specific amino acids evolving under positive selection outside of its mannose-binding pocket. Mutating two of these residues (A27V/V163A) had no effect on phase variation, pilus assembly, or mannose binding in vitro. However, compared to wild-type, this double mutant strain exhibited a 10,000-fold reduction in mouse bladder colonization 24 h after inoculation and was unable to form IBCs even though it bound normally to mannosylated receptors in the urothelium. In contrast, the single A62S mutation altered phase variation, reducing the proportion of piliated cells, reduced mannose binding 8-fold, and decreased bladder colonization 30-fold in vivo compared to wild-type. A phase-locked ON A62S mutant restored virulence to wild-type levels even though in vitro mannose binding remained impaired. Thus, positive selection analysis of FimH has separated mannose binding from in vivo fitness, suggesting that IBC formation is critical for successful infection of the mammalian bladder, providing support for more general use of in silico positive selection analysis to define the molecular underpinnings of bacterial pathogenesis.

Roadmap for future research on plant pathogen effectors

Bacterial and eukaryotic plant pathogens deliver effector proteins into plant cells to promote pathogenesis. Bacterial pathogens containing type III protein secretion systems are known to inject many of these effectors into plant cells. More recently, oomycete pathogens have been shown to possess a large family of effectors containing the RXLR motif, and many effectors are also being discovered in fungal pathogens. Although effector activities are largely unknown, at least a subset suppress plant immunity. A plethora of new plant pathogen genomes that will soon be available thanks to next-generation sequencing technologies will allow the identification of many more effectors. This article summarizes the key approaches used to identify plant pathogen effectors, many of which will continue to be useful for future effector discovery. Thus, it can be viewed as a 'roadmap' for effector and effector target identification. Because effectors can be used as tools to elucidate components of innate immunity, advances in our understanding of effectors and their targets should lead to improvements in agriculture.

The evolution of Pseudomonas syringae host specificity and type III effector repertoires

The discovery 45 years ago that many Pseudomonas syringae pathovars elicit the hypersensitive response in plant species other than their hosts fostered the use of these bacteria as experimental models. However, the basis for host specificity and the corresponding resistance of nonhosts remain unclear. Pseudomonas syringae is now known to inject into the host cytoplasm, via the type III secretion system, effector proteins that suppress basal innate immunity, but may be recognized by cognate resistance (R) proteins in a second level of defence. The identification and manipulation of complete repertoires of type III effectors have revealed the highly polymorphic nature of effector repertoires and their potential to limit the host range. However, the maintenance of compatible effector repertoires may be driven by adaptations to life in a given plant species involving many factors. Tools are now available to test several hypotheses for the nature and evolution of P. syringae host specificity and nonhost resistance.

Coevolution between a family of parasite virulence effectors and a class of LINE-1 retrotransposons

Parasites are able to evolve rapidly and overcome host defense mechanisms, but the molecular basis of this adaptation is poorly understood. Powdery mildew fungi (Erysiphales, Ascomycota) are obligate biotrophic parasites infecting nearly 10,000 plant genera. They obtain their nutrients from host plants through specialized feeding structures known as haustoria. We previously identified the AVR(k1) powdery mildew-specific gene family encoding effectors that contribute to the successful establishment of haustoria. Here, we report the extensive proliferation of the AVR(k1) gene family throughout the genome of B. graminis, with sequences diverging in formae speciales adapted to infect different hosts. Also, importantly, we have discovered that the effectors have coevolved with a particular family of LINE-1 retrotransposons, named TE1a. The coevolution of these two entities indicates a mutual benefit to the association, which could ultimately contribute to parasite adaptation and success. We propose that the association would benefit 1) the powdery mildew fungus, by providing a mechanism for amplifying and diversifying effectors and 2) the associated retrotransposons, by providing a basis for their maintenance through selection in the fungal genome.