key: cord-1027801-3q8x2u5y
authors: Durán, David; Bernal, Patricia; Vazquez-Arias, David; Blanco-Romero, Esther; Garrido-Sanz, Daniel; Redondo-Nieto, Miguel; Rivilla, Rafael; Martín, Marta
title: Pseudomonas fluorescens F113 type VI secretion systems mediate bacterial killing and adaption to the rhizosphere microbiome
date: 2021-03-11
journal: Sci Rep
DOI: 10.1038/s41598-021-85218-1
sha: 16854a2b563715b0fbad40660d55d3db0ce6fa4f
doc_id: 1027801
cord_uid: 3q8x2u5y

The genome of Pseudomonas fluorescens F113, a model rhizobacterium and a plant growth-promoting agent, encodes three putative type VI secretion systems (T6SSs); F1-, F2- and F3-T6SS. Bioinformatic analysis of the F113 T6SSs has revealed that they belong to group 3, group 1.1, and group 4a, respectively, similar to those previously described in Pseudomonas aeruginosa. In addition, in silico analyses allowed us to identify genes encoding a total of five orphan VgrG proteins and eight putative effectors (Tfe), some with their cognate immunity protein (Tfi) pairs. Genes encoding Tfe and Tfi are found in the proximity of P. fluorescens F113 vgrG, hcp, eagR and tap genes. RNA-Seq analyses in liquid culture and rhizosphere have revealed that F1- and F3-T6SS are expressed under all conditions, indicating that they are active systems, while F2-T6SS did not show any relevant expression under the tested conditions. The analysis of structural mutants in the three T6SSs has shown that the active F1- and F3-T6SSs are involved in interbacterial killing while F2 is not active in these conditions and its role is still unknown.. A rhizosphere colonization analysis of the double mutant affected in the F1- and F3-T6SS clusters showed that the double mutant was severely impaired in persistence in the rhizosphere microbiome, revealing the importance of these two systems for rhizosphere adaption.

Distribution of T6SSs in the Pseudomonas fluorescens species complex. In silico analyses of the genomes of 134 strains belonging to the P. fluorescens species complex revealed that only 20% of them encode T6SS clusters (Table S3 ). The number of T6SS clusters in a single strain fluctuated from zero in P. fluorescens UK4 to three in P. fluorescens F113, and most strains contained two or three clusters (Table S3) . Overall, we identified sixty-two complete T6SS gene clusters distributed mostly in three main phylogenetic clades. We referred to these three groups as 1.1, 3, and 4A ( Fig. 1) following the previous nomenclature 7, 21 . The distribution of the clusters in these three groups is the same as in P. aeruginosa 2, 11 . No clusters were found corresponding to group 5, typical of Agrobacterium spp. and very few belonged to clusters 1.2 and 4B as in P. putida 38 . Each of these groups contains distinguishable genetic architecture and features (Fig. 2) , as described in the next section.

Pseudomonas fluorescens F113 contains three putative T6SSs. An extensive genomic analysis of the strain F113 using bioinformatics approaches (e.g., BLASTP, SMART, or Phyre-version 2) allowed us to identify a large number of T6SS-related ORFs (Tables S4-S7) . Most of the genes are clustered in three genomic regions that we have named F1-, F2-and F3-T6SS (Fig. 2 , Tables S4-S6). We identified a total of three hcp and eight vgrG genes,one of the hcp and six of the vgrG are orphan genes that were found scattered on the chromosome and not within the main T6SS clusters ( Fig. 2 and Table S7 ). Phylogenetic analysis revealed that similarly to P. aeruginosa PAO1, the F1-T6SS belongs to phylogenetic group 3, F2-T6SS to group 1.1 and the F3-T6SS to group 4A (Fig. 1) .

Each of the F113 T6SS clusters (F1-, F2-and F3-T6SS) consists of two divergently transcribed gene clusters that contain the 13 predicted core-genes (tssABCDEFGHIJKLM) ( Fig. 2 and Tables S4-S6) and genes putatively encoding accessory components as described below.

The large F1 cluster (44 Kb) contains a set of genes putatively encoding regulatory proteins previously described in P. aeruginosa: fha, tagF, pppA, ppkA, and tagRSTQ. Interestingly, it does not contain the gene that encodes the accessory component TagJ (Fig. 2 ) present in P. aeruginosa H1-T6SS 24, 65 and conserved in other clusters from the same phylogenetic group (Group 3). In those clusters, tagJ is found located downstream hcp and upstream tssE, and it has been recently described as an accessory component that stabilizes the sheath from the baseplate 24 . Instead, in this location, P. fluorescens, contains the orthologues of tae4-tai4 from Salmonella typhimurium, which encode an EI pair (see below). A gene encoding a DUF1795 domain-containing protein is found downstream vgrG1a within the F1-cluster. The product of this gene is also known as an EagR chaperone for a downstream Rhs type effector first identified in Serratia marcescens 40 . Two Rhs-type effectors are encoded downstream vgrG1a and the eagR chaperone in F113 (Fig. 2) and are described in the following section.

