Genomic structure and phylogeny of the plant pathogen Ralstonia solanacearum inferred from gene distribution analysis

Comparative genomics and host range analysis of four Ralstonia pseudosolanacearum strains isolated from sunflower reveals genomic and phenotypic differences

BACKGROUND

Bacterial wilt caused by Ralstonia solanacearum species complex (RSSC) is one of the devastating diseases in crop production, seriously reducing the yield of crops. R. pseudosolanacearum, is known for its broad infrasubspecific diversity and comprises 36 sequevars that are currently known. Previous studies found that R. pseudosolanacearum contained four sequevars (13, 14, 17 and 54) isolated from sunflowers sown in the same field.

RESULTS

Here, we provided the complete genomes and the results of genome comparison of the four sequevars strains (RS639, RS642, RS647, and RS650). Four strains showed different pathogenicities to the same cultivars and different host ranges. Their genome sizes were about 5.84 ~ 5.94 Mb, encoding 5002 ~ 5079 genes and the average G + C content of 66.85% ~ 67%. Among the coding genes, 146 ~ 159 specific gene families (contained 150 ~ 160 genes) were found in the chromosomes and 34 ~ 77 specific gene families (contained 34 ~ 78 genes) in the megaplasmids from four strains. The average nucleotide identify (ANI) values between any two strains ranged from 99.05% ~ 99.71%, and the proportion of the total base length of collinear blocks accounts for the total gene length of corresponding genome was all more than 93.82%. Then, we performed a search for genomic islands, prophage sequences, the gene clusters macromolecular secretion systems, type III secreted effectors and other virulence factors in these strains, which provided detailed comparison results of their presence and distinctive features compared to the reference strain GMI1000. Among them, the number and types of T2SS gene clusters were different in the four strains, among which RS650 included all five types. T4SS gene cluster of RS639 and RS647 were missed. In the T6SS gene cluster, several genes were inserted in the RS639, RS647, and RS650, and gene deletion was also detected in the RS642. A total of 78 kinds of type III secreted effectors were found, which included 52 core and 9 specific effectors in four strains.

CONCLUSION

This study not only provided the complete genomes of multiple R. pseudosolanacearum strains isolated from a new host, but also revealed the differences in their genomic levels through comparative genomics. Furthermore, these findings expand human knowledge about the range of hosts that Ralstonia can infect, and potentially contribute to exploring rules and factors of the genetic evolution and analyzing its pathogenic mechanism.

Loss of mobile genomic islands in metal resistant, hydrogen-oxidizing Cupriavidus metallidurans

The genome of the metal resistant, hydrogen-oxidizing bacterium Cupriavidus metallidurans strain CH34 contains horizontally acquired plasmids and genomic islands. Metal-resistance determinants on the two plasmids may exert genetic dominance over other related determinants. To investigate whether these recessive determinants can be activated in the absence of the dominant ones, the transcriptome of the highly zinc-sensitive deletion mutant Δe4 (ΔcadA ΔzntA ΔdmeF ΔfieF) of the plasmid-free parent AE104 was characterized using gene arrays. As a consequence of some unexpected results, close examination by PCR and genomic re-resequencing of strains CH34, AE104, Δe4 and others revealed that the genomic islands CMGIs 2, 3, 4, D, E, but no other islands or recessive determinants, were deleted in some of these strains. Provided CH34 wild type was kept under alternating zinc and nickel selection pressure, no comparable deletions occurred. All current data suggest that genes were actually deleted and were not, as previously surmised, simply absent from the respective strain. As a consequence, a cured database was compiled from the newly generated and previously published gene array data. Analysis of data from this database indicated that some genes of recessive, no longer needed determinants were nevertheless expressed and up-regulated. Their products may interact with those of the dominant determinants to mediate a mosaic phenotype. The ability to contribute to such a mosaic phenotype may prevent deletion of the recessive determinant. The data suggest that the bacterium actively modifies its genome to deal with metal stress and the same time ensures metal homeostasis. Significance In their natural environment, bacteria continually acquire genes by horizontal gene transfer and newly acquired determinants may become dominant over related ones already present in the host genome. When a bacterium is taken into laboratory culture, it is isolated from the horizontal gene transfer network. It can no longer gain genes, but instead may lose them. This was indeed observed in Cupriavidus metallidurans for loss key metal-resistance determinants when no selection pressure was continuously kept. However, some recessive metal-resistance determinants were maintained in the genome. It is proposed that they might contribute some accessory genes to related dominant resistance determinants, for instance periplasmic metal-binding proteins or two-component regulatory systems. Alternatively, they may only remain in the genome because their DNA serves as a scaffold for the nucleoid. Using C. metallidurans as an example, this study sheds light on the fate and function of horizontally acquired genes in bacteria.

Ralstonia solanacearum Effectors Localize to Diverse Organelles in Solanum Hosts

Phytopathogenic bacteria secrete type III effector (T3E) proteins directly into host plant cells. T3Es can interact with plant proteins and frequently manipulate plant host physiological or developmental processes. The proper subcellular localization of T3Es is critical for their ability to interact with plant targets, and knowledge of T3E localization can be informative for studies of effector function. Here we investigated the subcellular localization of 19 T3Es from the phytopathogenic bacteria Ralstonia pseudosolanacearum and Ralstonia solanacearum. Approximately 45% of effectors in our library localize to both the plant cell periphery and the nucleus, 15% exclusively to the cell periphery, 15% exclusively to the nucleus, and 25% to other organelles, including tonoplasts and peroxisomes. Using tomato hairy roots, we show that T3E localization is similar in both leaves and roots and is not impacted by Solanum species. We find that in silico prediction programs are frequently inaccurate, highlighting the value of in planta localization experiments. Our data suggest that Ralstonia targets a wide diversity of cellular organelles and provides a foundation for developing testable hypotheses about Ralstonia effector function.

Molecular Diversity and Pathogenicity of Ralstonia solanacearum Species Complex Associated With Bacterial Wilt of Potato in Rwanda

Bacterial wilt (BW), caused by Ralstonia solanacearum species complex (RSSC), leads to substantial potato yield losses in Rwanda. Studies were conducted to (i) determine the molecular diversity of RSSC strains associated with BW of potato, (ii) generate an RSSC distribution map for epidemiological inferences, and (iii) test the pathogenicity of predominant RSSC phylotypes on six commercial potato cultivars. In surveys conducted in 2018 and 2019, tubers from wilting potato plants were collected for pathogen isolation. DNA was extracted from 95 presumptive RSSC strain colonies. The pathogen was phylotyped by multiplex PCR and typed at sequevar level. Phylotype II sequevar 1 strains were then haplotyped using multilocus tandem repeat sequence typing (TRST) schemes. Pathogenicity of one phylotype II strain and two phylotype III strains were tested on cultivars Kinigi, Kirundo, Victoria, Kazeneza, Twihaze, and Cruza. Two RSSC phylotypes were identified, phylotype II (95.79%, n = 91) and phylotype III (4.21%, n = 4). This is the first report of phylotype III strains from Rwanda. Phylotype II strains were identified as sequevar 1 and distributed across potato growing regions in the country. The TRST scheme identified 14 TRST haplotypes within the phylotype II sequevar 1 strains with moderate diversity index (HGDI = 0.55). Mapping of TRST haplotypes revealed that a single TRST '8-5-12-7-5' haplotype plays an important epidemiological role in BW of potato in Rwanda. None of the cultivars had complete resistance to the tested phylotypes; the level of susceptibility varied among cultivars. Cultivar Cruza, which is less susceptible to phylotype II and III strains, is recommended when planting potatoes in the fields with history of BW.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Leaf-to-Whole Plant Spread Bioassay for Pepper and Ralstonia solanacearum Interaction Determines Inheritance of Resistance to Bacterial Wilt for Further Breeding

Bacterial wilt (BW) disease from Ralstonia solanacearum is a serious disease and causes severe yield losses in chili peppers worldwide. Resistant cultivar breeding is the most effective in controlling BW. Thus, a simple and reliable evaluation method is required to assess disease severity and to investigate the inheritance of resistance for further breeding programs. Here, we developed a reliable leaf-to-whole plant spread bioassay for evaluating BW disease and then, using this, determined the inheritance of resistance to R. solanacearum in peppers. Capsicum annuum ‘MC4′ displayed a completely resistant response with fewer disease symptoms, a low level of bacterial cell growth, and significant up-regulations of defense genes in infected leaves compared to those in susceptible ‘Subicho’. We also observed the spreading of wilt symptoms from the leaves to the whole susceptible plant, which denotes the normal BW wilt symptoms, similar to the drenching method. Through this, we optimized the evaluation method of the resistance to BW. Additionally, we performed genetic analysis for resistance inheritance. The parents, F₁ and 90 F₂ progenies, were evaluated, and the two major complementary genes involved in the BW resistance trait were confirmed. These could provide an accurate evaluation to improve resistant pepper breeding efficiency against BW.

Genomic sequencing of different sequevars of Ralstonia solanacearum belonging to the Moko ecotype

