The Type VI Secretion System of&nbsp;Xanthomonas&nbsp;phaseoli pv. manihotis&nbsp;is involved in virulence and&nbsp;in&nbsp;vitro&nbsp;motility

doi:10.21203/rs.3.rs-48039/v2

Download PDF

Research article

The Type VI Secretion System of Xanthomonas phaseoli pv. manihotis is involved in virulence and in vitro motility

https://doi.org/10.21203/rs.3.rs-48039/v2

This work is licensed under a CC BY 4.0 License

Journal Publication

published 06 Jan, 2021

Read the published version in BMC Microbiology →

You are reading this older preprint version

Read the latest preprint version →

Background: The type VI protein secretion system (T6SS) is important in diverse cellular processes in Gram-negative bacteria, including interactions with other bacteria and with eukaryotic hosts. In this study we analyze the evolution of the T6SS in the genus Xanthomonas and evaluate its importance of the T6SS for virulence and in vitro motility in Xanthomonas phaseoli pv. manihotis (Xpm), the causal agent of bacterial blight in cassava (Manihot esculenta). We delineate the organization of the T6SS gene clusters in Xanthomonas and then characterize proteins of this secretion system in Xpm strain CIO151. Results: We describe the presence of three different clusters in the genus Xanthomonas that vary in their organization and degree of synteny between species. Using a gene knockout mutagenesis, we also found that vgrG and hcp are required for maximal aggressiveness of Xpm on cassava plants while clpV is important for both motility and maximal aggressiveness. Conclusion: We characterized the T6SS in 15 different strains in Xanthomonas and our phylogenetic analyses suggest that the T6SS might have been acquired by a very ancient event of horizontal gene transfer and maintained through evolution, hinting at their importance for the adaptation of Xanthomonas to their hosts. Finally, we demonstrated that the T6SS of Xpm is functional, and significantly contributes to virulence and motility. This is the first experimental study that demonstrates the role of the T6SS in the Xpm-cassava interaction and the T6SS organization in the genus Xanthomonas.

General Microbiology

Applied & Industrial Microbiology

type VI secretion

Hcp

Vgr

IcmF

ClpV

Xanthomonas

bacterial pathogenesis

A large number of Gram-negative bacteria use the type VI secretion system (T6SS) to transport proteins across the bacterial cell envelope. This versatile protein secretion system seems to be involved in a variety of cellular processes in bacteria, including antibacterial activity, biofilm formation and interactions with eukaryotic hosts. Thus, the T6SS may confer a competitive advantage in multi-species environments. The T6SS is responsible for antagonism towards potentially competing bacteria by direct injection of protein effectors in species such as Pseudomonas aeruginosa¹^,²^,, Salmonella typhimurium³ and Agrobacterium tumefaciens⁴. In Acidovorax citrulli⁵ and Burkholderia cenocepacia⁶ the T6SS has been implicated in biofilm formation as well.

In addition, the T6SS also participates in the interactions of pathogenic and commensal bacteria with their eukaryotic hosts. For example, in P. aeruginosa, two out of the three T6SS clusters are important in virulence against eukaryotic cells⁷^,⁸. This system is also involved in cell to cell signaling and communication. Vibrio cholerae uses the T6SS to induce changes in the host cellular behavior that reduce the population of other, potentially competing, bacteria⁹. In plant pathogens, such as Pantoea ananatis, the T6SS plays a key role in pathogenesis and bacterial competition¹⁰. Despite its importance, information about the functions of this system in plant pathogens remains scarce.

The T6SS injects diverse effector proteins into cells by contracting a spike-containing inner tube that perforates the membrane of target cells¹¹. The system is typically encoded by a set of fifteen to twenty genes. However, bacterial genomes encoding for T6SSs share a group of thirteen fundamental core genes¹². One of the most important proteins in this machinery is ClpV, an ATPase that forms a hexameric complex that provides the energy required for secreting T6SS substrates¹³^,¹⁴. Aditionally, the intracellular multiplication protein F (IcmF or TssM), is essential for the secretion of the haemolysis-corregulated protein (Hcp)^15–17. Hcp and VgrG (Valine-Glycine repeats G), are both effector proteins and important parts of the structural machinery of the T6SS in A. tumefaciens ¹⁸. Hcp and VgrG show structural homology to proteins found in the tail structures of bacteriophages, suggesting an evolutionary relationship between the T6SSs and the cell-puncturing machinery of family bacteriophages belonging to the Myoviridae^19,20.

The T6SS components are encoded in gene clusters that vary in organization and frequency. For example P. aeruginosa has three T6SS clusters⁸ while Burkholderia thailandesis has five different clusters²¹. Boyer and collaborators¹² found that Xanthomonas axonopodis and Xantohomonas campestris both have two complementary T6SS loci. Similarly, Xanthomonas euvesicatoria 85-10 (Xeu) has two T6SS loci or clusters with 15 conserved components²². The T6SS has been partially characterized, through bioinformatics, in Xeu, Xanthomonas vesicatoria strain 1111 (ATCC 35937: Xv) and Xanthomonas perforans strain 91-118 (Xp)²³. Mutants have been generated for genes vgrG and clpV in Xeu²². But, no change in virulence was reported for those mutant strains. Moreover, in Xanthomonas citri, the T6SSmediates resistance to Dictyostelium predation²⁴. Overall, the T6SS of members of the Xanthomonas genus need further characterization.

Xanthomonas phaseoli pv. manihotis (Xpm) is the causal agent of cassava bacterial blight, an economically important disease in Africa and South America, causing losses that may reach up to 100% after three cycles of cassava production²⁵. This local and systemic pathogen induces a wide combination of symptoms such as angular leaf spots, blight, wilting, dieback, gum exudation and vascular necrosis²⁵. Here, we report a bioinformatic study of the organization of the T6SS cluster in the genus Xanthomonas. We also report on the importance of this system for bacterial virulence and in vitro motility.

Determining the core components of the T6SS in Xanthomonads