The F2-T6SS is a shorter cluster (around 29 Kb) with four genes putatively encoding regulatory proteins previously defined in P. aeruginosa 66 (stp, stk, fha2 and sfa2). All the core components are present within the cluster except for hcp2 (PSF113_1976) that is found in a different genomic region. A gene encoding a DUF4123 domain-containing protein chaperone 37, 48 (Tap2) is found downstream vgrG2a followed by a hypothetical protein and a lipase, that are likely to correspond to an immunity-toxin pair ( Fig. 2 and Table S5 ).

The F3-T6SS cluster, with an intermediate size of 35 Kb, does not contain genes encoding regulatory proteins ( Fig. 2 and Table S6 ), as opposed to P. aeruginosa, which contains sfa3, a gene encoding a putative sigma-54 transcriptional regulator enhancer. Interestingly, further down the T6SS cluster and three no-T6SS genes, we identified genes putatively encoding a Toxin Complex (TC) composed of one TcdB and two TccC subunits. The absence of the TcdA subunit, which forms an injection-like structure, has been observed in some strains 67 and it opens the possibility for this complex to be secreted through the T6SS. www.nature.com/scientificreports/ A total of eight vgrG genes are found in the chromosome of P. fluorescens F113. Three of these vgrG genes are within the F1-, F2-and F3-T6SS, whereas the other five genes are scattered over the chromosome (Table S7) . A phylogenetic study of the eight VgrG proteins shows that VgrG1a and 1b clustered in the same phylogenetic group as P. aeruginosa VgrG1 proteins. On the other hand, P. fluorescens VgrG2a/b, VgrG4, and VGrG5a/b are closely related among them and cluster with P. aeruginosa VgrG4, 5 and 6 ( Figure S1 ).

Identification of eight putative type 6 effectors in P. fluorescens F113. As stated before, genes encoding T6SS effector and their cognate immunity proteins are frequently found genetically associated to hcp and vgrG genes and in some cases to genes encoding T6SS chaperones or adaptors (i.e., EagR and Tap/Tec proteins) 15, 36, 37, 48, 68 . We have performed an in silico analysis that allowed us to identify a total of eight putative EI pairs in the proximity of F113 vgrG, hcp, eagR, and tap genes ( Fig. 2 and Tables S4-S7). These EI pairs have been named Tfe and Tfi standing for Type six F113 effector and immunity, respectively (Figs. 2, 3 and Tables S4-S7). www.nature.com/scientificreports/ Within the F1 cluster and downstream hcp1, we have identified the EI pair tfe1-tfi1. Tfe1 is orthologue to Salmonella typhimurium Tae4. Tae4 is a type VI amidase effector, which degrades the peptidoglycan of the cell wall and whose immunity protein is Tai4 43 . Furthermore, a gene encoding an EagR chaperone (eagR1) is found downstream vgrG1a in the F1-cluster and followed by two genes (tfe2 and tfe3) encoding Rhs effectors in a row. These effectors both have the same domain architecture: an N-terminal PAAR domain 69, 70 , an intermediate www.nature.com/scientificreports/ conserved Rhs domain 71 confined by specific RVxxxxxxxxG and PxxxxDPxGL motifs, and a C-terminal region encoding the toxic domain. The toxic domain of Tfe2 is a putative RES domain, whose orthologue in P. putida is known to cause depletion of intracellular NAD + levels leading to inhibition of cell growth 72 (Fig. 3 ). In parallel, the C-terminal domain of Tfe3 carries a putative nuclease of the HNH/ENDO VII superfamily with a conserved WHH domain. tfe2 and tfe3 putative cognate immunity pairs, named tfi2 and tfi3, are found immediately downstream of the genes encoding their effectors and have 92 and 136 amino acids, respectively (Table S4 and Figs. 2,  3) . Interestingly, tfi3 is a newly identified ORF in the P. fluorescens F113 genome and we have named the locus PSF113_5811.1 to indicate that the locus is located downstream PSF113_5811 (Table S4 and Fig. 2 ). Whereas Tfi2 immunity protein has no homologues or recognizable features, Tfi3 belongs to the Smi1/Knr4 family, which has been previously related to immunity proteins of polymorphic toxins 73 . Genes encoding putative effectors were also found within clusters F2 and F3, downstream vgrG2a and vgrG3 respectively and downstream orphan vgrG genes including vgrG1b, vgrG2b and vgrG5a (Fig. 2) . The tfe4, tfe7, and tfe8 genes within the F2 and vgrG2b and vgrG5a operons, respectively, are genetically associated with genes encoding Tap chaperones (tap2, tap2b, and tap5a) (Fig. 2 ). Tfe4 is a lipase from class 3, putatively targeting eukaryotic and prokaryotic membranes by hydrolysing long-chain acyl-triglycerides from lipids into di-and monoglycerides, glycerol, and free fatty acids. Tfi4 is the cognate immunity protein of Tfe4 and it does not contain any known function or recognisable domain (Table S5 ). Tfe7 is a homologue of a ribosomal RNA large subunit methyltransferase D required for the full methylation of 23S ribosomal RNA. Interestingly, no gene encoding an immunity pair can be found linked to this effector gene, indicating that it might not be necessary for the producer strain. This effector could replace a functional methyltransferase and only be toxic in those strains where this enzyme is necessary for bacterial viability. A similar case has been described in P. aeruginosa PAO1, in which the effector Tse8 is not linked to an immunity pair and replaces a functional component of the transamidosome complex only present in some strains 74 . Tfe8 may represent a new type of effector identified for the first time in this work and contains a MatE domain. None of these three Tap-linked effectors contains a conserved N-terminal MIX motif considered a marker for T6SS effectors 75 and frequently found in Tap-associated effectors 38 .