Banana vascular wilt or Moko is a disease caused by Ralstonia solanacearum. This study aimed to sequence, assemble, annotate, and compare the genomes of R. solanacearum Moko ecotypes of different sequevar strains from Brazil. Average nucleotide identity analyses demonstrated a high correlation (> 96%) between the genome sequences of strains CCRMRs277 (sequevar IIA-24), CCRMRs287 (IIB-4), CCRMRs304 (IIA-24), and CCRMRsB7 (IIB-25), which were grouped into phylotypes IIA and IIB. The number of coding sequences present in chromosomes and megaplasmids varied from 3,070 to 3,521 and 1,669 to 1,750, respectively. Pangenome analysis identified 3,378 clusters in the chromosomes, of which 2,604 were shared by all four analyzed genomes and 2,580 were single copies. In megaplasmids, 1,834 clusters were identified, of which 1,005 were shared by all four genomes and 992 were identified as single copies. Strains CCRMRsB7 and CCRMRs287 differed from the others by having unique clusters in both their chromosomes and megaplasmids, and CCRMRsB7 possessed the largest genome among all Moko ecotype strains sequenced to date. Therefore, the genomic information obtained in this study provides a theoretical basis for the identification, characterization, and phylogenetic analysis of R. solanacearum Moko ecotypes.

Taxonomy and Phylogenetic Research on Ralstonia solanacearum Species Complex: A Complex Pathogen with Extraordinary Economic Consequences

The bacterial wilt pathogen, first known as Bacillus solanacearum, has undergone numerous taxonomic changes since its first description in 1896. The history and significance of this pathogen is covered in this review with an emphasis on the advances in technology that were used to support each reclassification that finally led to the current separation of Ralstonia solanacearum into three genomic species. Frequent name changes occurred as methodology transitioned from phenotypic, biochemical, and molecular studies, to genomics and functional genomics. The diversity, wide host range, and geographical distribution of the bacterial wilt pathogen resulted in its division into three species as genomic analyses elucidated phylogenetic relationships among strains. Current advances in phylogenetics and functional genomics now open new avenues for research into epidemiology and control of the devastating bacterial wilt disease.

A quick and efficient hydroponic potato infection method for evaluating potato resistance and Ralstonia solanacearum virulence

BACKGROUND

Potato, the third most important crop worldwide, plays a critical role in human food security. Brown rot, one of the most destructive potato diseases caused by Ralstonia solanacearum, results in huge economic losses every year. A quick, stable, low cost and high throughout method is required to meet the demands of identification of germplasm resistance to bacterial wilt in potato breeding programs.

RESULTS

Here we present a novel R. solanacearum hydroponic infection assay on potato plants grown in vitro. Through testing wilt symptom appearance and bacterial colonization in aerial part of plants, we found that the optimum conditions for in vitro potato infection were using an OD₆₀₀ 0.01 bacterial solution suspended with tap water for infection, broken potato roots and an open container. Infection using R. solanacearum strains with differential degree of aggressivity demonstrated that this infection system is equally efficient as soil-drench inoculation for assessment of R. solanacearum virulence on potato. A small-scale assessment of 32 potato germplasms identified three varieties highly resistant to the pathogen, which indicates this infection system is a useful method for high-throughout screening of potato germplasm for resistance. Furthermore, we demonstrate the utility of a strain carrying luminescence to easily quantify bacterial colonization and the detection of latent infections in hydroponic conditions, which can be efficiently used in potato breeding programs.

CONCLUSIONS

We have established a quick and efficient in vitro potato infection system, which may facilitate breeding for new potato cultivars with high resistance to R. solanacearum.

Complete Genome Sequence of Sequevar 14M Ralstonia solanacearum Strain HA4-1 Reveals Novel Type III Effectors Acquired Through Horizontal Gene Transfer

Ralstonia solanacearum, which causes bacterial wilt in a broad range of plants, is considered a "species complex" due to its significant genetic diversity. Recently, we have isolated a new R. solanacearum strain HA4-1 from Hong'an county in Hubei province of China and identified it being phylotype I, sequevar 14M (phylotype I-14M). Interestingly, we found that it can cause various disease symptoms among different potato genotypes and display different pathogenic behavior compared to a phylogenetically related strain, GMI1000. To dissect the pathogenic mechanisms of HA4-1, we sequenced its whole genome by combined sequencing technologies including Illumina HiSeq2000, PacBio RS II, and BAC-end sequencing. Genome assembly results revealed the presence of a conventional chromosome, a megaplasmid as well as a 143 kb plasmid in HA4-1. Comparative genome analysis between HA4-1 and GMI1000 shows high conservation of the general virulence factors such as secretion systems, motility, exopolysaccharides (EPS), and key regulatory factors, but significant variation in the repertoire and structure of type III effectors, which could be the determinants of their differential pathogenesis in certain potato species or genotypes. We have identified two novel type III effectors that were probably acquired through horizontal gene transfer (HGT). These novel R. solanacearum effectors display homology to several YopJ and XopAC family members. We named them as RipBR and RipBS. Notably, the copy of RipBR on the plasmid is a pseudogene, while the other on the megaplasmid is normal. For RipBS, there are three copies located in the megaplasmid and plasmid, respectively. Our results have not only enriched the genome information on R. solanacearum species complex by sequencing the first sequevar 14M strain and the largest plasmid reported in R. solanacearum to date but also revealed the variation in the repertoire of type III effectors. This will greatly contribute to the future studies on the pathogenic evolution, host adaptation, and interaction between R. solanacearum and potato.

Prediction of Host-Specific Genes by Pan-Genome Analyses of the Korean Ralstonia solanacearum Species Complex

The soil-borne pathogenic Ralstonia solanacearum species complex (RSSC) is a group of plant pathogens that is economically destructive worldwide and has a broad host range, including various solanaceae plants, banana, ginger, sesame, and clove. Previously, Korean RSSC strains isolated from samples of potato bacterial wilt were grouped into four pathotypes based on virulence tests against potato, tomato, eggplant, and pepper. In this study, we sequenced the genomes of 25 Korean RSSC strains selected based on these pathotypes. The newly sequenced genomes were analyzed to determine the phylogenetic relationships between the strains with average nucleotide identity values, and structurally compared via multiple genome alignment using Mauve software. To identify candidate genes responsible for the host specificity of the pathotypes, functional genome comparisons were conducted by analyzing pan-genome orthologous group (POG) and type III secretion system effectors (T3es). POG analyses revealed that a total of 128 genes were shared only in tomato-non-pathogenic strains, 8 genes in tomato-pathogenic strains, 5 genes in eggplant-non-pathogenic strains, 7 genes in eggplant-pathogenic strains, 1 gene in pepper-non-pathogenic strains, and 34 genes in pepper-pathogenic strains. When we analyzed T3es, three host-specific effectors were predicted: RipS3 (SKWP3) and RipH3 (HLK3) were found only in tomato-pathogenic strains, and RipAC (PopC) were found only in eggplant-pathogenic strains. Overall, we identified host-specific genes and effectors that may be responsible for virulence functions in RSSC in silico. The expected characters of those genes suggest that the host range of RSSC is determined by the comprehensive actions of various virulence factors, including effectors, secretion systems, and metabolic enzymes.

Transcriptome and Comparative Genomics Analyses Reveal New Functional Insights on Key Determinants of Pathogenesis and Interbacterial Competition in Pectobacterium and Dickeya spp

Soft-rot Enterobacteriaceae (SRE), typified by Pectobacterium and Dickeya genera, are phytopathogenic bacteria inflicting soft-rot disease in crops worldwide. By combining genomic information from 100 SRE with whole-transcriptome data sets, we identified novel genomic and transcriptional associations among key pathogenicity themes in this group. Comparative genomics revealed solid linkage between the type I secretion system (T1SS) and the carotovoricin bacteriophage (Ctv) conserved in 96.7% of Pectobacterium genomes. Moreover, their coactivation during infection indicates a novel functional association involving T1SS and Ctv. Another bacteriophage-borne genomic region, mostly confined to less than 10% of Pectobacterium strains, was found, presumably comprising a novel lineage-specific prophage in the genus. We also detected the transcriptional coregulation of a previously predicted toxin/immunity pair (WHH and SMI1_KNR4 families), along with the type VI secretion system (T6SS), which includes hcp and/or vgrG genes, suggesting a role in disease development as T6SS-dependent effectors. Further, we showed that another predicted T6SS-dependent endonuclease (AHH family) exhibited toxicity in ectopic expression assays, indicating antibacterial activity. Additionally, we report the striking conservation of the group 4 capsule (GFC) cluster in 100 SRE strains which consistently features adjacently conserved serotype-specific gene arrays comprising a previously unknown organization in GFC clusters. Also, extensive sequence variations found in gfcA orthologs suggest a serotype-specific role in the GfcABCD machinery.IMPORTANCE Despite the considerable loss inflicted on important crops yearly by Pectobacterium and Dickeya diseases, investigations on key virulence and interbacterial competition assets relying on extensive comparative genomics are still surprisingly lacking for these genera. Such approaches become more powerful over time, underpinned by the growing amount of genomic information in public databases. In particular, our findings point to new functional associations among well-known genomic themes enabling alternative means of neutralizing SRE diseases through disruption of pivotal virulence programs. By elucidating novel transcriptional and genomic associations, this study adds valuable information on virulence candidates that could be decisive in molecular applications in the near future. The utilization of 100 genomes of Pectobacterium and Dickeya strains in this study is unprecedented for comparative analyses in these taxa, and it provides novel insights on the biology of economically important plant pathogens.

Using single nucleotide variations in single-cell RNA-seq to identify subpopulations and genotype-phenotype linkage

