Population-genomic insights into emergence, crop-adaptation, and dissemination of Pseudomonas syringae pathogens

This is the author accepted manuscript. The final version is available from the publisher via the DOI in this record.Many bacterial pathogens are well characterized but, in some cases, relatively little is known about the populations from which they emerged. This limits understanding of the molecular mechanisms underlying disease. The crop pathogen Pseudomonas syringae sensu lato has been widely isolated from the environment, including wild plants and components of the water cycle, and causes disease in several economically important crops. Here, we compared genome sequences of 45 P. syringae crop pathogen outbreak strains with 69 closely related environmental isolates. Phylogenetic reconstruction revealed that crop pathogens emerged many times independently from environmental populations. Unexpectedly, differences in gene content between environmental populations and outbreak strains were minimal with most virulence genes present in both. However, a genome-wide association study identified a small number of genes, including the type III effector genes hopQ1 and hopD1, to be associated with crop pathogens, but not with environmental populations, suggesting that this small group of genes may play an important role in crop disease emergence. Intriguingly, genome-wide analysis of homologous recombination revealed that the locus Psyr 0346, predicted to encode a protein that confers antibiotic resistance, has been frequently exchanged among lineages and thus may contribute to pathogen fitness. Finally, we found that isolates from diseased crops and from components of the water cycle, collected during the same crop disease epidemic, form a single population. This provides the strongest evidence yet that precipitation and irrigation water are an overlooked inoculum source for disease epidemics caused by P. syringae.Caroline L. Monteil received support from INRA and the European Union, in the framework of the Marie-Curie FP7 COFUND People Programme, through the award of an AgreenSkills’ fellowship (under grant agreement n° 267196). Research in Boris A. Vinatzer’s laboratory and genome sequencing was funded by the National Science Foundation of the USA (grants IOS-1354215 and DEB-1241068). Funding for work in the Vinatzer laboratory was also provided in part by the Virginia Agricultural Experiment Station and the Hatch Program of the National Institute of Food and Agriculture, U.S. Department of Agriculture. Work carried out in the Sheppard laboratory was supported by the Medical Research Council (MRC) grant MR/L015080/1, and the Wellcome Trust grant 088786/C/09/Z. GM was supported by a NISCHR Health Research Fellowship (HF-14-13)

The plant pathogen Pseudomonas syringae pv. tomato is genetically monomorphic and under strong selection to evade tomato immunity

addresses: Department of Plant Pathology, Physiology, and Weed Science, Virginia Tech, Blacksburg, Virginia, United States of America.notes: PMCID: PMC3161960types: Journal Article; Research Support, U.S. Gov't, Non-P.H.S.This is an open-access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication.Recently, genome sequencing of many isolates of genetically monomorphic bacterial human pathogens has given new insights into pathogen microevolution and phylogeography. Here, we report a genome-based micro-evolutionary study of a bacterial plant pathogen, Pseudomonas syringae pv. tomato. Only 267 mutations were identified between five sequenced isolates in 3,543,009 nt of analyzed genome sequence, which suggests a recent evolutionary origin of this pathogen. Further analysis with genome-derived markers of 89 world-wide isolates showed that several genotypes exist in North America and in Europe indicating frequent pathogen movement between these world regions. Genome-derived markers and molecular analyses of key pathogen loci important for virulence and motility both suggest ongoing adaptation to the tomato host. A mutational hotspot was found in the type III-secreted effector gene hopM1. These mutations abolish the cell death triggering activity of the full-length protein indicating strong selection for loss of function of this effector, which was previously considered a virulence factor. Two non-synonymous mutations in the flagellin-encoding gene fliC allowed identifying a new microbe associated molecular pattern (MAMP) in a region distinct from the known MAMP flg22. Interestingly, the ancestral allele of this MAMP induces a stronger tomato immune response than the derived alleles. The ancestral allele has largely disappeared from today's Pto populations suggesting that flagellin-triggered immunity limits pathogen fitness even in highly virulent pathogens. An additional non-synonymous mutation was identified in flg22 in South American isolates. Therefore, MAMPs are more variable than expected differing even between otherwise almost identical isolates of the same pathogen strain

Genomic organisation of the Mal d 1 gene cluster on linkage group 16 in apple