A bioinformatic search for genes involved in the T6SS machinery was performed by selecting a group of genes of P. aeruginosa¹⁵, and performing a TBLASTN search with the BLOSUM62 matrix²⁶ against the genomes of Xanthomonas citri subsp. citri strain 306 (Xcc3), Xeu and Xpm (Additional file: Table S1). Additionally, genes identified in Xpm were confirmed as orthologs using reciprocal best hit. In order to identify all components of T6SS components in other xanthomonads, we used the same genes from P. aeruginosa and performed a BLASTP search with the BLOSUM62 matrix²⁶ against the genomes of Xanthomonas oryzae pv. oryzae (Xoo), Xanthomonas campestris pv. campestris (Xcac) and Xanthomonas albilineans (Xalb) (Additional file: Table S1). In both cases, a homolog was considered as significant if the BLAST e-value was <10^-20 and the amino acid identity was at least 30%. A subsequent search for orthologs with ORTHOMCL²⁷ confirmed the results^28,29. To detect signatures of Horizontal Gene Transfer, a search for genomic islands and insertion sequences was performed using Alien Hunter³⁰ and IS finder³¹, respectively. These results allowed the reconstruction of the T6SS clusters of Xcc3, Xeu and Xpm (Figure 1). The resulting T6SS gene clusters were used in BLAST searches against the other Xanthomonas genomes (listed in Table S1). In addition, conservation of the T6SS in Xpm was assessed by BLASTN searches default parameters²⁶ for the 65 Xpm strains reported by Bart and collaborators³².

Phylogenetic reconstruction

The orthologous gene sequences were aligned using MUSCLE 3.8³³. All core ortholog alignments were concatenated into a super matrix in Geneious³⁴ for phylogenetic reconstruction. The phylogenetic tree was constructed using a maximum likelihood approach in RAxML V.7.2.8³⁵ using the GTR+G+I model of nucleotide evolution and partitions per gene in the different clusters. Support values for phylogeny branches were estimated by means of 1000 replicates of bootstrap in RAxML³⁵.

Cellular localization, protein family and motif prediction for T6SS genes of Xpm

The cellular localization of each protein in the dataset was determined using PSORT³⁶, with options for Gram-negative bacteria. To determine transmembrane region and their orientation, TMpred was used³⁷. The databases InterPro³⁸, CDD³⁹, and ProDOM^40,41 were used for protein domain identification. The Pfam database⁴¹ was used to assign protein families. The presence of known protein motifs in the T6SS components was determined using Motif Finder⁴² and MOTIF Search in TRANSFAC databases⁴³.

Bacterial strains and growth conditions

The strains, plasmids and primer sequences used in this study are described in Tables S2 and S3. Escherichia coli strain DH5α was grown on LB at 37ºC and the Xpm strain CIO151 was grown on LPGA (5 gr yeast extract, 5 gr dextrose, 5 gr peptone and 15 gr agar per liter of distilled water) at 28ºC. For inoculation assays, Xpm cells were grown for two days in LPGA agar medium with the appropriate antibiotics. Cells were subsequently grown overnight in Phi broth (5 gr yeast extract, 5 gr dextrose, 5 gr casamino acids and 15 g agar per liter of distilled water) at 28ºC with the appropriate antibiotics. To adjust the cell suspensions to an OD_600nm of 0.2, overnight cultures were harvested by centrifugation at 14000 rpm for 2 min and re-suspended in 10mM MgCl₂ without antibiotics. In vitro growth of Xpm strains was measured in Phi broth with the appropriate antibiotics at 28ºC and shaking at 200 rpm for 38 hours. Three independent measurements were performed.

Swimming motility assays

Swimming motility of Xpm strains was evaluated by measuring the motility diameter of in LPGA medium with 0.3% agar. Swimming plates were inoculated with cultures grown overnight and adjusted to an OD_600nm of 0.3. They were subsequently incubated at 25ºC for 24 and 48 hours. Two independent experiments were performed with four replicates each.

Mutagenesis of vgrG, clpV, icmF and hcp

Homologous recombination by single crossing-over was used to generate mutants of these genes in Xpm. A fragment of nearly 400 bp from the central region of each gene was amplified by PCR. A 25 μl reaction was performed including 1X Buffer, 2 mM MgCl₂, 0.2 mM dNTPs, 0.2 mM of each primer, 10-50 ng of DNA of Xpm and 2U of Taq polymerase (Invitrogen Corp.). The amplification proceeded with an initial denaturation step of 5 min at 95°C, followed by 35 cycles of 45 seconds at 95°C, 45 seconds at annealing temperature (Additional file: Table S3) and an extension time of 1 min per Kb of expected product at 72°C. The obtained fragments were subcloned in pENTR^TM/ D-TOPO® (Invitrogen Corp., Grand Island, New York, USA). The clone was digested with the enzyme EcoRI and the resulting fragment was inserted in the suicide vector pAC3.1⁴⁴; and transformed into E. coli DH5α cells. The colonies were selected on LB plates + kanamycin (50 μg/ml) + chloramphenicol (25μg/ml). Insertion of the fragment in the vector was confirmed by PCR (Additional file: Table S3) and sequencing (Macrogen Inc. Korea). The insertion of each mutagenesis fragment in Xpm was performed by triparental mating and the resulting colonies were confirmed by PCR and sequencing (Macrogen Inc, Korea). For the PCR confirmation, extension and nested primers were designed (Additional file: Table S3).

In planta virulence assay

The virulence of the four generated mutants of Xpm was assayed on susceptible cassava plants (HMC-1 and MCOL2215). We employed two inoculation methodologies to determine differences between the pathogenicity of wild type and mutant strains: the first consisted on making 2mm diameter perforations on the leaflets, as previously described by Restrepo and collaborators⁴⁵. These perforations were inoculated with 10 μL of liquid bacterial suspension in 10 mM MgCl₂ at an OD_600nm of 0.2 (10⁸ CFU/ml). Three independent experiments were performed with five replicates each. To measure the lesion area, we used the ImageJ package⁴⁶ and a one-way analysis of variance (ANOVA) test was performed to validate statistically significant differences between wild type and mutant bacteria. The second inoculation method was the leaf clipping method⁴⁷ that consisted of cutting one to two cm from the tip of the leaves with scissors previously dipped in the inoculum at a concentration of 0.2 OD_600nm. Three independent measurements were performed. The onset of symptoms was monitored until day 15 post inoculation. Dilution plating was performed for each bacterial suspension in order to ensure the presence of the bacterium and to measure the concentration of the initial inoculum. Xpm CIO151ΔhrpX mutant, where secretion of effector proteins by the T3SS has been abrogated rendering the strain non-pathogenic, was used as a negative control.