Despite its popularity, characterization of subpopulations with transcript abundance is subject to a significant amount of noise. We propose to use effective and expressed nucleotide variations (eeSNVs) from scRNA-seq as alternative features for tumor subpopulation identification. We develop a linear modeling framework, SSrGE, to link eeSNVs associated with gene expression. In all the datasets tested, eeSNVs achieve better accuracies than gene expression for identifying subpopulations. Previously validated cancer-relevant genes are also highly ranked, confirming the significance of the method. Moreover, SSrGE is capable of analyzing coupled DNA-seq and RNA-seq data from the same single cells, demonstrating its value in integrating multi-omics single cell techniques. In summary, SNV features from scRNA-seq data have merits for both subpopulation identification and linkage of genotype-phenotype relationship.

Genetic Diversity of the Ralstonia solanacearum Species Complex in the Southwest Indian Ocean Islands

Epidemiological surveillance of plant pathogens based on genotyping methods is mandatory to improve disease management strategies. In the Southwest Indian Ocean (SWIO) islands, bacterial wilt (BW) caused by the Ralstonia solanacearum species complex (RSSC) is hampering the production of many sustainable and cash crops. To thoroughly analyze the genetic diversity of the RSSC in the SWIO, we performed a wide sampling survey (in Comoros, Mauritius, Reunion, Rodrigues, and Seychelles) that yielded 1,704 isolates from 129 plots, mainly from solanaceous crops. Classification of the isolates to the four major RSSC phylogenetic groups, named phylotypes, showed that 87% were phylotype I, representing the most prevalent strain in each of the SWIO islands. Additionally, 9.7% were phylotype II, and 3.3% were phylotype III; however, these isolates were found only in Reunion. Phylotype IV (2 isolates), known to be restricted to Indonesia-Australia-Japan, was reported in Mauritius, representing the first report of this group in the SWIO. Partial endoglucanase (egl) sequencing, based on the selection of 145 isolates covering the geographic and host diversity in the SWIO (also including strains from Mayotte and Madagascar), revealed 14 sequevars with Reunion and Mauritius displaying the highest sequevar diversity. Through a multilocus sequence analysis (MLSA) scheme based on the partial sequencing of 6 housekeeping genes (gdhA, gyrB, rplB, leuS, adk, and mutS) and 1 virulence-associated gene (egl), we inferred the phylogenetic relationships between these 145 SWIO isolates and 90 worldwide RSSC reference strains. Phylotype I was the most recombinogenic, although recombination events were detected among all phylotypes. A multilocus sequence typing (MLST) scheme identified 29 sequence types (STs) with variable geographic distributions in the SWIO. The outstanding epidemiologic feature was STI-13 (sequevar I-31), which was overrepresented in the SWIO and obviously reflected a lineage strongly adapted to the SWIO environment. A goeBURST analysis identified eight clonal complexes (CCs) including SWIO isolates, four CCs being geographically restricted to the SWIO, and four CCs being widespread beyond the SWIO. This work, which highlights notable genetic links between African and SWIO strains, provides a basis for the epidemiological surveillance of RSSC and will contribute to BW management in the SWIO.

Tube-Wise Diagnostic Microarray for the Multiplex Characterization of the Complex Plant Pathogen Ralstonia solanacearum

Ralstonia solanacearum is a well-known agricultural and ecological threat worldwide. The complexity of the R. solanacearum species complex (Rssc) represents a challenge for the accurate characterization of epidemiological strains by official services and research laboratories. The majority of protocols only focus on a narrow range of strains; however, this species complex includes strains that represent major constraints and are under strict regulation. The main drawback associated with the current methods of detecting and characterizing Rssc strains is their reliance on combining different protocols to properly characterize the strains at the ecotype level, which require time and money. Therefore, we used microarray technology (ArrayTube) to develop a standard protocol, which characterizes 17 major groups of interest in the Rssc, in a single multiplex reaction. These 17 majors groups are linked with a phylogenetic assignation (phylotypes, sequevars), but also with an ecotype assignation associated with a range of hosts (e.g., brown rot, Moko). Probes were designed with a 50-mer length constraint and thoroughly evaluated for any flaws or secondary structures. The strains are characterized based on a DNA extraction from pure culture. Validation data showed strong intra-repeatability, inter-repeatability, and reproducibility as well as good specificity. A hierarchical analysis of the probe groups is suitable for an accurate characterization. Compared with single marker detection tests, the method described in this paper addresses efficiently the issue of combining several tests by testing a large number of phylogenetic markers in a single reaction assay. This custom microarray (RsscAT) represents a significant improvement in the epidemiological monitoring of Rssc strains worldwide, and it has the potential to provide insights for phylogenetic incongruence of Rssc strains based on the host of isolation and may be used to indicate potentially emergent strains.

Comparative Analysis of Ralstonia solanacearum Methylomes

Ralstonia solanacearum is an important soil-borne plant pathogen with broad geographical distribution and the ability to cause wilt disease in many agriculturally important crops. Genome sequencing of multiple R. solanacearum strains has identified both unique and shared genetic traits influencing their evolution and ability to colonize plant hosts. Previous research has shown that DNA methylation can drive speciation and modulate virulence in bacteria, but the impact of epigenetic modifications on the diversification and pathogenesis of R. solanacearum is unknown. Sequencing of R. solanacearum strains GMI1000 and UY031 using Single Molecule Real-Time technology allowed us to perform a comparative analysis of R. solanacearum methylomes. Our analysis identified a novel methylation motif associated with a DNA methylase that is conserved in all complete Ralstonia spp. genomes and across the Burkholderiaceae, as well as a methylation motif associated to a phage-borne methylase unique to R. solanacearum UY031. Comparative analysis of the conserved methylation motif revealed that it is most prevalent in gene promoter regions, where it displays a high degree of conservation detectable through phylogenetic footprinting. Analysis of hyper- and hypo-methylated loci identified several genes involved in global and virulence regulatory functions whose expression may be modulated by DNA methylation. Analysis of genome-wide modification patterns identified a significant correlation between DNA modification and transposase genes in R. solanacearum UY031, driven by the presence of a high copy number of ISrso3 insertion sequences in this genome and pointing to a novel mechanism for regulation of transposition. These results set a firm foundation for experimental investigations into the role of DNA methylation in R. solanacearum evolution and its adaptation to different plants.

Identification and characterization of chemosensors for d-malate, unnatural enantiomer of malate, in Ralstonia pseudosolanacearum

Ralstonia pseudosolanacearum Ps29 is attracted by nonmetabolizable d-malate, an unnatural enantiomer. Screening of a complete collection of single-mcp-gene deletion mutants of Ps29 revealed that the RSc1156 homologue is a chemosensor for d-malate. An RSc1156 homologue deletion mutant of Ps29 showed decreased but significant responses to d-malate, suggesting the existence of another d-malate chemosensor. McpM previously had been identified as a chemosensor for l-malate. We constructed an RSc1156 homologue mcpM double deletion mutant and noted that this mutant failed to respond to d-malate; thus, the RSc1156 homologue and McpM are the major chemosensors for d-malate in this organism. To further characterize the ligand specificities of the RSc1156 homologue and McpM, we constructed a Ps29 derivative (designated K18) harbouring deletions in 18 individual mcp genes, including mcpM and RSc1156. K18 harbouring the RSc1156 homologue responded strongly to l-tartrate and d-malate and moderately to d-tartrate, but not to l-malate or succinate. K18 harbouring mcpM responded strongly to l-malate and d-tartrate and moderately to succinate, fumarate and d-malate. Ps29 utilizes l-malate and l-tartrate, but not d-malate. We therefore concluded that l-tartrate and l-malate are natural ligands of the RSc1156 homologue and McpM, respectively, and that chemotaxis toward d-malate is a fortuitous response by the RSc1156 homologue and McpM in Ps29. We propose re-designation of the RSc1156 homologue as McpT. In tomato plant infection assays, the mcpT deletion mutant of highly virulent R. pseudosolanacearum MAFF106611 was as infectious as wild-type MAFF106611, suggesting that McpT-mediated chemotaxis does not play an important role in tomato plant infection.

A Resource Allocation Trade-Off between Virulence and Proliferation Drives Metabolic Versatility in the Plant Pathogen Ralstonia solanacearum

Bacterial pathogenicity relies on a proficient metabolism and there is increasing evidence that metabolic adaptation to exploit host resources is a key property of infectious organisms. In many cases, colonization by the pathogen also implies an intensive multiplication and the necessity to produce a large array of virulence factors, which may represent a significant cost for the pathogen. We describe here the existence of a resource allocation trade-off mechanism in the plant pathogen R. solanacearum. We generated a genome-scale reconstruction of the metabolic network of R. solanacearum, together with a macromolecule network module accounting for the production and secretion of hundreds of virulence determinants. By using a combination of constraint-based modeling and metabolic flux analyses, we quantified the metabolic cost for production of exopolysaccharides, which are critical for disease symptom production, and other virulence factors. We demonstrated that this trade-off between virulence factor production and bacterial proliferation is controlled by the quorum-sensing-dependent regulatory protein PhcA. A phcA mutant is avirulent but has a better growth rate than the wild-type strain. Moreover, a phcA mutant has an expanded metabolic versatility, being able to metabolize 17 substrates more than the wild-type. Model predictions indicate that metabolic pathways are optimally oriented towards proliferation in a phcA mutant and we show that this enhanced metabolic versatility in phcA mutants is to a large extent a consequence of not paying the cost for virulence. This analysis allowed identifying candidate metabolic substrates having a substantial impact on bacterial growth during infection. Interestingly, the substrates supporting well both production of virulence factors and growth are those found in higher amount within the plant host. These findings also provide an explanatory basis to the well-known emergence of avirulent variants in R. solanacearum populations in planta or in stressful environments.