European populations exhibit progressive sensitisation to food allergens, and apples are one of the foods for which sensitisation is observed most frequently. Apple cultivars vary greatly in their allergenic characteristics, and a better understanding of the genetic basis of low allergenicity may therefore allow allergic individuals to increase their fruit intake. Mal d 1 is considered to be a major apple allergen, and this protein is encoded by the most complex allergen gene family. Not all Mal d 1 members are likely to be involved in allergenicity. Therefore, additional knowledge about the existence and characteristics of the different Mal d 1 genes is required. In the present study, we investigated the genomic organisation of the Mal d 1 gene cluster in linkage group 16 of apple through the sequencing of two bacterial artificial chromosome clones. The results provided new information on the composition of this family with respect to the number and orientation of functional and pseudogenes and their physical distances. The results were compared with the apple and peach genome sequences that have recently been made available. A broad analysis of the whole apple genome revealed the presence of new genes in this family, and a complete list of the observed Mal d 1 genes is supplied. Thus, this study provides an important contribution towards a better understanding of the genetics of the Mal d 1 family and establishes the basis for further research on allelic diversity among cultivars in relation to variation in allergenicity

Deletions in the Repertoire of Pseudomonas syringae pv. tomato DC3000 Type III Secretion Effector Genes Reveal Functional Overlap among Effectors

The γ-proteobacterial plant pathogen Pseudomonas syringae pv. tomato DC3000 uses the type III secretion system to inject ca. 28 Avr/Hop effector proteins into plants, which enables the bacterium to grow from low inoculum levels to produce bacterial speck symptoms in tomato, Arabidopsis thaliana, and (when lacking hopQ1-1) Nicotiana benthamiana. The effectors are collectively essential but individually dispensable for the ability of the bacteria to defeat defenses, grow, and produce symptoms in plants. Eighteen of the effector genes are clustered in six genomic islands/islets. Combinatorial deletions involving these clusters and two of the remaining effector genes revealed a redundancy-based structure in the effector repertoire, such that some deletions diminished growth in N. benthamiana only in combination with other deletions. Much of the ability of DC3000 to grow in N. benthamiana was found to be due to five effectors in two redundant-effector groups (REGs), which appear to separately target two high-level processes in plant defense: perception of external pathogen signals (AvrPto and AvrPtoB) and deployment of antimicrobial factors (AvrE, HopM1, HopR1). Further support for the membership of HopR1 in the same REG as AvrE was gained through bioinformatic analysis, revealing the existence of an AvrE/DspA/E/HopR effector superfamily, which has representatives in virtually all groups of proteobacterial plant pathogens that deploy type III effectors

Performance and long-term stability of the barley hordothionin gene in multiple transgenic apple lines

Introduction of sustainable scab resistance in elite apple cultivars is of high importance for apple cultivation when aiming at reducing the use of chemical crop protectants. Genetic modification (GM) allows the rapid introduction of resistance genes directly into high quality apple cultivars. Resistance genes can be derived from apple itself but genetic modification also opens up the possibility to use other, non-host resistance genes. A prerequisite for application is the long-term performance and stability of the gene annex trait in the field. For this study, we produced and selected a series of transgenic apple lines of two cultivars, i.e. ‘Elstar’ and ‘Gala’ in which the barley hordothionin gene (hth) was introduced. After multiplication, the GM hth-lines, non-GM susceptible and resistant controls and GM non-hth controls were planted in a random block design in a field trial in 40 replicates. Scab resistance was monitored after artificial inoculation (first year) and after natural infection (subsequent years). After the trial period, the level of expression of the hth gene was checked by quantitative RT-PCR. Four of the six GM hth apple lines proved to be significantly less susceptible to apple scab and this trait was found to be stable for the entire 4-year period. Hth expression at the mRNA level was also stable

Meta-analytic approach to the accurate prediction of secreted virulence effectors in gram-negative bacteria