Structural genes and phylogenetic tree of the T6SS in the genus Xanthomonas

The T6SS of P. aeruginosa has been well studied both at the bioinformatic and experimental levels, and its importance in pathogenicity and in the interaction with other bacteria has been demonstrated^1,15. We used the T6SS genes of Pseudomonas aeruginosa PAO1 as references to identify the T6SS genes in 44 out of the 60 evaluated genomes of the genus Xanthomonas. Our BLASTP and ORTHOMCL analyses (see methods) show that the reconstructed T6SS clusters in Xanthomonas contain between twelve and sixteen genes. For further analyses, we only considered the genomes of species that had at the minimum set of genes required for T6SS functionality⁴⁸. Therefore, we selected 14 representative genomes within the Xanthomonas genus (Additional file: Table S4 bold letter) to perform phylogenetic analyses that offer insights into the evolutionary paths and diversity in organization of the T6SS in this group of plant pathogens (Figure 1). We identified each of the genes encoded in the T6SS in these genomes by OrthoMCL and BLASTP (Additional file: Table S5), and detected a subset of the genes in additional genomes (Xanthomonas cassavae str. CFBP 4642, Xanthomonas perforans str 91-18 and Xanthomonas axonopodis str.29) using EDGAR2.0 ⁴⁹; Additional file: Table S6). These analyses suggest that the T6SS is present in a widespread array of species in Xanthomonas, pinpointing at the biological importance of this cluster in this genus.

Among the Xanthomonas species with putative functional T6SS, we found three different types of clusters (referred to I, II and II in Figure 1, additional file: Table S4, Table S5). Remarkably, the T6SS cluster are not the result of simple duplication events. For example, the Xoo strains have two clusters, one of which follows the phylogeny of the species, being distantly related form the phaseoli clade. However, the second cluster forms a monophyletic cluster with one of the copies from X. euvesicatoria and X. perforans. This implies that the clusters have a recent common origin and could have been horizontally acquired. We observed conserved synteny in the organization of Cluster III of Xv, Xeu, Xfa1 and Xcc3 (all members of axonopodis clade^28,29). This suggests that a common origin of the clusters precedes the divergence of the pathovars. In summary, there is evidence for both vertical and horizontal inheritance of the T6SS clusters in Xanthomonas.

In agreement with previous reports ²³^,²⁴, transcriptional regulators (LysR, TssB and TssA) were detected at the boundaries of some of the reconstructed T6SS clusters (I and III). The presence of noncoding RNAs in Cluster I suggests additional post-transcriptional regulation. Notably, we detected transcriptional regulators of the LysR family as a new class of regulators for Cluster I. The identified regulators may conditionally act individually or in combination to regulate their target clusters. No regulators were identified for type II clusters. Surprisingly, transcriptional regulators of the AraC family, previously propopsed as a characteristic feature of type III clusters²³, were not detected in the clusters of Xfa and Xcc.

Pseudogenes ⁵⁰ were noted in the T6SS clusters of some of the genomes analyzed (Figure 1, additional file: Table S5). We considered some proteins with an early stop codon that were divided in two parts as well as other genes that had shorter versions or frameshifts (Figure 1). fha, clpV and impH were among the truncated genes in the clusters I of Xvm0 and XvmN. Since Fha is a target for TagF involved in posttranslational regulation of P. aeruginosa and A. tumefaciens⁵¹, this process may be affected in these strains and they might have an alternative type of regulation, which would need to be experimentally determined. Other divided genes of the T6SS clusters include: icmF in Xvm0 and XvmN (cluster I), ompA in Xvm0 (cluster I), and clpV in Xpm (cluster III) (Figure 1). Experimental procedures are necessary to determine the function of these genes in the Xanthomonas T6SS. Sixteen Xanthomonas strains did not contain the core genes of the T6SS; some were Xanthomonas campestris pv. campestris (strain 8004, strain ATCC33913 and strain B100), Xanthomonas campestris pv. armoraciae str. 756C and Xanthomonas albilineans and were therefore considered as T6SS-depleted.

Xpm, Xeu and Xcc3 show similarity in the organization of the T6SS gene clusters

A comparative analysis was performed between Xpm and its closest relatives with completely sequenced genomes, Xeu and Xcc3²⁸. Figure 2 shows the presence of two T6SS clusters (containing the same genes in different order) of Xeu, In contrast, Xpm and Xcc3 only have a single T6SS cluster (cluster III), orthologous to cluster III of Xeu.

To determine if the T6SS is conserved among members of the Xpm pathovar, or is exclusively present in the CIO151 strain, we searched for the T6SS genes in the genomes available for this pathovar³². We performed a BLASTN analyses with each T6SS gene of CIO151 against the other 64 genomes of manihotis pathovar (data not shown). We found the 16 T6SS genes are conserved among the analyzed genomes, with 75% to 100% identity at the nucleotide level. We found that clpV is divided in two parts due to a stop codon in the middle of the two AAA-ATPase domains. Nonetheless, because the stop codon does not truncate any of the protein domains (Additional file: Table S7), the two ClpV fragments may interact to create a functional protein complex (Figure 3 and 4).

The cluster of orthologous groups (COG), together with the results obtained from our phylogenetic results discussed above, suggest Horizontal Gene Transfer (HGT) events as the probable origin of the T6SS in Xanthomonas. Thus we searched for genomic islands (GI), insertion sequences (IS), deviations in the GC content, and tRNA genes surrounding the T6SS clusters¹². The positions and types of identified GI and IS were not exactly the same in clusters I and clusters III of the three organisms (Xpm, Xeu and Xcc3, Figure 2 and additional file: Table S8). However, GIs ere consistently located between the vgrG and clpV genes. The IS within the T6SS clusters of Xeu and Xpm belong to the IS3 family, while the IS in the Xcc3 cluster is part of the IS4 family³¹ (Figure 2 and additional file: Table S8). In contrast, no significant difference in GC content (with respect to the average GC content of the corresponding genomes) was detected int the T6SS cluster of Xeu, Xcc3, and Xpm. Hence, there are some characteristics of HGT for these clusters, but if these events occurred, there has been enough time to adapt these clusters to the rest of the genome.