Metabolic versatility is a critical element for pathogen’s virulence and their ability to survive in the host. Beyond the necessity to collect resources during infection, pathogens face a resource allocation dilemma: they have to use nutritional resources to proliferate inside the host, and in the other hand they need to mobilize matter and energy for the production of essential virulence factors. In this study, we provide evidence of that such a trade-off constrains antagonistically bacterial proliferation and virulence in the bacterial plant pathogen Ralstonia solanacearum. We determined the energetic cost required by R. solanacearum to produce and secrete exopolysaccharide, which is a major virulence factor required for wilting symptom appearance. We validated this result by showing that bacterial mutants defective for exopolysaccharide production or other virulence factor indeed have an increased growth rate compared to the wild-type strain. We provide evidence that this trade-off mechanism is orchestrated by the phcA master regulatory gene, which directly connects quorum-sensing regulation to metabolic versatility and virulence. Our results also support the view that R. solanacearum specializes towards a restricted number of substrates used during in planta growth.

The effects of host age and spatial location on bacterial community composition in the English Oak tree (Quercus robur)

Drivers of bacterial community assemblages associated with plants are diverse and include biotic factors, such as competitors and host traits, and abiotic factors, including environmental conditions and dispersal mechanisms. We examine the roles of spatial distribution and host size, as an approximation for age, in shaping the microbiome associated with Quercus robur woody tissue using culture-independent 16S rRNA gene amplicon sequencing. In addition to providing a baseline survey of the Q. robur microbiome, we screened for the pathogen of acute oak decline. Our results suggest that age is a predictor of bacterial community composition, demonstrating a surprising negative correlation between tree age and alpha diversity. We find no signature of dispersal limitation within the Wytham Woods plot sampled. Together, these results provide evidence for niche-based hypotheses of community assembly and the importance of tree age in bacterial community structure, as well as highlighting that caution must be applied when diagnosing dysbiosis in a long-lived plant host.

Escaping Underground Nets: Extracellular DNases Degrade Plant Extracellular Traps and Contribute to Virulence of the Plant Pathogenic Bacterium Ralstonia solanacearum

Plant root border cells have been recently recognized as an important physical defense against soil-borne pathogens. Root border cells produce an extracellular matrix of protein, polysaccharide and DNA that functions like animal neutrophil extracellular traps to immobilize pathogens. Exposing pea root border cells to the root-infecting bacterial wilt pathogen Ralstonia solanacearum triggered release of DNA-containing extracellular traps in a flagellin-dependent manner. These traps rapidly immobilized the pathogen and killed some cells, but most of the entangled bacteria eventually escaped. The R. solanacearum genome encodes two putative extracellular DNases (exDNases) that are expressed during pathogenesis, suggesting that these exDNases contribute to bacterial virulence by enabling the bacterium to degrade and escape root border cell traps. We tested this hypothesis with R. solanacearum deletion mutants lacking one or both of these nucleases, named NucA and NucB. Functional studies with purified proteins revealed that NucA and NucB are non-specific endonucleases and that NucA is membrane-associated and cation-dependent. Single ΔnucA and ΔnucB mutants and the ΔnucA/B double mutant all had reduced virulence on wilt-susceptible tomato plants in a naturalistic soil-soak inoculation assay. The ΔnucA/B mutant was out-competed by the wild-type strain in planta and was less able to stunt root growth or colonize plant stems. Further, the double nuclease mutant could not escape from root border cells in vitro and was defective in attachment to pea roots. Taken together, these results demonstrate that extracellular DNases are novel virulence factors that help R. solanacearum successfully overcome plant defenses to infect plant roots and cause bacterial wilt disease.

Escaping Underground Nets: Extracellular DNases Degrade Plant Extracellular Traps and Contribute to Virulence of the Plant Pathogenic Bacterium Ralstonia solanacearum

Plant root border cells have been recently recognized as an important physical defense against soil-borne pathogens. Root border cells produce an extracellular matrix of protein, polysaccharide and DNA that functions like animal neutrophil extracellular traps to immobilize pathogens. Exposing pea root border cells to the root-infecting bacterial wilt pathogen Ralstonia solanacearum triggered release of DNA-containing extracellular traps in a flagellin-dependent manner. These traps rapidly immobilized the pathogen and killed some cells, but most of the entangled bacteria eventually escaped. The R. solanacearum genome encodes two putative extracellular DNases (exDNases) that are expressed during pathogenesis, suggesting that these exDNases contribute to bacterial virulence by enabling the bacterium to degrade and escape root border cell traps. We tested this hypothesis with R. solanacearum deletion mutants lacking one or both of these nucleases, named NucA and NucB. Functional studies with purified proteins revealed that NucA and NucB are non-specific endonucleases and that NucA is membrane-associated and cation-dependent. Single ΔnucA and ΔnucB mutants and the ΔnucA/B double mutant all had reduced virulence on wilt-susceptible tomato plants in a naturalistic soil-soak inoculation assay. The ΔnucA/B mutant was out-competed by the wild-type strain in planta and was less able to stunt root growth or colonize plant stems. Further, the double nuclease mutant could not escape from root border cells in vitro and was defective in attachment to pea roots. Taken together, these results demonstrate that extracellular DNases are novel virulence factors that help R. solanacearum successfully overcome plant defenses to infect plant roots and cause bacterial wilt disease.

Genomic and proteomic evidence supporting the division of the plant pathogen Ralstonia solanacearum into three species

Background

The increased availability of genome sequences has advanced the development of genomic distance methods to describe bacterial diversity. Results of these fast-evolving methods are highly correlated with those of the historically standard DNA-DNA hybridization technique. However, these genomic-based methods can be done more rapidly and less expensively and are less prone to technical and human error. They are thus a technically accessible replacement for species delineation. Here, we use several genomic comparison methods, supported by our own proteomic analyses and metabolic characterization as well as previously published DNA-DNA hybridization analyses, to differentiate members of the Ralstonia solanacearum species complex into three species. This pathogen group consists of diverse and widespread strains that cause bacterial wilt disease on many different plants.

Results

We used three different methods to compare the complete genomes of 29 strains from the R. solanacearum species complex. In parallel we profiled the proteomes of 73 strains using Matrix-Assisted Laser Desorption/Ionization-Time of Flight Mass Spectrometry (MALDI-TOF-MS). Proteomic profiles together with genomic sequence comparisons consistently and comprehensively described the diversity of the R. solanacearum species complex. In addition, genome-driven functional phenotypic assays excitingly supported an old hypothesis (Hayward et al. (J Appl Bacteriol 69:269–80, 1990)), that closely related members of the R. solanacearum could be identified through a simple assay of anaerobic nitrate metabolism. This assay allowed us to clearly and easily differentiate phylotype II and IV strains from phylotype I and III strains. Further, genomic dissection of the pathway distinguished between proposed subspecies within the current phylotype IV. The assay revealed large scale differences in energy production within the R. solanacearum species complex, indicating coarse evolutionary distance and further supporting a repartitioning of this group into separate species.

Conclusions

Together, the results of these studies support the proposed division of the R. solanacearum species complex into three species, consistent with recent literature, and demonstrate the utility of proteomic and genomic approaches to delineate bacterial species.

Electronic supplementary material

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

Microsatellites in Pursuit of Microbial Genome Evolution

Microsatellites or short sequence repeats are widespread genetic markers which are hypermutable 1-6 bp long short nucleotide motifs. Significantly, their applications in genetics are extensive due to their ceaseless mutational degree, widespread length variations and hypermutability skills. These features make them useful in determining the driving forces of evolution by using powerful molecular techniques. Consequently, revealing important questions, for example, what is the significance of these abundant sequences in DNA, what are their roles in genomic evolution? The answers of these important questions are hidden in the ways these short motifs contributed in altering the microbial genomes since the origin of life. Even though their size ranges from 1 -to- 6 bases, these repeats are becoming one of the most popular genetic probes in determining their associations and phylogenetic relationships in closely related genomes. Currently, they have been widely used in molecular genetics, biotechnology and evolutionary biology. However, due to limited knowledge; there is a significant gap in research and lack of information concerning hypermutational mechanisms. These mechanisms play a key role in microsatellite loci point mutations and phase variations. This review will extend the understandings of impacts and contributions of microsatellite in genomic evolution and their universal applications in microbiology.

TALE-Like Effectors Are an Ancestral Feature of the Ralstonia solanacearum Species Complex and Converge in DNA Targeting Specificity

Ralstonia solanacearum, a species complex of bacterial plant pathogens divided into four monophyletic phylotypes, causes plant diseases in tropical climates around the world. Some strains exhibit a broad host range on solanaceous hosts, while others are highly host-specific as for example some banana-pathogenic strains. Previous studies showed that transcription activator-like (TAL) effectors from Ralstonia, termed RipTALs, are capable of activating reporter genes in planta, if these are preceded by a matching effector binding element (EBE). RipTALs target DNA via their central repeat domain (CRD), where one repeat pairs with one DNA-base of the given EBE. The repeat variable diresidue dictates base repeat specificity in a predictable fashion, known as the TALE code. In this work, we analyze RipTALs across all phylotypes of the Ralstonia solanacearum species complex. We find that RipTALs are prevalent in phylotypes I and IV but absent from most phylotype III and II strains (10/12, 8/14, 1/24, and 1/5 strains contained a RipTAL, respectively). RipTALs originating from strains of the same phylotype show high levels of sequence similarity (>98%) in the N-terminal and C-terminal regions, while RipTALs isolated from different phylotypes show 47-91% sequence similarity in those regions, giving rise to four RipTAL classes. We show that, despite sequence divergence, the base preference for guanine, mediated by the N-terminal region, is conserved across RipTALs of all classes. Using the number and order of repeats found in the CRD, we functionally sub-classify RipTALs, introduce a new simple nomenclature, and predict matching EBEs for all seven distinct RipTALs identified. We experimentally study RipTAL EBEs and uncover that some RipTALs are able to target the EBEs of other RipTALs, referred to as cross-reactivity. In particular, RipTALs from strains with a broad host range on solanaceous hosts cross-react on each other's EBEs. Investigation of sequence divergence between RipTAL repeats allows for a reconstruction of repeat array biogenesis, for example through slipped strand mispairing or gene conversion. Using these studies we show how RipTALs of broad host range strains evolved convergently toward a shared target sequence. Finally, we discuss the differences between TALE-likes of plant pathogens in the context of disease ecology.