<p>Abstract</p> <p>Background</p> <p>Many pathogens use a type III secretion system to translocate virulence proteins (called effectors) in order to adapt to the host environment. To date, many prediction tools for effector identification have been developed. However, these tools are insufficiently accurate for producing a list of putative effectors that can be applied directly for labor-intensive experimental verification. This also suggests that important features of effectors have yet to be fully characterized.</p> <p>Results</p> <p>In this study, we have constructed an accurate approach to predicting secreted virulence effectors from Gram-negative bacteria. This consists of a support vector machine-based discriminant analysis followed by a simple criteria-based filtering. The accuracy was assessed by estimating the average number of true positives in the top-20 ranking in the genome-wide screening. In the validation, 10 sets of 20 training and 20 testing examples were randomly selected from 40 known effectors of <it>Salmonella enterica </it>serovar Typhimurium LT2. On average, the SVM portion of our system predicted 9.7 true positives from 20 testing examples in the top-20 of the prediction. Removal of the N-terminal instability, codon adaptation index and ProtParam indices decreased the score to 7.6, 8.9 and 7.9, respectively. These discrimination features suggested that the following characteristics of effectors had been uncovered: unstable N-terminus, non-optimal codon usage, hydrophilic, and less aliphathic. The secondary filtering process represented by coexpression analysis and domain distribution analysis further refined the average true positive counts to 12.3. We further confirmed that our system can correctly predict known effectors of <it>P. syringae </it>DC3000, strongly indicating its feasibility.</p> <p>Conclusions</p> <p>We have successfully developed an accurate prediction system for screening effectors on a genome-wide scale. We confirmed the accuracy of our system by external validation using known effectors of <it>Salmonella </it>and obtained the accurate list of putative effectors of the organism. The level of accuracy was sufficient to yield candidates for gene-directed experimental verification. Furthermore, new features of effectors were revealed: non-optimal codon usage and instability of the N-terminal region. From these findings, a new working hypothesis is proposed regarding mechanisms controlling the translocation of virulence effectors and determining the substrate specificity encoded in the secretion system.</p

Genome Sequence Analyses of Pseudomonas savastanoi pv. glycinea and Subtractive Hybridization-Based Comparative Genomics with Nine Pseudomonads

Bacterial blight, caused by Pseudomonas savastanoi pv. glycinea (Psg), is a common disease of soybean. In an effort to compare a current field isolate with one isolated in the early 1960s, the genomes of two Psg strains, race 4 and B076, were sequenced using 454 pyrosequencing. The genomes of both Psg strains share more than 4,900 highly conserved genes, indicating very low genetic diversity between Psg genomes. Though conserved, genome rearrangements and recombination events occur commonly within the two Psg genomes. When compared to each other, 437 and 163 specific genes were identified in B076 and race 4, respectively. Most specific genes are plasmid-borne, indicating that acquisition and maintenance of plasmids may represent a major mechanism to change the genetic composition of the genome and even acquire new virulence factors. Type three secretion gene clusters of Psg strains are near identical with that of P. savastanoi pv. phaseolicola (Pph) strain 1448A and they shared 20 common effector genes. Furthermore, the coronatine biosynthetic cluster is present on a large plasmid in strain B076, but not in race 4. In silico subtractive hybridization-based comparative genomic analyses with nine sequenced phytopathogenic pseudomonads identified dozens of specific islands (SIs), and revealed that the genomes of Psg strains are more similar to those belonging to the same genomospecies such as Pph 1448A than to other phytopathogenic pseudomonads. The number of highly conserved genes (core genome) among them decreased dramatically when more genomes were included in the subtraction, suggesting the diversification of pseudomonads, and further indicating the genome heterogeneity among pseudomonads. However, the number of specific genes did not change significantly, suggesting these genes are indeed specific in Psg genomes. These results reinforce the idea of a species complex of P. syringae and support the reclassification of P. syringae into different species

The Plant Pathogen Pseudomonas syringae pv. tomato Is Genetically Monomorphic and under Strong Selection to Evade Tomato Immunity