Mutations in vgrG, clpV and hcp decrease the virulence of Xpm on susceptible cassava plants

To determine whether the T6SS contributes to the virulence of Xpm, we constructed and then inoculated vgrG, clpV, icmF and hcp single deletion mutants into susceptible cassava plants. Initially, in vitro growth assays were performed for each mutant strain, to determine whether the introduced mutations decrease the general fitness of the pathogen (Figure 3). The mutants were complemented with the wild type genes in all cases to restore the genotype. The resulting vgrG, clpV and hcp mutant strains did not show significant differences in fitness with respect to the wild type control (paired t-test p-values=0.14, 0.48 and 0.15, respectively). icmF mutants showed a small decrease in growth in vitro at 18 hours with respect to the wild type strain (paired t-test p-value=0.02) but the effect was undetectable at later time points.

These results suggest that while mutations in vgrG, clpV and hcp do not affect the ability of Xpm to grow in vitro, mutation in icmF has a slight effect on the growth of this pathogen.

We then tested the mutants for their ability to cause disease in susceptible cassava plants. The vgrG, clpV and hcp mutant strains were able to produce symptoms on leaves of susceptible cassava. Therefore, we conclude that these genes are not required for full pathogenicity (Figure 4). However, a decrease in symptoms was observed for CIO151ΔvgrG, CIO151ΔclpV, and CIO151Δhcp strains at 15 days post-inoculation by two different methods of inoculation (Fig. 4a and Fig. 4b). For the hole-inoculation method, the lesion area was measured with the program ImageJ⁴⁶ and statistically significant differences when compared to the wild type were observed (ANOVA p-value>0.01), as show in Figure 4c. The phenotype was complemented when mutants were transformed back with their respective wild type gene (Fig. 4c). Together, these results suggest that the genes vgrG, clpV and hcp are required for maximal aggressiveness of Xpm on susceptible cassava plants.

clpV deletion decreases Xpm motility

The T6SS has been implicated in pathogenicity^15,52,53, motility^54–56 and interaction with other bacteria^21,57. Motility contributes to virulence in the genus Xanthomonas^58–60. We therefore tested the T6SS mutants for swimming motility on petri dish with a low proportion of agar (0,3% agar) after 24 h and 48 h. Motility was significantly different after 48 for the CIO151ΔclpV mutant with respect to the (p-value=0.0007). Notably, the motility of this strain was indistinguishable from that of the CIO151 strain, after transformation with the wild type clpV gene (Fig. 5c), which demonstrates that the observed differences was due to the truncation/deletion of this gene. For the other evaluated mutants, there were no differences in motility with respect to the wild type (p-value>0.05; Fig. 5a, 5b and 5d). In general, these results suggest that ClpV activity impacts motility, whereas VgrG, IcmF and Hcp are dispensable for this function in Xpm.

The T6SS is a versatile macromolecular assembly that has not been comprehensively characterized in plant pathogens. We have provided insights into the gene cluster organization and phylogeny of 15 genes of the T6SS machinery in the genus Xanthomonas. We found that the T6SS may have suffered duplications and HGT events in Xanthomonas. Finally, we show the importance of the T6SS for aggressiveness in susceptible cassava plants and in vitro motility of Xpm.

We detected the presence of three distinct T6SS clusters in the genus Xanthomonas that differ in the order in which the orthologous genes are arranged (Figure 1 and Figure S1). The number of T6SS clusters has also been found to be variable in other taxa¹². The organization of the clusters I and III is similar, containing 16 genes, contrary to cluster II, with only 12 genes. Because the T6SS cluster II always co-occurs with one of the other clusters, cluster II may act as a complementary system in Xoo. This has previously been reported for other plant pathogens such as B. glumae, B gladioli and B. plantari strains, which generally possess at least two potential functional T6SS clusters whilst most representatives of beneficial (plant growth-promoting, symbionts or nodule-forming) Burkholderia species showed one or two clusters⁶¹.

The presence of two T6SS in X. oryzae, Xeu and Xp could indicate that both clusters together are required to function properly. This is also the case of P. aeruginosa, where the cluster I, which is involved in virulence¹⁵ and toxicity¹ requires the activity of clusters II and III⁸. Similarly, Vibrio cholerae has two auxiliary clusters with important activity during competition between V. cholerae strains⁶². Hence, it would be interesting to define the environments where individual T6SS clusters are active and conditionally essential in Xanthomonads.

Several lines of evidence suggested horizontal gene transfer events in the T6SS regions of the analyzed Xanthomonas. This observation agrees with previous phylogenetic studies, in other taxa, reporting the presence of T6SS gene clusters in horizontally acquired pathogenicity islands^12,65. The presence of insertion sequences and genomic islands in the T6SS cluster III, together with the disagreement between the genomes-based and T6SS-based phylogenetic relationships of the genus Xanthomonas²⁸ support the independent acquisition of the cluster III, by HGT. No differences were found in GC content were detected between the T6SS clusters and their corresponding genomes. This pattern of homogeneity in GC content was previously found in the T6SS of Klebsiella spp. and V. cholerae^63,64. Thus, we hypothesize that if the horizontal acquisition of the T6SS indeed happened in Xanthomonas, the horizontally acquired T6SS regions have adapted to the characteristics (GC content, codon usage, etc.) of the receiving genome to optimize its expression.