Towards the Identification of Type III Effectors Associated with Ralstonia solanacearum Virulence on Tomato and Eggplant

For the development of pathogen-informed breeding strategies, identifying the microbial genes involved in interactions with the plant is a critical step. To identify type III effector (T3E) repertoires associated with virulence of the bacterial wilt pathogen Ralstonia solanacearum on Solanaceous crops, we used an original association genetics approach combining DNA microarray data and pathogenicity data on resistant eggplant, pepper, and tomato accessions. From this first screen, 25 T3Es were further full-length polymerase chain reaction-amplified within a 35-strain field collection, to assess their distribution and allelic diversity. Six T3E repertoire groups were identified, within which 11 representative strains were chosen to challenge the bacterial wilt-resistant egg plants 'Dingras multiple Purple' and 'AG91-25', and tomato Hawaii 7996. The virulence or avirulence phenotypes could not be explained by specific T3E repertoires, but rather by individual T3E genes. We identified seven highly avirulence-associated genes, among which ripP2, primarily referenced as conferring avirulence to Arabidopsis thaliana. Interestingly, no T3E was associated with avirulence to both egg-plants. Highly virulence-associated genes were also identified: ripA5_2, ripU, and ripV2. This study should be regarded as a first step toward investigating both avirulence and virulence function of the highlighted genes, but also their evolutionary dynamics in natural R. solanacearum populations.

Identification of the mcpA and mcpM genes, encoding methyl-accepting proteins involved in amino acid and l-malate chemotaxis, and involvement of McpM-mediated chemotaxis in plant infection by Ralstonia pseudosolanacearum (formerly Ralstonia solanacearum phylotypes I and III)

Sequence analysis has revealed the presence of 22 putative methyl-accepting chemotaxis protein (mcp) genes in the Ralstonia pseudosolanacearum GMI1000 genome. PCR analysis and DNA sequencing showed that the highly motile R. pseudosolanacearum strain Ps29 possesses homologs of all 22 R. pseudosolanacearum GMI1000 mcp genes. We constructed a complete collection of single mcp gene deletion mutants of R. pseudosolanacearum Ps29 by unmarked gene deletion. Screening of the mutant collection revealed that R. pseudosolanacearum Ps29 mutants of RSp0507 and RSc0606 homologs were defective in chemotaxis to l-malate and amino acids, respectively. RSp0507 and RSc0606 homologs were designated mcpM and mcpA. While wild-type R. pseudosolanacearum strain Ps29 displayed attraction to 16 amino acids, the mcpA mutant showed no response to 12 of these amino acids and decreased responses to 4 amino acids. We constructed mcpA and mcpM deletion mutants of highly virulent R. pseudosolanacearum strain MAFF106611 to investigate the contribution of chemotaxis to l-malate and amino acids to tomato plant infection. Neither single mutant exhibited altered virulence for tomato plants when tested by root dip inoculation assays. In contrast, the mcpM mutant (but not the mcpA mutant) was significantly less infectious than the wild type when tested by a sand soak inoculation assay, which requires bacteria to locate and invade host roots from sand. Thus, McpM-mediated chemotaxis, possibly reflecting chemotaxis to l-malate, facilitates R. pseudosolanacearum motility to tomato roots in sand.

Comparative Transcriptome Analysis Reveals Cool Virulence Factors of Ralstonia solanacearum Race 3 Biovar 2

While most strains of the plant pathogenic bacterium Ralstonia solanacearum are tropical, the race 3 biovar 2 (R3bv2) subgroup attacks plants in cooler climates. To identify mechanisms underlying this trait, we compared the transcriptional profiles of R. solanacearum R3bv2 strain UW551 and tropical strain GMI1000 at 20°C and 28°C, both in culture and during tomato pathogenesis. 4.2% of the ORFs in the UW551 genome and 7.9% of the GMI1000 ORFs were differentially expressed by temperature in planta. The two strains had distinct transcriptional responses to temperature change. GMI1000 up-regulated several stress response genes at 20°C, apparently struggling to cope with plant defenses. At the cooler temperature, R3bv2 strain UW551 up-regulated a cluster encoding a mannose-fucose binding lectin, LecM; a quorum sensing-dependent protein, AidA; and a related hypothetical protein, AidC. The last two genes are absent from the GMI1000 genome. In UW551, all three genes were positively regulated by the adjacent SolI/R quorum sensing system. These temperature-responsive genes were required for full virulence in R3bv2. Mutants lacking lecM, aidA, or aidC were each significantly more reduced in virulence on tomato at 20°C than at 28°C in both a naturalistic soil soak inoculation assay and when they were inoculated directly into tomato stems. The lecM and aidC mutants also survived poorly in potato tubers at the seed tuber storage temperature of 4°C, and the lecM mutant was defective in biofilm formation in vitro. Together, these results suggest novel mechanisms, including a lectin, are involved in the unique temperate epidemiology of R3bv2.

Genome-Enabled Phylogeographic Investigation of the Quarantine Pathogen Ralstonia solanacearum Race 3 Biovar 2 and Screening for Sources of Resistance Against Its Core Effectors

Phylogeographic studies inform about routes of pathogen dissemination and are instrumental for improving import/export controls. Genomes of 17 isolates of the bacterial wilt and potato brown rot pathogen Ralstonia solanacearum race 3 biovar 2 (R3bv2), a Select Agent in the United States, were thus analyzed to get insight into the phylogeography of this pathogen. Thirteen of fourteen isolates from Europe, Africa, and Asia were found to belong to a single clonal lineage while isolates from South America were genetically diverse and tended to carry ancestral alleles at the analyzed genomic loci consistent with a South American origin of R3bv2. The R3bv2 isolates share a core repertoire of 31 type III-secreted effector genes representing excellent candidates to be targeted with resistance genes in breeding programs to develop durable disease resistance. Toward this goal, 27 R3bv2 effectors were tested in eggplant, tomato, pepper, tobacco, and lettuce for induction of a hypersensitive-like response indicative of recognition by cognate resistance receptors. Fifteen effectors, eight of them core effectors, triggered a response in one or more plant species. These genotypes may harbor resistance genes that could be identified and mapped, cloned, and expressed in tomato or potato, for which sources of genetic resistance to R3bv2 are extremely limited.

A Multiplex PCR Assay to Detect and Differentiate Select Agent Strains of Ralstonia solanacearum

Ralstonia solanacearum race 3 biovar 2 strains are considered select agents by the U.S. government because they are not endemic to the United States and have the potential to cause brown rot in our potato production fields. Simple and accurate methods are needed for quick identification prior to more discriminating but time-consuming verification methods. We developed a multiplex PCR assay that identifies R. solanacearum species complex strains, signals whether the strain detected is a select agent, and controls for false negatives associated with PCR inhibition or unsuccessful DNA extractions in one reaction. We identified unique sequences of non-phage-related DNA for the R. solanacearum species complex strains, and for select agent strains, using in silico genome subtraction. We also designed and included an internal plant DNA control assay. Our multiplex PCR assay correctly identified 90 R. solanacearum species complex strains and 34 select agent strains, while not recognizing five out-group bacterial species. Additionally, the multiplex PCR assay facilitated the detection of plant DNA and R. solanacearum from infected tomato, potato, geranium, and tobacco plants. Our rapid, accurate, and reliable detection assay can help government officials make timely and appropriate recommendations to exclude this bacterium from the United States.

The C-terminal extension of PrhG impairs its activation of hrp expression and virulence in Ralstonia solanacearum

Ralstonia solanacearum is the second most destructive bacterial plant pathogens worldwide and HrpG is the master regulator of its pathogenicity. PrhG is a close paralogue of HrpG and both belong to OmpR/PhoB family of two-component response regulators. Despite a high similarity (72% global identity and 96% similarity in helix-loop-helix domain), they display distinct roles in pathogenicity. HrpG is necessary for the bacterial growth in planta and pathogenicity, while PrhG is dispensable for bacterial growth in planta and contributes little to pathogenicity. The main difference between HrpG and PrhG is the 50-amino-acid-long C-terminal extension in PrhG (amino-acid residues 230-283), which is absent in HrpG. When this extension is deleted, truncated PrhGs (under the control of its native promoter) allowed complete recovery of bacterial growth in planta and wild-type virulence of hrpG mutant. This novel finding demonstrates that the extension region in PrhG is responsible for the functional difference between HrpG and PrhG, which may block the binding of PrhG to target promoters and result in impaired activation of hrp expression by PrhG and reduced virulence of R. solanacearum.