Gene tfe5 is located downstream vgrG3 and linked to a small PAAR-encoding gene. Tfe5 is a homologue of the TseF toxin from P. aeruginosa, an iron scavenging effector that interacts with PQS vesicles to help bacteria grow under very limiting conditions 76 . A small gene located between clpv3 and vgrG3 encodes a protein of 188 amino acids with no recognizable features but has a C-terminal domain of 50 amino acids that resembles a coronavirus RNA-binding domain according to a Phyre prediction (confidence: 80.5). The tfe6 is linked to the orphan VgrG1b cluster and has a similar genetic architecture to P. aeruginosa PAO1 VgrG1b cluster 53 . In PAO1 cluster, PA0095, PA0096, PA0097, PA0098, PA0099, PA0100, and PA0101 encode VgrG1b, an OB-fold, an immunoglobulin-like and a thiolase-like protein, a PAAR protein with a C-terminal cytotoxic domain, an immunity protein, and a heat repeat-containing protein respectively (Table S7) . Pseudomonas fluorescens F113 vgrG1b cluster (PSF113_2885-PSF113_2890) only differs from the PAO1 one in the sequence of the toxic domain and the immunity pair, a common characteristic of this genetic island previously described by 53 . The specific function of Tfe6 and the mechanistic of the cognate immunity pair Tfi7 remains unknown.

No putative effector-encoding genes were identified linked to hcp2, vgrG4, and vgr5b clusters and we hypothesised that they could be found somewhere else in the chromosome as orphan T6SS effectors as described before in other T6SS clusters, including Vibrio proteolyticus 77 among others.

In summary, we identified eight putative T6SS effectors in the F113 genome. Three of them, Tfe2, Tfe3, and Tfe6, contain an N-terminal PAAR-domain (Figs. 2, 3) and are considered "specialised" effectors, whether the others that are not fused to any T6SS component are considered "cargo" effectors. F1-and F3-T6SSs genes are expressed in P. fluorescens F113. In order to determine the expression of genes encoding the components of the three T6SSs in P. fluorescens F113, we used a new dataset from an RNA-Seq study of F113 and derivatives grown under different culture conditions and in the alfalfa (Medicago sativa) rhizosphere. For this analysis, we chose the structural genes from each of the three systems (tssABCDEFGHI-JKLM). As shown in Fig. 4A , F1-and F3-T6SS cluster genes were significantly expressed under all tested conditions: growing in liquid cultures of Minimal Sucrose-Asparagine (SA) at exponential and stationary phases and colonising the alfalfa rhizosphere. On the contrary, F2-T6SS genes showed little, if any, activity under the same conditions. Expression of F1-T6SS genes was similar in all culture conditions regardless of the growth phase analysed, and these genes were also expressed during rhizosphere colonization, suggesting that the F1-T6SS is constitutively expressed in F113. The expression of F3-T6SS genes is lower than that of the F1-T6SS genes in all conditions. This situation is similar to P. putida, which contains three clusters named K1-, K2-and K3-T6SS, that do not belong to the same phylogenetic groups of P. fluorescens F113. In P. putida, the K1-T6SS is constitutively expressed under laboratory conditions and is active in vitro and in planta assays 38 ; however, no activity has been reported for K2-and K3-T6SS clusters to date. Conversely, in the case of P. aeruginosa that harbours the same groups as F113 and share regulatory components within the T6SS clusters, T6SS genes are not constitutively expressed and all of them are tightly regulated at transcriptional, post-transcriptional, and post-translational levels 11, [78] [79] [80] . In PAO1, the three systems have the common regulator Rsm 52 and the widely studied H1-T6SS is known to also be regulated by many other factors including RetS, GacA/S, TagF, PpkA/PppA, TagQRST, and FHA 2, 33, 34, [80] [81] [82] .