Recently, genome sequencing of many isolates of genetically monomorphic bacterial human pathogens has given new insights into pathogen microevolution and phylogeography. Here, we report a genome-based micro-evolutionary study of a bacterial plant pathogen, Pseudomonas syringae pv. tomato. Only 267 mutations were identified between five sequenced isolates in 3,543,009 nt of analyzed genome sequence, which suggests a recent evolutionary origin of this pathogen. Further analysis with genome-derived markers of 89 world-wide isolates showed that several genotypes exist in North America and in Europe indicating frequent pathogen movement between these world regions. Genome-derived markers and molecular analyses of key pathogen loci important for virulence and motility both suggest ongoing adaptation to the tomato host. A mutational hotspot was found in the type III-secreted effector gene hopM1. These mutations abolish the cell death triggering activity of the full-length protein indicating strong selection for loss of function of this effector, which was previously considered a virulence factor. Two non-synonymous mutations in the flagellin-encoding gene fliC allowed identifying a new microbe associated molecular pattern (MAMP) in a region distinct from the known MAMP flg22. Interestingly, the ancestral allele of this MAMP induces a stronger tomato immune response than the derived alleles. The ancestral allele has largely disappeared from today's Pto populations suggesting that flagellin-triggered immunity limits pathogen fitness even in highly virulent pathogens. An additional non-synonymous mutation was identified in flg22 in South American isolates. Therefore, MAMPs are more variable than expected differing even between otherwise almost identical isolates of the same pathogen strain

Dynamic Evolution of Pathogenicity Revealed by Sequencing and Comparative Genomics of 19 Pseudomonas syringae Isolates

Closely related pathogens may differ dramatically in host range, but the molecular, genetic, and evolutionary basis for these differences remains unclear. In many Gram- negative bacteria, including the phytopathogen Pseudomonas syringae, type III effectors (TTEs) are essential for pathogenicity, instrumental in structuring host range, and exhibit wide diversity between strains. To capture the dynamic nature of virulence gene repertoires across P. syringae, we screened 11 diverse strains for novel TTE families and coupled this nearly saturating screen with the sequencing and assembly of 14 phylogenetically diverse isolates from a broad collection of diseased host plants. TTE repertoires vary dramatically in size and content across all P. syringae clades; surprisingly few TTEs are conserved and present in all strains. Those that are likely provide basal requirements for pathogenicity. We demonstrate that functional divergence within one conserved locus, hopM1, leads to dramatic differences in pathogenicity, and we demonstrate that phylogenetics-informed mutagenesis can be used to identify functionally critical residues of TTEs. The dynamism of the TTE repertoire is mirrored by diversity in pathways affecting the synthesis of secreted phytotoxins, highlighting the likely role of both types of virulence factors in determination of host range. We used these 14 draft genome sequences, plus five additional genome sequences previously reported, to identify the core genome for P. syringae and we compared this core to that of two closely related non-pathogenic pseudomonad species. These data revealed the recent acquisition of a 1 Mb megaplasmid by a sub-clade of cucumber pathogens. This megaplasmid encodes a type IV secretion system and a diverse set of unknown proteins, which dramatically increases both the genomic content of these strains and the pan-genome of the species

Functional analysis and expression profiling of HcrVf1 and HcrVf2 for development of scab resistant cisgenic and intragenic apples

Apple scab resistance genes, HcrVf1 and HcrVf2, were isolated including their native promoter, coding and terminator sequences. Two fragment lengths (short and long) of the native gene promoters and the strong apple rubisco gene promoter (PMdRbc) were used for both HcrVf genes to test their effect on expression and phenotype. The scab susceptible cultivar ‘Gala’ was used for plant transformations and after selection of transformants, they were micrografted onto apple seedling rootstocks for scab disease tests. Apple transformants were also tested for HcrVf expression by quantitative RT-PCR (qRT-PCR). For HcrVf1 the long native promoter gave significantly higher expression that the short one; in case of HcrVf2 the difference between the two was not significant. The apple rubisco gene promoter proved to give the highest expression of both HcrVf1 and HcrVf2. The top four expanding leaves were used initially for inoculation with monoconidial isolate EU-B05 which belongs to race 1 of V. inaequalis. Later six other V. inaequalis isolates were used to study the resistance spectra of the individual HcrVf genes. The scab disease assays showed that HcrVf1 did not give resistance against any of the isolates tested regardless of the expression level. The HcrVf2 gene appeared to be the only functional gene for resistance against Vf avirulent isolates of V. inaequalis. HcrVf2 did not provide any resistance to Vf virulent strains, even not in case of overexpression. In conclusion, transformants carrying the apple-derived HcrVf2 gene in a cisgenic as well as in an intragenic configuration were able to reach scab resistance levels comparable to the Vf resistant control cultivar obtained by classical breeding, cv. ‘Santana’