Moko Disease-Causing Strains of Ralstonia solanacearum from Brazil Extend Known Diversity in Paraphyletic Phylotype II

The epidemic situation of Moko disease-causing strains in Latin America and Brazil is unclear. Thirty-seven Ralstonia solanacearum strains from Brazil that cause the Moko disease on banana and heliconia plants were sampled and phylogenetically typed using the endoglucanase (egl) and DNA repair (mutS) genes according to the phylotype and sequevar classification. All of the strains belonged to phylotype II and a portion of the strains was typed as the Moko disease-related sequevars IIA-6 and IIA-24. Nevertheless, two unsuspected sequevars also harbored the Moko disease-causing strains IIA-41 and IIB-25, and a new sequevar was described and named IIA-53. All of the strains were pathogenic to banana and some of the strains of sequevars IIA-6, IIA-24, and IIA-41 were also pathogenic to tomato. The Moko disease-causing strains from sequevar IIB-25 were pathogenic to potato but not to tomato. These results highlight the high diversity of strains of Moko in Brazil, reinforce the efficiency of the egl gene to reveal relationships among these strains, and contribute to a better understanding of the diversity of paraphyletic Moko disease-causing strains of the R. solanacearum species complex, where the following seven distinct genetic clusters have been described: IIA-6, IIA-24, IIA-41, IIA-53, IIB-3, IIB-4, and IIB-25.

Multihost experimental evolution of the pathogen Ralstonia solanacearum unveils genes involved in adaptation to plants

Ralstonia solanacearum, the causal agent of a lethal bacterial wilt plant disease, infects an unusually wide range of hosts. These hosts can further be split into plants where R. solanacearum is known to cause disease (original hosts) and those where this bacterium can grow asymptomatically (distant hosts). Moreover, this pathogen is able to adapt to many plants as supported by field observations reporting emergence of strains with enlarged pathogenic properties. To investigate the genetic bases of host adaptation, we conducted evolution experiments by serial passages of a single clone of the pathogen on three original and two distant hosts over 300 bacterial generations and then analyzed the whole-genome of nine evolved clones. Phenotypic analysis of the evolved clones showed that the pathogen can increase its fitness on both original and distant hosts although the magnitude of fitness increase was greater on distant hosts. Only few genomic modifications were detected in evolved clones compared with the ancestor but parallel evolutionary changes in two genes were observed in independent evolved populations. Independent mutations in the regulatory gene efpR were selected for in three populations evolved on beans, a distant host. Reverse genetic approaches confirmed that these mutations were associated with fitness gain on bean plants. This work provides a first step toward understanding the within-host evolutionary dynamics of R. solanacearum during infection and identifying bacterial genes subjected to in planta selection. The discovery of EfpR as a determinant conditioning host adaptation of the pathogen illustrates how experimental evolution coupled with whole-genome sequencing is a potent tool to identify novel molecular players involved in central life-history traits.

New type IV pili-related genes involved in early stages of Ralstonia solanacearum potato infection

This study provides insights into the pathogenesis of Ralstonia solanacearum, in particular with regards to strains belonging to phylotype IIB, sequevar 1 (IIB-1) and their interaction with potato, its natural host. We performed a comparative genomic analysis among IIB-1 R. solanacearum strains with different levels of virulence in order to identify candidate virulence genes. With this approach, we identified a 33.7-kb deletion in a strain showing reduced virulence on potato. This region contains a cluster of six genes putatively involved in type IV pili (Tfp) biogenesis. Functional analysis suggests that these proteins contribute to several Tfp-related functions such as twitching motility and biofilm formation. In addition, this genetic cluster was found to contribute to early bacterial wilt pathogenesis and colonization fitness of potato roots.

Polyphasic taxonomic revision of the Ralstonia solanacearum species complex: proposal to emend the descriptions of Ralstonia solanacearum and Ralstonia syzygii and reclassify current R. syzygii strains as Ralstonia syzygii subsp. syzygii subsp. nov., R. solanacearum phylotype IV strains as Ralstonia syzygii subsp. indonesiensis subsp. nov., banana blood disease bacterium strains as Ralstonia syzygii subsp. celebesensis subsp. nov. and R. solanacearum phylotype I and III strains as Ralstonia pseudosolanacearum sp. nov

The Ralstonia solanacearum species complex has long been recognized as a group of phenotypically diverse strains that can be subdivided into four phylotypes. Using a polyphasic taxonomic approach on an extensive set of strains, this study provides evidence for a taxonomic and nomenclatural revision of members of this complex. Data obtained from phylogenetic analysis of 16S-23S rRNA ITS gene sequences, 16S-23S rRNA intergenic spacer (ITS) region sequences and partial endoglucanase (egl) gene sequences and DNA-DNA hybridizations demonstrate that the R. solanacearum species complex comprises three genospecies. One of these includes the type strain of Ralstonia solanacearum and consists of strains of R. solanacearum phylotype II only. The second genospecies includes the type strain of Ralstonia syzygii and contains only phylotype IV strains. This genospecies is subdivided into three distinct groups, namely R. syzygii, the causal agent of Sumatra disease on clove trees in Indonesia, R. solanacearum phylotype IV strains isolated from different host plants mostly from Indonesia, and strains of the blood disease bacterium (BDB), the causal agent of the banana blood disease, a bacterial wilt disease in Indonesia that affects bananas and plantains. The last genospecies is composed of R. solanacearum strains that belong to phylotypes I and III. As these genospecies are also supported by phenotypic data that allow the differentiation of the three genospecies, the following taxonomic proposals are made: emendation of the descriptions of Ralstonia solanacearum and Ralstonia syzygii and descriptions of Ralstonia syzygii subsp. nov. (type strain R 001(T) = LMG 10661(T) = DSM 7385(T)) for the current R. syzygii strains, Ralstonia syzygii subsp. indonesiensis subsp. nov. (type strain UQRS 464(T) = LMG 27703(T) = DSM 27478(T)) for the current R. solanacearum phylotype IV strains, Ralstonia syzygii subsp. celebesensis subsp. nov. (type strain UQRS 627(T) = LMG 27706(T) = DSM 27477(T)) for the BDB strains and Ralstonia pseudosolanacearum sp. nov. (type strain UQRS 461(T) = LMG 9673(T) = NCPPB 1029(T)) for the strains of R. solanacearum phylotypes I and III.

Repertoire, unified nomenclature and evolution of the Type III effector gene set in the Ralstonia solanacearum species complex

Background

Ralstonia solanacearum is a soil-borne beta-proteobacterium that causes bacterial wilt disease in many food crops and is a major problem for agriculture in intertropical regions. R. solanacearum is a heterogeneous species, both phenotypically and genetically, and is considered as a species complex. Pathogenicity of R. solanacearum relies on the Type III secretion system that injects Type III effector (T3E) proteins into plant cells. T3E collectively perturb host cell processes and modulate plant immunity to enable bacterial infection.

Results

We provide the catalogue of T3E in the R. solanacearum species complex, as well as candidates in newly sequenced strains. 94 T3E orthologous groups were defined on phylogenetic bases and ordered using a uniform nomenclature. This curated T3E catalog is available on a public website and a bioinformatic pipeline has been designed to rapidly predict T3E genes in newly sequenced strains. Systematical analyses were performed to detect lateral T3E gene transfer events and identify T3E genes under positive selection. Our analyses also pinpoint the RipF translocon proteins as major discriminating determinants among the phylogenetic lineages.

Conclusions

Establishment of T3E repertoires in strains representatives of the R. solanacearum biodiversity allowed determining a set of 22 T3E present in all the strains but provided no clues on host specificity determinants. The definition of a standardized nomenclature and the optimization of predictive tools will pave the way to understanding how variation of these repertoires is correlated to the diversification of this species complex and how they contribute to the different strain pathotypes.

Ralstonia solanacearum requires PopS, an ancient AvrE-family effector, for virulence and To overcome salicylic acid-mediated defenses during tomato pathogenesis

During bacterial wilt of tomato, the plant pathogen Ralstonia solanacearum upregulates expression of popS, which encodes a type III-secreted effector in the AvrE family. PopS is a core effector present in all sequenced strains in the R. solanacearum species complex. The phylogeny of popS mirrors that of the species complex as a whole, suggesting that this is an ancient, vertically inherited effector needed for association with plants. A popS mutant of R. solanacearum UW551 had reduced virulence on agriculturally important Solanum spp., including potato and tomato plants. However, the popS mutant had wild-type virulence on a weed host, Solanum dulcamara, suggesting that some species can avoid the effects of PopS. The popS mutant was also significantly delayed in colonization of tomato stems compared to the wild type. Some AvrE-type effectors from gammaproteobacteria suppress salicylic acid (SA)-mediated plant defenses, suggesting that PopS, a betaproteobacterial ortholog, has a similar function. Indeed, the popS mutant induced significantly higher expression of tomato SA-triggered pathogenesis-related (PR) genes than the wild type. Further, pretreatment of roots with SA exacerbated the popS mutant virulence defect. Finally, the popS mutant had no colonization defect on SA-deficient NahG transgenic tomato plants. Together, these results indicate that this conserved effector suppresses SA-mediated defenses in tomato roots and stems, which are R. solanacearum’s natural infection sites. Interestingly, PopS did not trigger necrosis when heterologously expressed in Nicotiana leaf tissue, unlike the AvrE homolog DspE_Pcc from the necrotroph Pectobacterium carotovorum subsp. carotovorum. This is consistent with the differing pathogenesis modes of necrosis-causing gammaproteobacteria and biotrophic R. solanacearum.