To study the regulation of T6SSs during the process of rhizosphere colonization, we also analysed the T6SSs genes with differential expression in RNA-Seq data from the P. fluorescens F113 amrZ and fleQ mutants compared to the wild-type strain grown in the rhizosphere of alfalfa. AmrZ and FleQ are two transcription factors (TFs) crucial during competitive rhizosphere colonization in this bacterium [83] [84] [85] [86] [87] www.nature.com/scientificreports/ and it has been previously linked to the control of T6SS in P. aeruginosa 52 . FleQ belongs to the NtrC/NifA family and has been related to the regulation of T6SS and c-di-GMP metabolism in P. putida 89 . The differential gene expression analyses shown in Fig. 4B has revealed that in the rhizospheric environment, AmrZ functions as a negative transcriptional regulator of F1-and F3-T6SS, while FleQ acts as a positive regulator. This antagonist role of both transcriptional regulators fits with the proposed model for the AmrZ/FleQ hub, which has been proposed in F113 to act as an oscillator with opposing effects in gene expression, in order to integrate the bacterial responses to the environment 83 . In addition to the vrgG genes related to the F1-and F3-T6SS, most orphan vrgG genes are also negatively regulated by AmrZ under the tested conditions. However, FleQ can act both as a positive and negative regulator of certain vgrG genes. These results show that both AmrZ and FleQ regulate T6SSs in F113, as previously shown for AmrZ in P. aeruginosa 52 and FleQ in P. putida 89 .

in the antibacterial activity of some strains due to the injection of T6SS toxins into competitor target cells 5, 7, 30, 42 . Therefore, to analyse the role of T6SS in P. fluorescens F113 in inter-bacterial competition, we performed bacterial killing assays as previously described 38 . In these assays, P. fluorescens F113 and its isogenic mutants were used as predators, and Escherichia coli harbouring a pK18mobsacB plasmid, containing the lacZ gene that confers blue colour to the colony in the presence of X-gal, was used as prey. Bacteria were mixed in a 1:1 ratio (predator:prey), co-cultured for 5 h and a sample of each mix was grown on selective plates for 48 h. Additionally, serial dilutions of the different assays were plated on selective media for predator and prey quantification.

The tssA genes of P. fluorescens F113 were selected as targets for the construction of insertional mutants. Mutants affected in each of the three systems were used as predators: tssA1 − , tssA2 − and tssA3 − for F1-, F2-and F3-T6SS, respectively. Additionally, a double mutant tssA1 − /tssA3 − was constructed and used as predator. In the bacterial competition assays (Fig. 5) , the wild-type strain was able to kill E. coli cells efficiently as observed for the significant prey survival reduction. However, single and double mutants in the structural genes tssA1 and tssA3 (F1-and F3-T6SSs) were affected in E. coli killing (Fig. 5) , showing a survival rate of the prey similar to the control without predator. These results indicate that both systems are functional and have bactericidal activity, as www.nature.com/scientificreports/ has been shown before for other T6SSs in pseudomonads 4, 90, 91 . By contrast, the tssA2 mutant (F2-T6SS) exhibited the same capacity to outcompete E. coli that the wild type strain, indicating that the F2-T6SS is not involved in the antibacterial activity of P. fluorescens F113 against E. coli under our experimental conditions. It is interesting to note that we have not observed expression of the genes encoding this system in any of the tested conditions, suggesting that the lack of function might be a consequence of lack of expression and therefore, the F2-T6SS could be functional under other yet-unknown conditions. Lack of expression and/or antibacterial activity has been observed in other pseudomonads in laboratory conditions, like P. aeruginosa wild type strain 2,52 and P. putida K2-and K3-T6SS 38 .

Pseudomonas fluorescens F113 T6SSs are important for the adaption to the rhizosphere microbiome. It has been described that the T6SSs have a role in modulating and shaping the natural microbiota and, in the case of plant-associated bacteria, these weapons are of interest for bacterial persistence in the plant niche, reviewed in 92 . Also Vacheron et al., 2019, showed that T6SS of P. protegens contributed to the invasion of the gut microbiome of an insect. Therefore, we wanted to know whether F1 and F3 T6SS could play a role in F113 competence in the rhizosphere. We inoculated 7-day-old tomato plants growing in agricultural, non-sterile soil microcosms, with the wild type strain and the double mutant F1 − /F3 − strain. Bacteria from microcosms were isolated 2 weeks after inoculation and F113 derivatives were selected by using their rifampicin marker resistance. As shown in Fig. 6 , when the inoculation is done with the double mutant strain, there is a significant (90%) reduction in the bacterial recovery after rhizosphere colonization compared with the wild-type strain. This result shows that the two tested T6SSs play a relevant role in invading, establishing and/or persisting in the tomato rhizosphere microbiome. The effect of these T6SSs in microbiome adaption is likely due to their role in the interbacterial competition that might confer F113 with the capacity to outcompete foes.