T6SS gene expression is tightly regulated at the transcriptional [by transcription factors (TFs)] and post-transcriptional levels (e.g. phosporylation). Cluster I showed a phosphorylation-type regulator (Kinase / Phosphatase / Forkhead). Cluster III has both a phosphorylation-type regulator and an AraC-type TF. Cluster II does not have either cis-acting regulator. Miyata et al. recently proposed that the presence of multiple regulators may help a pathogen coordinate T6SS gene expression to avoid identification by the host immune system⁶⁵. A similar phenomenon could be occurring in Xanthomonads. Understanding the regulatory mechanisms of the T6SS in Xanthomonas will be a key step to uncover its role in diverse ecological interactions.

A single T6SS (with 16 genes and classified as type III) was identified in the Xpm CIO151 genome. We showed that the T6SS is an important factor for virulence and motility in this pathovar. clpV is divided in two contigous sequences in both Xp and Xpm. The importance of ClpV is species-specific. For example, ClpV supplies the energy for the assembly of the T6SS external machinery, composed of VipA (ImpB) and VipB (ImpD) in V. cholerae and P. aeruginosa species^14,66. However, the lack of ClpV does not affect the functionality of the T6SS of Campylobacter jejuni⁶⁷. In Xpm CIO151, each clpV fragment contain a potentially functional ATPase domain and P-loop containing nucleoside triphosphate hydrolase, corresponding to a common ClpV protein. Using a ΔclpV mutant, we demonstrated that clpV is important for in vitro motility of Xpm. This finding suggests that the fragmented clpV gene is still functional. Motility plays a role predominantly in the early phases of infection; thus it helps in the development of plant disease in Xanthomonas⁶⁸. Xpm probably requires motility for the pre-entry processes and the spread of the pathogen inside the plant, but no experimental tests have been performed to test this. To fully understand the role of clpV in motility, it would be important to fully characterize the regulation of motility in Xpm and what molecular role ClpV plays in it.

We demonstrate a decrease in virulence for the CIO151ΔvgrG, CIO151ΔclpV and CIO151Δhcp mutants. Hcp has been reported as a substrate of the T6SS in P. aeruginosa¹⁵and Burkholderia mallei²¹. More importantly, Hcp was required for full tumorigenesis efficiency in A. tumefaciens¹⁶, in agreement with our results for Xpm. Mougous and collaborators (2006)¹⁵ demonstrated that ClpV is necessary for the secretion of Hcp in P. aeruginosa. In Xpm, we have demonstrated that ClpV has pleiotropic effects on two processes – motility and virulence. It is therefore possible that Hcp is a substrate of ClpV, affecting other aspects of virulence different from motility, and that other ClpV-dependent, Hcp- and VgrG-independent activities are required for motility in Xpm. In addition to Hcp- and VgrG, a few T6SS effectors have been reported until now, namely EvpP from Edwardsiella tarda⁶⁹, RbsB in Rhizobium leguminosarum^70,71, TssM in B. mallei²¹, and Tse1, Tse2, and Tse3 in P. aeruginosa^1,72. Recently, Bayer-Santos and collaborators (2019) in silico predicted a series of T6SS effectors from Xanthomonadales using an in silico analysis and suggested the presence of a high number of putative antibacterial toxins⁷³ which need to be experimentally tested. Bioinformatic analyses (using the Bastion6 software⁷⁴) could be used to identify effector candidates for the T6SS in Xpm. Those predictions could be then experimentally evaluated in follow up experiments to better elucidate the distinct roles of each effector in the phenotypes observed here.

Our computational analyses identify 16 proteins of the T6SS in the genus Xanthomonas. T6SS presents three different T6SS-associated gene clusters are present in the genus Xanthomonas that vary principally in the organization and the synteny of orthologous genes between species. Clusters III and I have the same number of genes and organization, while cluster II only has 12 genes, and it is restricted to X. oryzae strains. Phylogenetic analyses suggest that the T6SS might have been acquired by a very ancient event of horizontal gene transfer and maintained through evolution, hinting at their importance for the adaptation of Xanthomonas to their hosts. Finally, we demonstrated that the T6SS of Xpm is functional, and significantly contributes to virulence and motility.

T6SS: type VI secretion system; Xpm: Xanthomonas phaseoli pv. manihotis; IcmF: Intracellular multiplication protein F; Hcp: Haemolysis-corregulated protein; VgrG: Valine-Glycine repeats G; COG: Cluster of orthologous groups; HGT: Horizontal Gene Transfer; IS: Insertion Sequence: Xcc3: Xanthomonas citri subsp. citri strain 306; Xeu: Xanthomonas euvesicatoria 85-10; Xv: Xanthomonas vesicatoria strain 1111; Xp: Xanthomonas perforans strain 91-118; Xoo: Xanthomonas oryzae pv. oryzae; Xcac: Xanthomonas campestris pv. campestris; Xalb: Xanthomonas albilineans.

Acknowledgments

We thank the College of Science of the Universidad de los Andes for providing the infrastructure and all the necessary support for carrying out this work.

Funding

This work was funded by Facultad de Ciencias, Universidad de los Andes (INV-2019-84-1855) and Departamento Administrativo de Ciencia, Tecnología e Innovación (521-2010), Colombia. The funders had no role in the study design, data collection, analysis and interpretation of the data, decision to publish, or preparation of the manuscript.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Authors’ contributions