The type III-secreted AvrE effector family is widely distributed in high-impact plant-pathogenic bacteria and is known to suppress plant defenses for virulence. We characterized the biology of PopS, the only AvrE homolog made by the bacterial wilt pathogen Ralstonia solanacearum. To our knowledge, this is the first study of R. solanacearum effector function in roots and stems, the natural infection sites of this pathogen. Unlike the functionally redundant R. solanacearum effectors studied to date, PopS is required for full virulence and wild-type colonization of two natural crop hosts. R. solanacearum is a biotrophic pathogen that causes a nonnecrotic wilt. Consistent with this, PopS suppressed plant defenses but did not elicit cell death, unlike AvrE homologs from necrosis-causing plant pathogens. We propose that AvrE family effectors have functionally diverged to adapt to the necrotic or nonnecrotic lifestyle of their respective pathogens.

Ralstonia solanacearum, a widespread bacterial plant pathogen in the post-genomic era

Ralstonia solanacearum is a soil-borne bacterium causing the widespread disease known as bacterial wilt. Ralstonia solanacearum is also the causal agent of Moko disease of banana and brown rot of potato. Since the last R. solanacearum pathogen profile was published 10 years ago, studies concerning this plant pathogen have taken a genomic and post-genomic direction. This was pioneered by the first sequenced and annotated genome for a major plant bacterial pathogen and followed by many more genomes in subsequent years. All molecular features studied now have a genomic flavour. In the future, this will help in connecting the classical field of pathology and diversity studies with the gene content of specific strains. In this review, we summarize the recent research on this bacterial pathogen, including strain classification, host range, pathogenicity determinants, regulation of virulence genes, type III effector repertoire, effector-triggered immunity, plant signalling in response to R. solanacearum, as well as a review of different new pathosystems.

TAXONOMY

Bacteria; Proteobacteria; β subdivision; Ralstonia group; genus Ralstonia.

DISEASE SYMPTOMS

Ralstonia solanacearum is the agent of bacterial wilt of plants, characterized by a sudden wilt of the whole plant. Typically, stem cross-sections will ooze a slimy bacterial exudate. In the case of Moko disease of banana and brown rot of potato, there is also visible bacterial colonization of banana fruit and potato tuber.

DISEASE CONTROL

As a soil-borne pathogen, infected fields can rarely be reused, even after rotation with nonhost plants. The disease is controlled by the use of resistant and tolerant plant cultivars. The prevention of spread of the disease has been achieved, in some instances, by the application of strict prophylactic sanitation practices.

USEFUL WEBSITES

Stock centre: International Centre for Microbial Resources-French Collection for Plant-associated Bacteria CIRM-CFBP, IRHS UMR 1345 INRA-ACO-UA, 42 rue Georges Morel, 49070 Beaucouzé Cedex, France, http://www.angers-nantes.inra.fr/cfbp/. Ralstonia Genome browser: https://iant.toulouse.inra.fr/R.solanacearum. GMI1000 insertion mutant library: https://iant.toulouse.inra.fr/R.solanacearumGMI1000/GenomicResources. MaGe Genome Browser: https://www.genoscope.cns.fr/agc/microscope/mage/viewer.php?

Comparative genomic and transcriptome analyses of pathotypes of Xanthomonas citri subsp. citri provide insights into mechanisms of bacterial virulence and host range

BACKGROUND

Citrus bacterial canker is a disease that has severe economic impact on citrus industries worldwide and is caused by a few species and pathotypes of Xanthomonas. X. citri subsp. citri strain 306 (XccA306) is a type A (Asiatic) strain with a wide host range, whereas its variant X. citri subsp. citri strain A(w)12879 (Xcaw12879, Wellington strain) is restricted to Mexican lime.

RESULTS

To characterize the mechanism for the differences in host range of XccA and Xcaw, the genome of Xcaw12879 that was completed recently was compared with XccA306 genome. Effectors xopAF and avrGf1 are present in Xcaw12879, but were absent in XccA306. AvrGf1 was shown previously for Xcaw to cause hypersensitive response in Duncan grapefruit. Mutation analysis of xopAF indicates that the gene contributes to Xcaw growth in Mexican lime but does not contribute to the limited host range of Xcaw. RNA-Seq analysis was conducted to compare the expression profiles of Xcaw12879 and XccA306 in Nutrient Broth (NB) medium and XVM2 medium, which induces hrp gene expression. Two hundred ninety two and 281 genes showed differential expression in XVM2 compared to in NB for XccA306 and Xcaw12879, respectively. Twenty-five type 3 secretion system genes were up-regulated in XVM2 for both XccA and Xcaw. Among the 4,370 common genes of Xcaw12879 compared to XccA306, 603 genes in NB and 450 genes in XVM2 conditions were differentially regulated. Xcaw12879 showed higher protease activity than XccA306 whereas Xcaw12879 showed lower pectate lyase activity in comparison to XccA306.

CONCLUSIONS

Comparative genomic analysis of XccA306 and Xcaw12879 identified strain specific genes. Our study indicated that AvrGf1 contributes to the host range limitation of Xcaw12879 whereas XopAF contributes to virulence. Transcriptome analyses of XccA306 and Xcaw12879 presented insights into the expression of the two closely related strains of X. citri subsp. citri. Virulence genes including genes encoding T3SS components and effectors are induced in XVM2 medium. Numerous genes with differential expression in Xcaw12879 and XccA306 were identified. This study provided the foundation to further characterize the mechanisms for virulence and host range of pathotypes of X. citri subsp. citri.

Constraints on genome dynamics revealed from gene distribution among the Ralstonia solanacearum species

Because it is suspected that gene content may partly explain host adaptation and ecology of pathogenic bacteria, it is important to study factors affecting genome composition and its evolution. While recent genomic advances have revealed extremely large pan-genomes for some bacterial species, it remains difficult to predict to what extent gene pool is accessible within or transferable between populations. As genomes bear imprints of the history of the organisms, gene distribution pattern analyses should provide insights into the forces and factors at play in the shaping and maintaining of bacterial genomes. In this study, we revisited the data obtained from a previous CGH microarrays analysis in order to assess the genomic plasticity of the R. solanacearum species complex. Gene distribution analyses demonstrated the remarkably dispersed genome of R. solanacearum with more than half of the genes being accessory. From the reconstruction of the ancestral genomes compositions, we were able to infer the number of gene gain and loss events along the phylogeny. Analyses of gene movement patterns reveal that factors associated with gene function, genomic localization and ecology delineate gene flow patterns. While the chromosome displayed lower rates of movement, the megaplasmid was clearly associated with hot-spots of gene gain and loss. Gene function was also confirmed to be an essential factor in gene gain and loss dynamics with significant differences in movement patterns between different COG categories. Finally, analyses of gene distribution highlighted possible highways of horizontal gene transfer. Due to sampling and design bias, we can only speculate on factors at play in this gene movement dynamic. Further studies examining precise conditions that favor gene transfer would provide invaluable insights in the fate of bacteria, species delineation and the emergence of successful pathogens.

The role of prophage in plant-pathogenic bacteria