Three T6SS clusters (F1-, F2-and F3-T6SS) and eight effectors (Tfe1, 2, 3, 4, 5, 6, 7, 8) together with 6 immunity proteins (Tfi1, 2, 3, 4, 6 y 8) have been identified in P. fluorescens F113. At least two of these systems, F1-and F3-T6SS are functional in laboratory conditions and more importantly in the rhizosphere and possess bactericidal activity. The systems are crucial elements in the colonization of the rhizosphere, most likely providing by providing F113 with the capacity to fight competitors from the rhizosphere microbiome. coli growth for each competition experiment. Blue colour in bacterial patches on LB plates supplemented with X-gal and kanamycin indicates E. coli growth. Error bars indicate the mean ± s.d. of three biological replicates, and significance was calculated using ANOVA test (*P < 0.01); ns, not significant differences when compared to non-competing E. coli. Bioinformatic analyses. Pseudomonas gene sequences were obtained from the Pseudomonas Genome database 95 . BLASTP analyses were performed at the NCBI website 96 and amino acid sequence searches using SMART 97, 98 . The Protein Homology/analogy Recognition Engine (Phyre 2 ) server was used to perform structural-base homology prediction 99 . The phylogenetic tree was constructed using MEGA X 100 . PSORTb v 3.0.2 software was used to predict subcellular location of proteins 101 , TMHMM server v.2.0 to predict transmembrane domains 102 , and SignalP to predict signal peptides 103 . The molecular biology tools included in Benchling platform were used to identify novel open reading frames (ORFs) (Benchling, Inc., https ://bench ling.com/acade mic) and the putative coding proteins were run using NCBI BLASTP tool to determine the degree of conservation.

Construction of T6SS mutants. Insertional single mutants in PSF113_2422 (tssA3), 5797 (tssA1) and 5829 (tssA2) genes were constructed by electroporation of the plasmid vector pCR2.1-TOPO (Invitrogen) carrying an internal region of the corresponding gene (400 bp ca.), in the case of the double mutant tssA1 − /A3 the plasmid vector pG18mob2 was employed. Single and double recombinant mutants were selected by kanamycin (25 μg ml −1 ) and gentamycin (4 μg ml −1 ) resistance in SA medium and checked by PCR and Southern Blot. Primers used are listed in Supplementary Table S2 .

Bacterial RNA isolation from cultures. RNA was isolated from P. fluorescens F113 cultures grown in Minimal Sucrose-Asparagine (SA) medium at exponential (OD 600 = 0.6) and stationary phase (OD 600 = 1.2) as indicated in 84 .

Bacterial RNA isolation from the rhizosphere. Rhizosphere colonization of alfalfa (Medicago sativa var. Resis) plants was performed essentially as described in 104 . Seven-day-old seedlings growing in Falcon tubes, using vermiculite as substrate, were inoculated with 1 mL (1 × 10 8 CFU) of P. fluorescens F113 or derivatives amrZand fleQmutants. RNA was extracted seven days post-inoculation. Aerial parts of alfalfa plants were removed, 4 mL of Phosphate-buffered saline (PBS) and 6 mL RNAlater 105 added and tubes were vortexed for 2 min to resuspend bacterial cells. The mix of the 32 preparations per sample was filtered through four layers of sterile muslin cloth in pyrex funnels and separated into six 50 mL tubes. The filtrate was centrifuge 1 min at 1000 rpm and 4 °C. The supernatant was transferred to a sterile tube and centrifuged at 10,000 rpm for 20 min and 4 °C. Supernatants were discarded and pelleted cells dried before liquid nitrogen freezing. RNA isolation was performed as indicated in 84 .

RNA sequencing. Qubit fluorometer quality assessment, rRNA depletion, strand-specific library construction, and sequencing were performed by Novogene Co., Ltd. (Beijing, China, and Cambridge, UK) using Illumina HiSeq paired-end, 2 × 150 bp.

Bioinformatic RNA-Seq data processing. After sequencing the RNA, the quality of the raw reads was checked using FastQC. Then, sequence reads were clipped and filtered using Trimmomatic v 0.35 106 to remove www.nature.com/scientificreports/ chaperones and low-quality nucleotides, defining a four nts sliding window with an average phred quality of 15 and 50 nts as minimum read length. High-quality reads were directly used for transcript-level quantification using Salmon software 107 , which performs a quasi-mapping to the new annotation of the P. fluorescens F113 CDSs (GenBank: NC_016830) and transcript quantification. Normalization of counts and differential gene expression was calculated with DESeq2 1.24.0 R package 108 . Differential gene expression comparisons were made setting a threshold for log 2 fold change (mutant/wild-type) of ≤ − 1/≥ 1 and a stringent p-adjusted value cutoff ≤ 0.001. RNA-Seq reads have been deposited in to the NCBI Sequence Read Archive database and it is available under the BioProject PRJNA419480: BioSamples accessions: F113 culture in exponential growth phase (SAMN17839758 and SAMN17839763), F113 culture in stationary growth phase (SAMN17839759 and SAMN17839764), F113 from the rhizosphere (SAMN17839757 and SAMN17839762), F113 amrZ mutant from the rhizosphere (SAMN17839760 and SAMN17839765) and F113 fleQ mutant from the rhizosphere (SAMN17839761 and SAMN17839766).

Interbacterial competition assays. Competition assays to assess T6SSs were performed on solid media plates according to the previously reported protocol 51 . Pseudomonas fluorescens F113 or isogenic derivatives (predator) and E. coli DH5α containing pK18mobsacB (prey) were grown overnight in LB medium. Next morning, each culture was adjusted to OD 600 of 1.0. 100 μL of predator and 100 μL of prey strains were mixed and co-cultured for 5 h at 200 rpm and 28 °C. 20 μL of each culture was spotted onto LB-agar supplemented with 5-bromo-4-chloro-3-indolyl-D-galactopyranoside (X-gal) and kanamycin and incubated at 28 °C. Additionally, serial dilutions of the different assays were plated to quantify the number of CFUs in each of them. Three biologically independent experiments were performed.

Rhizosphere colonization analysis. Tomato seeds (Rebelion F1, Vilmorin, France) were sterilized with 70 % ethanol and 5 % hypochlorite, washed and germinated on sterile 1.0 % agar plates for 24 h at 28 °C. Seedlings (one per tube) were then transferred into sterile 50 mL Falcon tubes containing 16 g of a mix 1:1 of agricultural soil and sterile sand (Merck KGaA, Darmstadt, Germany). Each tube was embedded with 3 mL of autoclaved water, and incubated in a controlled environment room (25 °C, 16 h light cycle). After seven days, each tube/plant was inoculated with 1 mL (1 × 10 6 CFU/mL) culture of either the P. fluorescens F113 wild type strain or the double mutant (tssA1 − /A3 − ) strain. Plants were grown for another two weeks, after which aerial parts were removed; 10 mL of saline solution (0.85 %) was added to each tube and vortexed thoroughly to resuspend the bacteria. After decantation of the soil and sand matrix, serial dilutions of the supernatant were plated onto SA plates supplemented with rifampicin (100 μg mL −1 ) and nystatin (50 μg mL −1 ). Assays were performed with six plants per condition tested and non-inoculated plants were used as a negative control.

Statistical analyses. The normal distribution of data was checked with Shapiro-Wilk test. When distribution was normal, the data were analysed with one-way ANOVA test, where multiple pairwise-comparisons between strains were performed with Tukey HSD test. When normality was not met, data were analysed with the Kruskal-Wallis non-parametric test. All data were analysed with R package version 3.6.3 109 . 

Identification of a conserved bacterial protein secretion system in Vibrio cholerae using the Dictyostelium host model system

A virulence locus of Pseudomonas aeruginosa encodes a protein secretion apparatus

Type VI secretion system translocates a phage tail spikelike protein into target cells where it cross-links actin

Type VI secretion delivers bacteriolytic effectors to target cells

Type VI secretion systems in plant-associated bacteria

Type VI secretion: a beginner's guide

Dissecting the bacterial type VI secretion system by a genome wide in silico analysis: what can be learned from available microbial genomic resources?

Evaluation of a Yersinia pestis mutant impaired in a thermoregulated type VI-like secretion system in flea, macrophage and murine models

The cluster 1 type VI secretion system is a major virulence determinant in Burkholderia pseudomallei

The Type VI secretion system of Rhizobium etli Mim1 has a positive effect in symbiosis

Quorum sensing differentially regulates Pseudomonas aeruginosa type VI secretion locus I and homologous loci II and III, which are required for pathogenesis

Quorum sensing coordinates brute force and stealth modes of infection in the plant pathogen Pectobacterium atrosepticum

Secretome analysis uncovers an Hcp-family protein secreted via a type VI secretion system in Agrobacterium tumefaciens

A type VI secretion system of Pseudomonas aeruginosa targets a toxin to bacteria

Identification of T6SS-dependent effector and immunity proteins by Tn-seq in Vibrio cholerae

A type VI secretion system is involved in Pseudomonas fluorescens bacterial competition

The type VI secretion system deploys antifungal effectors against microbial competitors

In vivo expression technology identifies a type VI secretion system locus in Burkholderia pseudomallei that is induced upon invasion of macrophages

The bacterial type VI secretion machine: yet another player for protein transport across membranes

The type VI secretion toolkit

Genomic analysis of the type VI secretion systems in Pseudomonas spp.: novel clusters and putative effectors uncovered

Priming and polymerization of a bacterial contractile tail structure

TssA forms a gp6-like ring attached to the type VI secretion sheath

A novel stabilization mechanism for the type VI secretion system sheath

Cryo-EM reconstruction of Type VI secretion system baseplate and sheath distal end

Diverse roles of TssA-like proteins in the assembly of bacterial type VI secretion systems

The type VI secretion system sheath assembles at the end distal from the membrane anchor

A view to a kill: the bacterial type VI secretion system

VgrG and PAAR proteins define distinct versions of a functional type VI secretion system

Aim, load, fire: the type VI secretion system, a bacterial nanoweapon

The Type VI secretion system: a versatile bacterial weapon

Structure and activity of the type VI secretion system

TagR promotes PpkA-catalysed type VI secretion activation in Pseudomonas aeruginosa

An ABC transporter and an outer membrane lipoprotein participate in posttranslational activation of type VI secretion in Pseudomonas aeruginosa

In vivo TssA proximity labelling during type VI secretion biogenesis reveals TagA as a protein that stops and holds the sheath

Agrobacterium tumefaciens deploys a superfamily of type VI secretion DNase effectors as weapons for interbacterial competition in planta

Chimeric adaptor proteins translocate diverse type VI secretion system effectors in Vibrio cholerae

The Pseudomonas putida T6SS is a plant warden against phytopathogens

The opportunistic pathogen Serratia marcescens utilizes type VI secretion to target bacterial competitors

Intraspecies competition in Serratia marcescens is mediated by type VI-secreted Rhs effectors and a conserved effector-associated accessory protein

Structure and regulation of the type VI secretion system

Type VI secretion system effectors: poisons with a purpose

A widespread bacterial type VI secretion effector superfamily identified using a heuristic approach

Effector(-)immunity pairs provide the T6SS nanomachine its offensive and defensive capabilities

Comparative genomics of the Type VI secretion systems of Pantoea and Erwinia species reveals the presence of putative effector islands that may be translocated by the VgrG and Hcp proteins

Identification of protein secretion systems in bacterial genomes

Genetically distinct pathways guide effector export through the type VI secretion system

Identification of divergent type VI secretion effectors using a conserved chaperone domain

The T6SSs of Pseudomonas aeruginosa strain PAO1 and their effectors: beyond bacterial-cell targeting

Haemolysin coregulated protein is an exported receptor and chaperone of type VI secretion substrates

A visual assay to monitor T6SS-mediated bacterial competition

RsmA and AmrZ orchestrate the assembly of all three type VI secretion systems in Pseudomonas aeruginosa

The Pseudomonas aeruginosa T6SS-VgrG1b spike is topped by a PAAR protein eliciting DNA damage to bacterial competitors

An interbacterial NAD(P)(+) glycohydrolase toxin requires elongation factor Tu for delivery to target cells

Involvement of type VI secretion system in secretion of iron chelator pyoverdine in Pseudomonas taiwanensis

The extracellular matrix protects Bacillus subtilis colonies from Pseudomonas invasion and modulates plant co-colonization

Contribution of the Pseudomonas fluorescens MFE01 Type VI secretion system to biofilm formation

Genomic analysis of the biocontrol strain Pseudomonas fluorescens Pf29Arp with evidence of T3SS and T6SS gene expression on plant roots

T6SS contributes to gut microbiome invasion and killing of an herbivorous pest insect by plant-beneficial Pseudomonas protegens

Diverse NADase effector families mediate interbacterial antagonism via the type VI secretion system

Identification, structure, and function of a novel type VI secretion peptidoglycan glycoside hydrolase effectorimmunity pair

Isolation of 2,4-diacetylphloroglucinol from a fluorescent pseudomonad and investigation of physiological parameters influencing its production

Genome sequence of the biocontrol strain Pseudomonas fluorescens F113

Genome sequence reveals that Pseudomonas fluorescens F113 possesses a large and diverse array of systems for rhizosphere function and host interaction

Coevolution of the ATPase ClpV, the sheath proteins TssB and TssC, and the accessory protein TagJ/HsiE1 distinguishes type VI secretion classes

An rhs gene linked to the second type VI secretion cluster is a feature of the Pseudomonas aeruginosa strain PA14

The role of TcdB and TccC subunits in secretion of the Photorhabdus Tcd toxin complex

The VgrG proteins are "a la carte" delivery systems for bacterial type VI effectors

PAAR-repeat proteins sharpen and diversify the type VI secretion system spike

Structural biology of type VI secretion systems

The BC component of ABC toxins is an RHS-repeatcontaining protein encapsulation device

The RES domain toxins of RES-Xre toxin-antitoxin modules induce cell stasis by degrading NAD+

A novel immunity system for bacterial nucleic acid degrading toxins and its recruitment in various eukaryotic and DNA viral systems

Discovery of a Pseudomonas aeruginosa Type VI secretion system toxin targeting bacterial protein synthesis using a global genomics approach

Marker for type VI secretion system effectors

A Pseudomonas T6SS effector recruits PQS-containing outer membrane vesicles for iron acquisition

Type VI secretion system MIX-effectors carry both antibacterial and anti-eukaryotic activities

The GacS/GacA signal transduction system of Pseudomonas aeruginosa acts exclusively through its control over the transcription of the RsmY and RsmZ regulatory small RNAs

The Pseudomonas aeruginosa sensor RetS switches type III and type VI secretion via c-di-GMP signalling

Threonine phosphorylation post-translationally regulates protein secretion in Pseudomonas aeruginosa

Separate inputs modulate phosphorylation-dependent and -independent type VI secretion activation

TagF-mediated repression of bacterial type VI secretion systems involves a direct interaction with the cytoplasmic protein Fha

Genome-wide analysis of the FleQ direct regulon in Pseudomonas fluorescens F113 and Pseudomonas putida KT2440. Sci. Rep. 8, 13145

AmrZ is a major determinant of c-di-GMP levels in Pseudomonas fluorescens F113

AmrZ is a global transcriptional regulator implicated in iron uptake and environmental adaption in P. fluorescens F113

Analysis of Pseudomonas fluorescens F113 genes implicated in flagellar filament synthesis and their role in competitive root colonization

Efficient rhizosphere colonization by Pseudomonas fluorescens f113 mutants unable to form biofilms on abiotic surfaces

Pseudomonas aeruginosa AlgZ, a ribbon-helix-helix DNA-binding protein, is essential for alginate synthesis and algD transcriptional activation

FleQ regulates both the type VI secretion system and flagella in Pseudomonas putida

The Pseudomonas aeruginosa T6SS delivers a Periplasmic toxin that disrupts bacterial cell morphology

Tit-for-tat: type VI secretion system counterattack during bacterial cell-cell interactions

Causalities of war: the connection between type VI secretion system and microbiota

Molecular Cloning: A Laboratory Manual

Effect of Pseudomonas putida and a synthetic iron chelator on induction of soil suppressiveness to Fusarium wilt pathogens

Enhanced annotations and features for comparing thousands of Pseudomonas genomes in the Pseudomonas genome database

BLAST: a more efficient report with usability improvements

20 years of the SMART protein domain annotation resource

SMART: recent updates, new developments and status in 2015

The Phyre2 web portal for protein modeling, prediction and analysis

molecular evolutionary genetics analysis across computing platforms

PSORTb 3.0: improved protein subcellular localization prediction with refined localization subcategories and predictive capabilities for all prokaryotes

Predicting transmembrane protein topology with a hidden Markov model: application to complete genomes

SignalP 4.0: discriminating signal peptides from transmembrane regions

Adaptive remodeling of the bacterial proteome by specific ribosomal modification regulates Pseudomonas infection and niche colonisation

Biological stability of RNA isolated from RNAlater-treated brain tumor and neuroblastoma xenografts

Trimmomatic: a flexible trimmer for Illumina sequence data

Salmon provides fast and bias-aware quantification of transcript expression

Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2

A Language and Environment for Statistical Computing

This work has been funded by Ministerio de Ciencia, Innovación y Universidades FEDER/EU Grant RTI2018-093991-B-I00. EB-R was granted by FPU-MECD program (FPU16/05513). D.G.-S. was funded by FPU-MECD fellowship program (FPU14/03965). PB was funded by an InterTalentum Fellowship at UAM (reference 010101_ T6SS-SUPERBUG). "Research Challenges" 2018 R+D+i Project RTI2018-096936-J-I00 (funded by the Spanish Ministry of Science, Innovation and Universities) and FEDER18-1264962 grant funded by the Regional Government of Andalusia.

The authors declare no competing interests.

The online version contains supplementary material available at https ://doi. org/10.1038/s4159 8-021-85218 -1.Correspondence and requests for materials should be addressed to M.M.Reprints and permissions information is available at www.nature.com/reprints.Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.