NMB performed all laboratory experiments, analyzed data and wrote manuscript. AAB obtained the vgrG mutant and helped draft the manuscript. LMR, MLAO, JFT, DB performed bioinformatic analysis and helped draft the manuscript. LC and CT helped design study and helped draft manuscript. AB contributed to experimental design, helped in data analysis and helped write the manuscript. SR help draft the manuscript. All authors have read and approve the manuscript.

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Hood, R. et al. System of Pseudomonas aeruginosa Targets a Toxin to Bacteria. Cell Host Microbe 7, 25–31 (2010).
Jiang, F., Waterfield, N. R., Yang, J., Yang, G. & Jin, Q. A Pseudomonas aeruginosa type VI secretion phospholipase D effector targets both prokaryotic and eukaryotic cells. Cell Host Microbe 15, 600–610 (2014).
Sana, T. G. et al. Salmonella Typhimurium utilizes a T6SS-mediated antibacterial weapon to establish in the host gut. Proc. Natl. Acad. Sci. U. S. A. 113, E5044-51 (2016).
Ma, L.-S., Hachani, A., Lin, J.-S., Filloux, A. & Lai, E.-M. Agrobacterium tumefaciens deploys a superfamily of type VI secretion DNase effectors as weapons for interbacterial competition in planta. Cell Host Microbe 16, 94–104 (2014).
Tian, Y. Y. et al. The type VI protein secretion system contributes to biofilm formation and seedtoseedling transmission of Acidovorax citrulli on melon. Mol. Plant Pathol. 16, 38–47 (2014).
Aubert, D. F., Flannagan, R. S. & Valvano, M. A. A novel sensor kinase-response regulator hybrid controls biofilm formation and type VI secretion system activity in Burkholderia cenocepacia. Infect. Immun. 76, 1979–1991 (2008).
Sana, T. G. et al. The second type VI secretion system of Pseudomonas aeruginosa strain PAO1 is regulated by quorum sensing and Fur and modulates internalization in epithelial cells. J. Biol. Chem. 287, 27095–27105 (2012).
Lesic, B., Starkey, M., He, J., Hazan, R. & Rahme, L. G. Quorum sensing differentially regulates Pseudomonas aeruginosa type VI secretion locus I and homologous loci II and III, which are required for pathogenesis. Microbiology 155, 2845–2855 (2009).
Logan, S. L. et al. The Vibrio cholerae type VI secretion system can modulate host intestinal mechanics to displace gut bacterial symbionts. Proc. Natl. Acad. Sci. U. S. A. 115, E3779–E3787 (2018).
Shyntum, D. Y. et al. Pantoea ananatis Utilizes a Type VI Secretion System for Pathogenesis and Bacterial Competition. Mol. Plant. Microbe. Interact. 28, 420–431 (2015).
Cianfanelli, F. R., Monlezun, L. & Coulthurst, S. J. Aim, Load, Fire: The Type VI Secretion System, a Bacterial Nanoweapon. Trends in Microbiology vol. 24 51–62 (2016).
Boyer, F., Fichant, G., Berthod, J., Vandenbrouck, Y. & Attree, I. Dissecting the bacterial type VI secretion system by a genome wide in silico analysis: what can be learned from available microbial genomic resources. BMC Genomics 10, 1115–1127 (2009).
Schlieker, C., Zentgraf, H., Dersch, P. & Mogk, A. ClpV, a unique Hsp100/Clp member of pathogenic proteobacteria. Biol. Chem. 386, 1115–1127 (2005).
Bönemann, G., Pietrosiuk, A., Diemand, A., Zentgraf, H. & Mogk, A. Remodelling of VipA/VipB tubules by ClpV‐mediated threading is crucial for type VI protein secretion. EMBO J. 28, 315–325 (2009).
Mougous, J. D. et al. A virulence locus of Pseudomonas aeruginosa encodes a protein secretion apparatus. Science 312, 1526–1530 (2006).
Wu, H.-Y., Chung, P.-C., Shih, H.-W., Wen, S.-R. & Lai, E.-M. Secretome analysis uncovers an Hcp-family protein secreted via a type VI secretion system in Agrobacterium tumefaciens. J. Bacteriol. 190, 2841–2850 (2008).
Mattinen, L. et al. Microarray profiling of host-extract-induced genes and characterization of the type VI secretion cluster in the potato pathogen Pectobacterium atrosepticum. Microbiology 154, 2387–2396 (2008).
Salomon, D. et al. Marker for type VI secretion system effectors. Proc. Natl. Acad. Sci. U. S. A. 111, 9271–9276 (2014).
Douzi, B. et al. Crystal structure and self-interaction of the type VI secretion tail-tube protein from enteroaggregative Escherichia coli. PLoS One 9, e86918 (2014).
Pukatzki, S., Ma, A. T., Revel, A. T., Sturtevant, D. & Mekalanos, J. J. Type VI secretion system translocates a phage tail spike-like protein into target cells where it cross-links actin. Proc. Natl. Acad. Sci. U. S. A. 104, 15508–15513 (2007).
Schwarz, S. et al. Burkholderia Type VI Secretion Systems Have Distinct Roles in Eukaryotic and Bacterial Cell Interactions. PLOS Pathog. 6, e1001068 (2010).
Abendroth, U. et al. Identification of new protein-coding genes with a potential role in the virulence of the plant pathogen Xanthomonas euvesicatoria. BMC Genomics 18, (2017).
Potnis, N. et al. Comparative genomics reveals diversity among xanthomonads infecting tomato and pepper. BMC Genomics 12, 146 (2011).
Bayer-Santos, E. et al. Xanthomonas citri T6SS mediates resistance to Dictyostelium predation and is regulated by an ECF σ factor and cognate Ser/Thr kinase. Environ. Microbiol. (2018) doi:10.1111/1462-2920.14085.
Lozano, E. A. Etiologic significance of bacterial isolates from rams with palpable epididymitis. Am. J. Vet. Res. 47, 1153–1156 (1986).
Altschul, S. F. et al. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 25, 3389–3402 (1997).
Li, L., Stoeckert, C. J. & Roos, D. S. OrthoMCL: Identification of ortholog groups for eukaryotic genomes. Genome Res. 13, 2178–2189 (2003).
Rodriguez-R, L. M. et al. Genomes-based phylogeny of the genus Xanthomonas. BMC Microbiol. 12, 43 (2012).
Arrieta-Ortiz, M. L. et al. Genomic survey of pathogenicity determinants and VNTR markers in the cassava bacterial pathogen Xanthomonas axonopodis pv. manihotis strain CIO151. PLoS One 8, (2013).
Vernikos, G. S. & Parkhill, J. Interpolated variable order motifs for identification of horizontally acquired DNA: Revisiting the Salmonella pathogenicity islands. Bioinformatics 22, 2196–2203 (2006).
Siguier, P., Perochon, J., Lestrade, L., Mahillon, J. & Chandler, M. ISfinder: the reference centre for bacterial insertion sequences. Nucleic Acids Res. 34, D32-6 (2006).
Bart, R. et al. High-throughput genomic sequencing of cassava bacterial blight strains identifies conserved effectors to target for durable resistance. Proc. Natl. Acad. Sci. U. S. A. 109, E1972-9 (2012).
Edgar, R. C. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 32, 1792–1797 (2004).
Kearse, M. et al. Geneious Basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics 28, 1647–1649 (2012).
Stamatakis, A. RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics 22, 2688–2690 (2006).
Yu, N. Y. et al. PSORTb 3.0: Improved protein subcellular localization prediction with refined localization subcategories and predictive capabilities for all prokaryotes. Bioinformatics (2010) doi:10.1093/bioinformatics/btq249.
Hofmann, K., Stoffel, W., Hoffman, K. & Stoffel, W. TMBASE - A database of membrane spanning protein segments. Biol Chem HoppeSeyler Volume p 166 vols. 374 SRC- (1993).
Zdobnov, E. M. & Apweiler, R. InterProScan--an integration platform for the signature-recognition methods in InterPro. Bioinformatics 17, 847–848 (2001).
Marchler-Bauer, A. et al. CDD: a Conserved Domain Database for protein classification. Nucleic Acids Res. 33, D192-6 (2005).
Corpet, F., Gouzy, J. & Kahn, D. Recent improvements of the ProDom database of protein domain families. Nucleic Acids Res. 27, 263–267 (1999).
Sonnhammer, E. L., Eddy, S. R. & Durbin, R. Pfam: A comprehensive database of protein domain families based on seed alignments. Proteins 7 December pp 28, 405–420 (1998).
An, J., Nakama, T., Kubota, Y. & Sarai, A. 3DinSight: an integrated relational database and search tool for the structure, function and properties of biomolecules. Bioinformatics 14, 188–195 (1998).
Wingender, E. et al. TRANSFAC: an integrated system for gene expression regulation. Nucleic Acids Res. 28, 316–319 (2000).
Castañeda, A., Reddy, J. D., El-Yacoubi, B. & Gabriel, D. W. Mutagenesis of all eight avr genes in Xanthomonas campestris pv. campestris had no detected effect on pathogenicity, but one avr gene affected race specificity. Mol. Plant. Microbe. Interact. 18, 1306–1317 (2005).
Restrepo, S., Duque, M. C. & Verdier, V. Characterization of pathotypes among isolates of Xanthomonas axonopodis pv. manihotis in Colombia. Plant Pathol. 49, 680–687 (2000).
Schneider, C. A. et al. NIH Image to ImageJ: 25 years of image analysis. Nat. Methods 9, 671–675 (2012).
Kauffman, H. E. An improved technique for evaluat-ing resistance of rice varieties to Xanthomonas oryzae. Plant Dis Rep 57, 537–541 (1973).
Shrivastava, S. & Mande, S. S. Identification and functional characterization of gene components of type VI secretion system in bacterial genomes. PLoS One 3, (2008).
Blom, J. et al. EDGAR 2.0: an enhanced software platform for comparative gene content analyses. Nucleic Acids Res. 44, W22–W28 (2016).
Tutar, Y. Pseudogenes. Comp. Funct. Genomics 2012, 424526 (2012).
Lin, J.-S. et al. TagF-mediated repression of bacterial type VI secretion systems involves a direct interaction with the cytoplasmic protein Fha. J. Biol. Chem. 293, 8829–8842 (2018).
Schell, M. A. et al. Type VI secretion is a major virulence determinant in Burkholderia mallei. Mol. Microbiol. 64, 1466–1485 (2007).
Pukatzki, S. et al. Identification of a conserved bacterial protein secretion system in Vibrio cholerae using the Dictyostelium host model system. Proc. Natl. Acad. Sci. 103, 1528–1533 (2006).
Records, A. R. & Gross, D. C. Sensor kinases RetS and LadS regulate Pseudomonas syringae Type VI secretion and virulence factors. J. Bacteriol. 192, 3584–3596 (2010).
Burtnick, M. N., DeShazer, D., Nair, V., Gherardini, F. C. & Brett, P. J. Burkholderia mallei Cluster 1 Type VI Secretion Mutants Exhibit Growth and Actin Polymerization Defects in RAW 264.7 Murine Macrophages. Infect. Immun. 78, 88 LP – 99 (2010).
Liu, L. et al. The type VI secretion system modulates flagellar gene expression and secretion in <em>Citrobacter freundii</em> and contributes to the adhesion and cytotoxicity to host cells. Infect. Immun. (2015).
Hood, R. D. et al. A Type VI Secretion System of Pseudomonas aeruginosa Targets a Toxin to Bacteria. Cell Host Microbe 7, 25–37 (2010).
Kamoun, S. & Kado, C. I. Phenotypic Switching Affecting Chemotaxis, Xanthan Production, and Virulence in Xanthomonas campestris. Appl. Environ. Microbiol. 56, 3855–3860 (1990).
Ryan, R. P. et al. Cyclic di-GMP signalling in the virulence and environmental adaptation of Xanthomonas campestris. Mol. Microbiol. 63, 429–442 (2007).
Ryan, R. P. et al. Pathogenomics of Xanthomonas: understanding bacterium-plant interactions. Nat. Rev. Microbiol. 9, 344–355 (2011).
Nguyen, T. T., Lee, H.-H., Park, I. & Seo, Y.-S. Genome-Wide Analysis of Type VI System Clusters and Effectors in Burkholderia Species. plant Pathol. J. 34, 11–22 (2018).
Kirchberger, P. C., Unterweger, D., Provenzano, D., Pukatzki, S. & Boucher, Y. Sequential displacement of Type VI Secretion System effector genes leads to evolution of diverse immunity gene arrays in Vibrio cholerae. Sci. Rep. 7, 45133 (2017).
Das, S. & Chaudhuri, K. Identification of a unique IAHP (IcmF associated homologous proteins) cluster in Vibrio cholerae and other proteobacteria through in silico analysis. In Silico Biol. 3, 287–300 (2003).
Sarris, P. F., Zoumadakis, C., Panopoulos, N. J. & Scoulica, E. V. Distribution of the putative type VI secretion system core genes in Klebsiella spp. Infect. Genet. Evol. 11, 157–166 (2011).
Miyata, S. T., Bachmann, V. & Pukatzki, S. Type VI secretion system regulation as a consequence of evolutionary pressure. J. Med. Microbiol. 62, 663–676 (2013).
Basler, M., Ho, B. T. & Mekalanos, J. J. Tit-for-tat: Type VI secretion system counterattack during bacterial cell-cell interactions. Cell 152, 884–894 (2013).
Lertpiriyapong, K. et al. Campylobacter jejuni Type VI Secretion System: Roles in Adaptation to Deoxycholic Acid, Host Cell Adherence, Invasion, and In Vivo Colonization. PLoS One 7, e42842 (2012).
Sena-Vélez, M. et al. Biofilm formation and motility of Xanthomonas strains with different citrus host range. Plant Pathol. 64, (2014).
Zheng, J. & Leung, K. Y. Dissection of a type VI secretion system in Edwardsiella tarda. Mol. Microbiol. 66, 1192–1206 (2007).
Bladergroen, M. R., Badelt, K. & Spaink, H. P. Infection-Blocking Genes of a Symbiotic Rhizobium leguminosarum Strain That Are Involved in Temperature-Dependent Protein Secretion. Mol. Plant-Microbe Interact. 16, 53–64 (2003).
James, D., Shao, H., Lamont, R. J. & Demuth, D. R. The Actinobacillus actinomycetemcomitans Ribose Binding Protein RbsB Interacts with Cognate and Heterologous Autoinducer 2 Signals. Infect. Immun. 74, 4021 LP – 4029 (2006).
Russell, A. B. et al. Type VI secretion delivers bacteriolytic effectors to target cells. Nature 475, 343–349 (2011).
Bayer-Santos, E., Ceseti, L. de M., Farah, C. S. & Alvarez-Martinez, C. E. Distribution, Function and Regulation of Type 6 Secretion Systems of Xanthomonadales . Frontiers in Microbiology vol. 10 1635 (2019).
Wang, J. et al. Bastion6: a bioinformatics approach for accurate prediction of type VI secreted effectors. Bioinformatics 34, 2546–2555 (2018).