A diverse set of phage lineages is associated with the bacterial plant-pathogen genomes sequenced to date. Analysis of 37 genomes revealed 5,169 potential genes (approximately 4.3 Mbp) of phage origin, and at least 50% had no function assigned or are nonessential to phage biology. Some phytopathogens have transcriptionally active prophage genes under conditions that mimic plant infection, suggesting an association between plant disease and prophage transcriptional modulation. The role of prophages within genomes for cell biology varies. For pathogens such as Pectobacterium, Pseudomonas, Ralstonia, and Streptomyces, involvement of prophage in disease symptoms has been demonstrated. In Xylella and Xanthomonas, prophage activity is associated with genome rearrangements and strain differentiation. For other pathogens, prophage roles are yet to be established. This review integrates available information in a unique interface ( http://propnav.esalq.usp.br ) that may be assessed to improve research in prophage biology and its association with genome evolution and pathogenicity.

The in planta transcriptome of Ralstonia solanacearum: conserved physiological and virulence strategies during bacterial wilt of tomato

Plant xylem fluid is considered a nutrient-poor environment, but the bacterial wilt pathogen Ralstonia solanacearum is well adapted to it, growing to 10⁸ to 10⁹ CFU/g tomato stem. To better understand how R. solanacearum succeeds in this habitat, we analyzed the transcriptomes of two phylogenetically distinct R. solanacearum strains that both wilt tomato, strains UW551 (phylotype II) and GMI1000 (phylotype I). We profiled bacterial gene expression at ~6 × 10⁸ CFU/ml in culture or in plant xylem during early tomato bacterial wilt pathogenesis. Despite phylogenetic differences, these two strains expressed their 3,477 common orthologous genes in generally similar patterns, with about 12% of their transcriptomes significantly altered in planta versus in rich medium. Several primary metabolic pathways were highly expressed during pathogenesis. These pathways included sucrose uptake and catabolism, and components of these pathways were encoded by genes in the scrABY cluster. A UW551 scrA mutant was significantly reduced in virulence on resistant and susceptible tomato as well as on potato and the epidemiologically important weed host Solanum dulcamara. Functional scrA contributed to pathogen competitive fitness during colonization of tomato xylem, which contained ~300 µM sucrose. scrA expression was induced by sucrose, but to a much greater degree by growth in planta. Unexpectedly, 45% of the genes directly regulated by HrpB, the transcriptional activator of the type 3 secretion system (T3SS), were upregulated in planta at high cell densities. This result modifies a regulatory model based on bacterial behavior in culture, where this key virulence factor is repressed at high cell densities. The active transcription of these genes in wilting plants suggests that T3SS has a biological role throughout the disease cycle.

Ralstonia solanacearum is a widespread plant pathogen that causes bacterial wilt disease. It inflicts serious crop losses on tropical farmers, with major economic and human consequences. It is also a model for the many destructive microbes that colonize the water-conducting plant xylem tissue, which is low in nutrients and oxygen. We extracted bacteria from infected tomato plants and globally identified the biological functions that R. solanacearum expresses during plant pathogenesis. This revealed the unexpected presence of sucrose in tomato xylem fluid and the pathogen’s dependence on host sucrose for virulence on tomato, potato, and the common weed bittersweet nightshade. Further, R. solanacearum was highly responsive to the plant environment, expressing several metabolic and virulence functions quite differently in the plant than in pure culture. These results reinforce the utility of studying pathogens in interaction with hosts and suggest that selecting for reduced sucrose levels could generate wilt-resistant crops.

Sequencing of K60, type strain of the major plant pathogen Ralstonia solanacearum

Ralstonia solanacearum is a widespread and destructive plant pathogen. We present the genome of the type strain, K60 (phylotype IIA, sequevar 7). Sequevar 7 strains cause ongoing tomato bacterial wilt outbreaks in the southeastern United States. K60 generally resembles R. solanacearum CFBP2957, a Caribbean tomato isolate, but has almost 360 unique genes.

The awr gene family encodes a novel class of Ralstonia solanacearum type III effectors displaying virulence and avirulence activities

We present here the characterization of a new gene family, awr, found in all sequenced Ralstonia solanacearum strains and in other bacterial pathogens. We demonstrate that the five paralogues in strain GMI1000 encode type III-secreted effectors and that deletion of all awr genes severely impairs its capacity to multiply in natural host plants. Complementation studies show that the AWR (alanine-tryptophan-arginine tryad) effectors display some functional redundancy, although AWR2 is the major contributor to virulence. In contrast, the strain devoid of all awr genes (Δawr1-5) exhibits enhanced pathogenicity on Arabidopsis plants. A gain-of-function approach expressing AWR in Pseudomonas syringae pv. tomato DC3000 proves that this is likely due to effector recognition, because AWR5 and AWR4 restrict growth of this bacterium in Arabidopsis. Transient overexpression of AWR in nonhost tobacco species caused macroscopic cell death to varying extents, which, in the case of AWR5, shows characteristics of a typical hypersensitive response. Our work demonstrates that AWR, which show no similarity to any protein with known function, can specify either virulence or avirulence in the interaction of R. solanacearum with its plant hosts.

Genome-enabled perspectives on the composition, evolution, and expression of virulence determinants in bacterial plant pathogens

Genome sequence analyses of bacterial plant pathogens are revealing important insights into the molecular determinants of pathogenicity and, through transcript characterization, responses to environmental conditions, evidence for small RNAs, and validation of uncharacterized genes. Genome comparison sheds further light on the processes impacting pathogen evolution and differences in gene repertoire among isolates contributing to niche specialization. Information derived from pathogen genome analysis is providing tools for use in diagnosis and interference with host-pathogen interactions for the purpose of disease control. However, the existing information infrastructure fails to adequately integrate the increasing numbers of sequence data sets, bioinformatic analyses, and experimental characterization, as required for effective systems-level analysis. Enhanced standardization of data formats at the point of publication is proposed as a possible solution.

Contrasting recombination patterns and demographic histories of the plant pathogen Ralstonia solanacearum inferred from MLSA

We used multilocus sequence analysis (MLSA) on a worldwide collection of the plant pathogenic Ralstonia solanacearum (Betaproteobacteria) to retrace its complex evolutionary history. Using genetic imprints left during R. solanacearum evolution, we were able to delineate distinct evolutionary complex displaying contrasting dynamics. Among the phylotypes already described (I, IIA, IIB, III, IV), eight groups of strains with distinct evolutionary patterns, named clades, were identified. From our recombination analysis, we identified 21 recombination events that occurred within and across these lineages. Although appearing the most divergent and ancestral phylotype, phylotype IV was inferred as a gene donor for the majority of the recombination events that we detected. Whereas this phylotype apparently fuelled the species diversity, ongoing diversification was mainly detected within phylotype I, IIA and III. These three groups presented a recent expanding population structure, a high level of homologous recombination and evidences of long-distance migrations. Factors such as adaptation to a specific host or intense trading of infected crops may have promoted this diversification. Whether R. solanacearum lineages will eventually evolve in distinct species remains an open question. The intensification of cropping and increase of geographical dispersion may favour situations of phylotype sympatry and promote higher exchange of key factors for host adaptation from their common genetic pool.

A chromosomal insertion toolbox for promoter probing, mutant complementation, and pathogenicity studies in Ralstonia solanacearum

We describe here the construction of a delivery system for stable and directed insertion of gene constructs in a permissive chromosomal site of the bacterial wilt pathogen Ralstonia solanacearum. The system consists of a collection of suicide vectors-the Ralstonia chromosome (pRC) series-that carry an integration element flanked by transcription terminators and two sequences of homology to the chromosome of strain GMI1000, where the integration element is inserted through a double recombination event. Unique restriction enzyme sites and a GATEWAY cassette enable cloning of any promoter::gene combination in the integration element. Variants endowed with different selectable antibiotic resistance genes and promoter::gene combinations are described. We show that the system can be readily used in GMI1000 and adapted to other R. solanacearum strains using an accessory plasmid. We prove that the pRC system can be employed to complement a deletion mutation with a single copy of the native gene, and to measure transcription of selected promoters in monocopy both in vitro and in planta. Finally, the system has been used to purify and study secretion type III effectors. These novel genetic tools will be particularly useful for the construction of recombinant bacteria that maintain inserted genes or reporter fusions in competitive situations (i.e., during plant infection).

Phylogeny and population structure of brown rot- and Moko disease-causing strains of Ralstonia solanacearum phylotype II

The ancient soilborne plant vascular pathogen Ralstonia solanacearum has evolved and adapted to cause severe damage in an unusually wide range of plants. In order to better describe and understand these adaptations, strains with very similar lifestyles and host specializations are grouped into ecotypes. We used comparative genomic hybridization (CGH) to investigate three particular ecotypes in the American phylotype II group: (i) brown rot strains from phylotypes IIB-1 and IIB-2, historically known as race 3 biovar 2 and clonal; (ii) new pathogenic variants from phylotype IIB-4NPB that lack pathogenicity for banana but can infect many other plant species; and (iii) Moko disease-causing strains from phylotypes IIB-3, IIB-4, and IIA-6, historically known as race 2, that cause wilt on banana, plantain, and Heliconia spp. We compared the genomes of 72 R. solanacearum strains, mainly from the three major ecotypes of phylotype II, using a newly developed pangenomic microarray to decipher their population structure and gain clues about the epidemiology of these ecotypes. Strain phylogeny and population structure were reconstructed. The results revealed a phylogeographic structure within brown rot strains, allowing us to distinguish European outbreak strains of Andean and African origins. The pangenomic CGH data also demonstrated that Moko ecotype IIB-4 is phylogenetically distinct from the emerging IIB-4NPB strains. These findings improved our understanding of the epidemiology of important ecotypes in phylotype II and will be useful for evolutionary analyses and the development of new DNA-based diagnostic tools.

Pathogenomics of the Ralstonia solanacearum species complex

Ralstonia solanacearum is a major phytopathogen that attacks many crops and other plants over a broad geographical range. The extensive genetic diversity of strains responsible for the various bacterial wilt diseases has in recent years led to the concept of an R. solanacearum species complex. Genome sequencing of more than 10 strains representative of the main phylogenetic groups has broadened our knowledge of the evolution and speciation of this pathogen and led to the identification of novel virulence-associated functions. Comparative genomic analyses are now opening the way for refined functional studies. The many molecular determinants involved in pathogenicity and host-range specificity are described, and we also summarize current understanding of their roles in pathogenesis and how their expression is tightly controlled by an intricate virulence regulatory network.

The Comprehensive Phytopathogen Genomics Resource: a web-based resource for data-mining plant pathogen genomes

The Comprehensive Phytopathogen Genomics Resource (CPGR) provides a web-based portal for plant pathologists and diagnosticians to view the genome and trancriptome sequence status of 806 bacterial, fungal, oomycete, nematode, viral and viroid plant pathogens. Tools are available to search and analyze annotated genome sequences of 74 bacterial, fungal and oomycete pathogens. Oomycete and fungal genomes are obtained directly from GenBank, whereas bacterial genome sequences are downloaded from the A Systematic Annotation Package (ASAP) database that provides curation of genomes using comparative approaches. Curated lists of bacterial genes relevant to pathogenicity and avirulence are also provided. The Plant Pathogen Transcript Assemblies Database provides annotated assemblies of the transcribed regions of 82 eukaryotic genomes from publicly available single pass Expressed Sequence Tags. Data-mining tools are provided along with tools to create candidate diagnostic markers, an emerging use for genomic sequence data in plant pathology. The Plant Pathogen Ribosomal DNA (rDNA) database is a resource for pathogens that lack genome or transcriptome data sets and contains 131 755 rDNA sequences from GenBank for 17 613 species identified as plant pathogens and related genera. Database URL: http://cpgr.plantbiology.msu.edu.