Figure S1. Genomic organization of the characterized T6SS cluster. Alignment orthologous sequences and the hits of T6SS genes are represented as colored arrows. The colors follow the same guidelines as in Fig. 1. A unique color was assigned to highly conserved orthologs of Xanthomonas euvesicatoria str. 85-10, Xanthomonas citri subsp. citri str. 306 and Xanthomonas phaseoli pv. manihotis str. CIO 151. The lines represent the orthologous genes found with ORTHOMCL.

Additional file:

Table S1. TBLASTN and BLASTP of homologous structural genes of Pseudomonas aeruginosa T6SS selected. Homologous structural genes of Pseudomonas aeruginosa T6SS selected.by TBLASTN (Altschul et al., 1997) against the genomes of Xanthomonas citri subsp. citri str. 306 (Xcc), Xanthomonas euvesicatoria str. 85-10 (Xeu) and Xanthomonas phaseoli pv. manihotis str. CIO151 (Xpm). Similarity of structural genes of Pseudomonas aeruginosa T6SS selected by BLASTP against the genomes of Xanthomonas oryzae pv. oryzae str. PXO 099 (Xoo), Xanthomonas campestris pv. campestris str. ATCC 33913 (Xcac) and Xanthomonas albilineans GPE PC73 (Xalb). The colors follow the same guidelines as in Fig. 1.

Table S2. Bacterial strains and plasmids. Strains of Xanthomonas phaseoli pv. manihotis str. CIO151 (Xpm) tested for the presence of T6SS gene clusters.

Table S3: List of primers. List of primers used for PCR amplification of type VI secretion system genes from different strains of Xanthomonas phaseoli pv. manihotis str. CIO151 (Xpm).

Table S4. Xanthomonas genomes analysis with or without T6SS genes. Results of BLASTP and ORTHOMCL analyses show that the reconstructed T6SS clusters in Xanthomonas contain between twelve and sixteen genes.

Table S5: Structural genes of type VI secretion system of Xanthomonas.

Distribution of type VI secretion system gene clusters in 14 strains of Xanthomonas. For each strain of Xanthomonas, the genes found in T6SS divided for clusters (I, II and II) and are shown with gene, protein and product sizes.

Table S6: Results of the search for T6SS-encoding genes in the genomes of Xanthomonas perforans (Xp 91--18), Xanthomonas cassavae str. CFBP4642 and Xanthomonas axonopodis str. 29 with EDGAR 2.0. EDGAR displays the orthologous genes in each genome. For each gene, the nomenclature and annotation are denoted.

Table S7. Localization, protein family and motif prediction in Xpm. The proteins core of T6SS presents different localization in the cell, predicted with PSRORT. Also, the proteins present different transmembrane regions (TMpred), domains (INTERPRO, CDD and PRODOM) families (Pfam) and motifs (Motif finder and MOTIF search). The colors follow the same guidelines as in Fig. 1.

Table S8. Insertion sequences in the T6SS. The clusters of the T6SS of Xeu, Xcc3 and Xpm were examined for IS families and groups, defined in the database ISfinder ³¹

Download PDF

Journal Publication

published 06 Jan, 2021

Read the published version in BMC Microbiology →

Editorial decision: Minor revision
21 Nov, 2020
Review #1 received at journal
20 Nov, 2020
Review #2 received at journal
19 Nov, 2020
Reviewer #2 agreed at journal
18 Nov, 2020
Reviewers invited by journal
16 Nov, 2020
Reviewer #1 agreed at journal
16 Nov, 2020
Editor assigned by journal
15 Nov, 2020
Submission checks completed at journal
15 Nov, 2020
Editor invited by journal
15 Nov, 2020

You are reading this older preprint version

Read the latest preprint version →

The Type VI Secretion System of Xanthomonas phaseoli pv. manihotis is involved in virulence and in vitro motility

Status:

Journal Publication

Version 2

Abstract

Figures

Background

Methods

Results

Discussion

Conclusions

Abbreviations

Declarations

References

Supplementary Information

Supplementary Files

Status:

Journal Publication

Version 2