Transcriptional organization and regulation of the Pseudomonas putida K1 type VI secretion system gene cluster

1 Transcriptional organization and regulation of the Pseudomonas putida K1 type VI secretion system genecluster PatriciaBernal1,2,3,*, CristinaCivantos1, DanielPacheco- Sánchez1, José M.Quesada1, AlainFilloux2,4,* and María A.Llamas1,* RESEARCH ARTICLE Bernal etal., Microbiology 2023;169:001295 DOI 10.1099/mic.0.001295 Received 11 October 2022; Accepted 20 December 2022; Published 27 January 2023 Author affiliations: 1Department of Environmental Protection, Estación Experimental del Zaidín (CSIC), Granada, Spain; 2MRC Centre for Molecular Bacteriology and Infection, Department of Life Sciences, Imperial College London, London, UK; 3Departamento de Microbiología, Facultad de Biología, Universidad de Sevilla, 41012 Seville, Spain; 4Singapore Centre for Environmental Life Sciences Engineering. Nanyang Technological University, Singapore. *Correspondence: Patricia Bernal, pbernal@ us. es; Alain Filloux, a. filloux@ imperial. ac. uk; María A. Llamas, marian. llamas@ eez. csic. es Keywords: type VI secretion system; gene regulation; RetS; GacS–GacA; RpoN; RpoS; Pseudomonas; FleQ. Abbreviations: bEBP, bacterial enhancerbinding protein; bp, base pair; CSS, cell surface signalling; ECF, extracytoplasmic function; H- NS, histonelike proteins; QS, quorum sensing; RNAP, RNA polymerase; TCS, twocomponent system; TF, transcriptional factor; T6SS, type VI secretion system; TSS, transcription start site. Two supplementary figures and four supplementary tables are available with the online version of this article. 001295 © 2023 The Authors This is an openaccess article distributed under the terms of the Creative Commons Attribution License. This article was made open access via a Publish and Read agreement between the Microbiology Society and the corresponding author’s institution. Abstract The type VI secretion system (T6SS) is an antimicrobial molecular weapon that is widespread in Proteobacteria and offers competitive advantages to T6SS- positive microorganisms. Three T6SSs have recently been described in Pseudomonas putida KT2440 and it has been shown that one, K1- T6SS, is used to outcompete a wide range of phytopathogens, protecting plants from pathogen infections. Given the relevance of this system as a powerful and innovative mechanism of biological control, it is critical to understand the processes that govern its expression. Here, we experimentally defined two transcriptional units in the K1- T6SS cluster. One encodes the structural components of the system and is transcribed from two adjacent promoters. The other encodes two hypothetical proteins, the tip of the system and the associated adapters, and effectors and cognate immunity proteins, and it is also transcribed from two adjacent promoters. The four identified promoters contain the typical features of σ70- dependent promoters. We have studied the expression of the system under different conditions and in a number of mutants lacking global regulators. P. putida K1- T6SS expression is induced in the stationary phase, but its transcription does not depend on the stationary σ factor RpoS. In fact, the expression of the system is indirectly repressed by RpoS. Furthermore, it is also repressed by RpoN and the transcriptional regulator FleQ, an enhancerbinding protein typically acting in conjunction with RpoN. Importantly, expression of the K1- T6SS gene cluster is positively regulated by the GacS–GacA twocomponent regulatory system (TCS) and repressed by the RetS sensor kinase, which inhibits this TCS. Our findings identified a complex regulatory network that governs T6SS expression in general and P. putida K1- T6SS in particular, with implications for controlling and manipulating a bacterial agent that is highly relevant in biological control. INTRODUCTION The type VI secretion system (T6SS) is a bacterial contractile nanomachine found in approximately 25 % of Gramnegative bacteria. The T6SS is used to deliver effectors/toxins into neighbouring cells in a contactdependent manner. Since its discovery in 2006, most T6SS studies have been performed in bacterial human pathogens such as Pseudomonas aeruginosa, Vibrio cholerae and enteropathogenic Escherichia coli strains, and it was initially assumed that T6SS’s main role was in delivering effectors inside eukaryotic cells [1–3]. However, more recent studies revealed that most T6SS effectors are antibacterial toxins and are produced and secreted by both pathogenic and nonpathogenic bacteria [4]. It is now generally recognized that the T6SS is important for interbacterial competition and that in nonpathogenic environmental bacteria this system is mainly used to outcompete foes [5–7]. The T6SS is made up of 12 core components that form a membrane complex (TssJLM) that anchors the system to the cell envelope and holds a baseplate platform (TssEFGK). From this platform and primed by a ring of TssA proteins, a tail polymerizes from the tip (VgrG and PAAR) and is followed by an inner tube of Hcp rings wrapped within a contractile sheath of TssBC proteins [8] (Fig.1a). At the structural level, these proteins and the ATPase that recycles the sheath (ClpV) are highly conserved among OPEN ACCESS 2 Bernal etal., Microbiology 2023;169:001295 the T6SS clusters, with few exceptions. Despite the general conservation of the structural components of the apparatus, different phylogenetic groups have been described [9, 10] and the T6SS clusters have been shown to differ in the organization of the genes encoding the core components (Fig.1b) [10], the presence of genes encoding a variable number of accessory proteins (TagP, TagJ, TagL, TagA, TagB) [11–13] (Fig.1b, dark blue), genes encoding effector–immunity pairs and, more importantly, genes encoding specific regulatory elements that control the T6SS at different levels (Fig.1b, purple). The activity of the T6SSs is tightly regulated at each possible stage, from transcriptional to translational and posttranslational levels. It is common to find genes that encode T6SS regulatory elements within the T6SS clusters, with the most common elements being PpkA/PppA serine–threonine kinase/phosphatase systems, bacterial enhancerbinding proteins (bEBPs) and other regulators, such as TagF and forkheadassociated (FHA) domaincontaining proteins [14–18]. For example, the P. aeruginosa H1- T6SS cluster encodes seven different regulatory components, including tagF1, pppA, ppkA and tagQRST [19–21] (Fig.1b, purple). The PpkA/PppA system regulates the activity of the P. aeruginosa H1- T6SS at the posttranslational level by phosphorylating/ dephosphorylating an FHA domaincontaining protein termed FHA1 that is required for the assembly of the H1- T6SS and the secretion of Hcp1 [14]. The accessory proteins TagOPQRST form a transmembrane signalling pathway that acts upstream of the PpkA/PppA, ultimately activating the H1- T6SS under surface growth conditions [19, 20]. This additional level of regulation is not widespread, but it can be found in other species of Pseudomonas such as P. fluorescens [6]. Proteins containing the FHA domain are found encoded in many other T6SS clusters, including those from Serratia marcescens [22] and Agrobacterium tumefaciens [23]. In A. tumefaciens, FHA interacts with a phosphorylated (p-) membrane complex component p- TssL, a process crucial for the interaction of TssL with Hcp and activation of type VI subassembly and secretion [23]. The phosphorylationindependent pathway of P. aeruginosa H1- T6SS is governed by TagF [24], a posttranslational repressor that can be found in many other T6SS clusters (Fig.1b). In P. aeruginosa, A. tumefaciens and S. marcescens, TagF interacts with p-FHA, impacting on T6SS assembly [16, 25]. T6SSs are also regulated by global regulatory mechanisms, including quorum sensing (QS), twocomponent systems (TCSs), histonelike proteins (H- NS), transcriptional factors (TFs), σ factors and small regulatory RNAs [26]. These global regulatory mechanisms control the expression of many processes in the cell and are widespread in bacteria. Importantly, although T6SS clusters have been horizontally transferred, the production of these systems has been incorporated into former regulatory networks. This observation suggests that these master regulators will allow coordination of T6SS expression, assembly and secretion, necessary to outcompete foes, with other important and equally energetically expensive processes to minimize energy costs. QS, which controls gene expression in response to bacterial cell population, is a major regulatory pathway for T6SS gene expression [26]. In P. aeruginosa, production of the three T6SSs of this bacterium is controlled by the LasRI and the MvfR (also known as PqsR) QS systems [27], whereas the CviR QS regulator controls the expression of the Chromobacterium violaceum T6SS cluster [28]. Furthermore, the most common regulatory mechanism governing the T6SS expression is TCSs. TCSs couple the sensing of an environmental signal to an adaptive response by a sensor kinase and a response regulator, allowing bacteria to adapt to different environmental conditions. Many T6SS loci are regulated by the GacS–GacA TCS, which Fig. 1. The T6SS apparatus. (a) Schematic representation of a T6SS after being fired. The VgrG/PAAR tip complex (beige/cyan), the Hcp tube (cyan) and the coupled effectors (grey) are ejected from the attacker into a target cell (pink). The membrane complex (TssJLM) and baseplate (TssEFGK) structures are shown in gold, and the sheath (TssBC) and the cap protein TssA are coloured in dark red and beige, respectively. (b) Schematic representation of the genomic organization of the structural genes of the H1- T6SS cluster from P. aeruginosa, V. cholerae and the K1-, K2- and K3- T6SS clusters from P. putida. Genes encoding core components are shown in light blue, accessory components are depicted in dark blue and regulatory elements are coloured in purple. 3 Bernal etal., Microbiology 2023;169:001295 is additionally modulated by other sensor kinases, such as RetS and/or LadS [29]. This signalling cascade converges on small RNAs (e.g. rsmZ and rsmY in P. aeruginosa), which control the function of RNA- binding proteins (e.g. RsmA) that often act as translational repressors of T6SSs [30]. In several bacteria, histonelike proteins repress the transcription of T6SS genes by their general mechanism, i.e. binding to AT- rich regions. For instance, MvaT, a P. aeruginosa H- NS- like protein, represses the H2- and H3- T6SS loci [31] and the TurA H- NS protein seems to silence the P. putida K2- T6SS [32]. Transcriptional factors that function as global regulators, such as the ferric uptake regulator (Fur), have also been involved in the regulation of T6SS loci [26, 33]. Moreover, although transcription of most T6SS genes depends on the housekeeping sigma factor σ70 (also known as RpoD or σA), other sigma factors modulate transcription of T6SS loci. This includes the alternative σ factor RpoS (also known as σ38 or σS), which belongs to the σ70 family of σ factors. RpoS regulates the bacterial general stress response and is activated in the stationary phase of growth [34]. Although RpoS does not appear to be a major regulatory player for T6SS expression, it has been shown to positively regulate the Klebsiella pneumoniae and Yersinia pseudotuberculosis T6SS loci [35, 36]. In the plant pathogen Xanthomonas citri, an alternative σ factor of the extracytoplasmic function sigma (σECF) factor group [37], the σEcfK factor, promotes T6SS expression during the interaction of the bacterium with its amoeba predator, Dictyostelium [38]. Further, RpoN (also known as σN or σ54), which is the only member of a different family of σ factors, the σ54 family [39], also modulates T6SS expression [26]. RpoN- dependent promoters require bacterial enhancerbinding proteins (bEBPs) to activate transcription initiation [40]. Interestingly, genes encoding these proteins are often found within the T6SS clusters. For example, vasH in the V. cholerae T6SS (Fig.1b) [2, 18, 41] and sfa2 and sfa3 encoded in the P. aeruginosa H2 and H3- T6SSs, respectively [15]. Curiously, the bEBP FleQ, which is a master regulator of flagella biosynthesis, has been associated with T6SS regulation in different species of Pseudomonas [42–44]. The biocontrol agent P. putida KT2440, an environmental strain able to promote the growth of crop plants and defend plants from the attack of phytopathogens, contains three T6SS gene clusters [5, 45, 46] (Fig.1b). We have recently reported that the K1- T6SS of this bacterium is an effective mechanism of biocontrol. K1- T6SS is a powerful molecular weapon used by P. putida to outcompete plant pathogens, including A. tumefaciens, Pseudomonas syringae, Xanthomonas campestris, Pectobacterium carotovorum and Pseudomonas savastanoi, among others [5, 11]. The structural components and associated effectors of the P. putida K1- T6SS have been studied in previous works [5, 10, 11]; however, the regulation of the system is still unknown. In this work, we have studied the transcriptional regulatory network that governs the P. putida K1- T6SS cluster. We have identified two transcriptional units within the K1- T6SS locus and four σ70- dependent promoters driving transcription of the K1- T6SS operons. We found that the expression of the system is induced at the stationary phase of growth by a still unknown mechanism that does not involve the stationary phase σ factor RpoS, which represses instead of induces expression of these genes. Analysis of K1- T6SS expression in several global regulators mutants showed that the system is also repressed by RpoN and its activator FleQ. A potential RpoN- binding box was found in the K1- T6SS promoter. However, mutation of this box showed that RpoN did not bind to this site. On the contrary, the K1- T6SS expression is activated by the GacS–GacA TCS, while its associated sensor protein RetS represses transcription of the system. Our results show that production of the P. putida K1- T6SS is governed by a complex regulatory network that in this way manages the biocontrol capabilities of this bacterium. METHODS Bacterial strains and growth conditions Bacterial strains are listed in Table S1 (available with the online version of this article). Unless otherwise stated, chemicals and reagents, including antibiotics, were purchased from Sigma- Aldrich and components to prepare growth media were purchased from Oxoid. All strains were grown in lysogeny broth (LB) (10 g l−1 NaCl) and agar (1.5 % w/v) [47] for routine growth with shaking at 200 r.p.m., as appropriate. E. coli strains were incubated at 37 °C and P. putida strains at 30 °C. Antibiotics were used at (µg ml−1): rifampicin (Rif), 20 for Pseudomonas; tetracycline (Tc), 20 for P. putida and 10 for E. coli; kanamycin (Km), 50 for P. putida and 25 for E. coli; ampicillin (Ap) 100 for E. coli; piperacillin (Pip), 25 for P. putida. Plasmids and cloning Plasmids and primers used in this study are listed in Tables S2 and S3. DNA manipulations were performed using standard methods [47]. KOD Hot Start DNA polymerase (Merck) or Phusion High- Fidelity DNA polymerase (New England Biolabs) were used for PCR reactions according to the manufacturer’s instructions. Primers were synthesized by Sigma- Aldrich and restriction enzymes were purchased from New England Biolabs. All DNA constructs were sequenced and verified to be correct before use. The 343, 507 and 456 bp DNA regions containing the putative K1- T6SS promoter regions PtagB1, Phcp1 and PvrgG1, respectively, were amplified from genomic DNA extracted from the P. putida KT2440 strain using primers P1–P6 (Table S3). Promoters were cloned into the broadhost- range, lowcopy- number vector pMP220 at the EcoRI/KpnI/PstI sites to produce transcriptional fusions to the lacZ gene. Recombinant plasmids were sequenced and transferred to P. putida by electroporation [48]. 4 Bernal etal., Microbiology 2023;169:001295 A DNA region upstream of PP3084 (173 bp) was amplified from genomic DNA extracted from the P. putida KT2440 strain using primers P7–P8 (Table S3) and cloned into pMP220 at EcoRI/XbaI sites to produce a transcriptional fusion to the lacZ gene. The gene PP3086 and its promoter region (711 bp) were amplified using primers P9–P10 (Table S3) and cloned into pMMB67EH at the EcoRI/HindIII sites. Recombinant plasmids (pMP220- based and pMMB67_PP3086) were sequenced and cotransferred to P. putida by electroporation [48]. β-galactosidase assays Overnight cultures of the P. putida strains carrying the pMP220 plasmid and derivates were diluted to a final turbidity (A600) of 0.05 in fresh LB medium containing tetracycline and cultures were grown at 30 °C and 200 r.p.m. for 4 and 24 h. Turbidity reached around 0.7 after 4 h and 5.5 after 24 h, and at this point aliquots were taken to measure βgalactosidase activity in permeabilized whole cells as described by Miller [49]. At least three independent assays were performed in each case, and standard errors of the means were calculated. The P. putida strains cotransformed with the plasmids derived from pMP220 and pMMB67 were handled in the same manner and the samples were prepared as explained above (exponential phase of growth). RNA purification P. putida KT2440 was grown at 30 °C in LB at 200 r.p.m. until turbidity reached approximately 0.7 (exponential phase cultures) or 5–6 (stationary phase cultures). Cells were harvested by centrifugation (8,000 g for 5 min at 4 °C) in disposable plastic tubes precooled in liquid nitrogen and the pellets were kept at −80 °C until use. Total RNA was extracted using the TRI reagent method (Ambion, ref. 9738, Austin, TX, USA) as recommended by the manufacturer, except that TriPure Isolation Reagent was preheated at 70 °C followed by purification with RNeasy columns (Qiagen). A final digestion step with RNasefree Ambion TURBO DNase was added at the end of the process. The RNA concentration was determined using a spectrophotometer and RNA integrity was assessed by 2 % agarose gel electrophoresis. Reverse transcription-PCR To determine the transcriptional unit of the K1- T6SS cluster, reverse transcription PCR was performed as previously described by Duque et al. [50] with primers P11–P36 (Table S3) based on adjacent genes. When the genes have a small size, the primers cover the intergenic regions between three adjacent genes to obtain an optimal PCR product in the range of 0.5 to 1.5 kb. We have used mRNA isolated from P. putida KT2440 cells in the stationary phase of growth as indicated above and the Roche One- Step RTPCR kit. Quantitative real-time PCR Total RNA was retrotranscribed to cDNA with Invitrogen Superscript II reverse transcriptase using random hexamers as primers. The specific primer pairs used in this study to amplify cDNA were: P37–P38 for the gene encoding Hcp1; P39–P40 for the gene encoding TssB1; P41–P42 for the gene encoding PP3084; and P43–P44 for the gene encoding 16S rRNA (Table S3). Realtime PCR was performed using the iCycler iQ Real- Time PCR Detection System (Bio- Rad). The target cDNA of the experimental and reference samples was amplified in triplicate. Each 25 µl reaction contained 2 µl of a dilution of target cDNA (1 : 10–1 : 10 000), 0.5 µl of each primer (20 µM), 9.5 µl milliQ water and 12.5 µl iQ SYBR Green Supermix (Bio- Rad, reference 170–8882). Samples were denatured by heating at 95 °C for 10 min prior to a 40- cycle amplification and quantification programme (95 °C for 15 s, 60 °C for 30 s and 72 °C for 20 s) with a single fluorescence measurement per cycle according to the manufacturer’s recommendations. A final 1 min extension cycle at 72 °C was performed. The PCR products were between 103 and 300 bp in length. To confirm the amplification of a single PCR product, a melting curve was obtained by slow heating from 55–95 °C at a rate of 0.5 °C every 10 s for 80 cycles with continuous fluorescence scanning. The results were normalized relative to those obtained for 16S rRNA. Quantification was based on the analysis of threshold cycle (Ct) values as described by Pfaffl [51]. Primer extension analysis Primer extension analyses were performed as previously described by Pacheco- Sánchez et al. [52]. An oligonucleotide complementary to the coding strand of the tagB1 gene (P45, Table S3) was 32P labelled at its 5′ ends in a 10 µl final volume that contained 1 µl 10× buffer, 10 pmol oligonucleotides, 1 µl [γ-32P]ATP (6000 mCi mmol−1), and 1 U phage T4 polynucleotide kinase. The reaction mixture was incubated for 1 h at 37 °C and 10 min at 70 °C to inactivate the kinase, and the labelled oligonucleotide was filtered through a Bio- Rad Micro Bio- Spin column to eliminate unbound nucleotide. Labelled primers were annealed to total RNA isolated as described above in a 10 µl annealing mixture that contained 2 µl 5× annealing buffer, 10 5 c.p.m. of 5′-end- labelled primer and 50 µg total RNA template. The mixture was heated at 95 °C for 3 min, incubated at 65 °C for 5 min, and then slowly cooled to 44 °C. cDNA was synthesized by the addition of 40 µl of reverse transcriptase buffer, 1 mM of dNTPs, 0.4 U µl−1 of RNase 5 Bernal etal., Microbiology 2023;169:001295 inhibitor and 8 U of AMV reverse transcriptase. The mixture was incubated for 1 h at 44 °C and the reaction was terminated by adding 5 µl 3 M sodium acetate and 150 µl ethanol. The product of reverse transcription was analysed in urea–polyacrylamide sequencing gels. The gel was exposed to a GS- 525 Molecular Imager (Bio- Rad). 5′ rapid amplification of cDNA Ends (5′-RACE) To identify the transcription start site (TSS) of the K1- T6SS promoters, 5′ rapid amplification of cDNA ends (RACE) [53] was performed on the tagB1 and PP3104 genes with the SMARTer RACE cDNA Amplification kit (cat. no. 634924) from Clontech. Following the manufacturer’s instructions, total RNA was ligated to an adapter at the 5′ end, and cDNA synthesis was initiated from an oligo(dT) primer containing a 3′ adapter sequence. Genespecific primers (P46 and P47, Table S3) were used in conjunction with the adapter primers provided to amplify 5′ and 3′ cDNA ends. The PCR products were extracted from a 1.5 % agarose gel and purified with the QIAquick PCR purification kit following the manufacturer’s instructions (Qiagen). The purified PCR products resulting from RACE were cloned into the pCR 2.1- TOPO vector (Invitrogen) as recommended. RACE products were sequenced by GATC Biotech services. Statistical analyses Statistical analyses are based on the ttest in which two conditions are compared independently. P- values from raw data were calculated by twotailed ttest and from ratio data to the control by onesample ttest using GraphPad Prism version 8.3.0 and are represented in the graphs by ns, nonsignificant; *P<0.05; **P<0.01; and ***P<0.001. RESULTS The K1-T6SS gene clusters are arranged in operons The 44 kb K1- T6SS gene cluster of P. putida KT2440 encoding structural and accessory components comprises two predicted operons that we previously named structural and vgrG operon, and an intermediate region [5]. The structural operon harbours a total of 15 genes, 12 of which encode the T6SS core components (TssABCDEFGHJKLM), while 3 encode accessory elements, including TagF1, TagP1 and the recently characterized TagB1 stabilizer protein [11]. The vgrG operon encompasses genes encoding proteins with unknown function (PP3104 and PP3105), the VgrG1 protein and associated adapters EagR1a and EagR1b, the Tke2 nuclease toxin and its cognate immunity protein Tki2, and an additional putative effector–immunity pair (Tke3–Tki3) (Fig.2a). The intermediate region mostly contains short ORFs that are not in the same transcriptional orientation [5], and thus it has not been considered for this study. Sequence analysis showed that most genes of the structural operon overlap and that the intergenic regions between the genes that do not overlap are smaller than 30 nucleotides (Table S4). This organization previously described by Normark et al. [54] indicates that these genes are arranged in a single transcriptional unit. Exceptions include the intergenic region between hcp1- tssA1 and tagP1- hcp1, which are 95 and 119 bp in length, respectively. The intergenic regions between the genes of the vgrG operon that do not overlap are also short, ranging from 0 to 25 nucleotides (Table S4), which suggests that these genes are also arranged in a single transcriptional unit. To test this prediction, we analysed the cotranscription of adjacent genes within the two clusters by RTPCR using primers designed every 1–1.5 kb (Fig. S1a and Table S3) and mRNA isolated from P. putida cells at the stationary phase. All sets of primers were functional as determined using genomic DNA (Fig. S1b). When using mRNA, cDNA fragments of the expected sizes were obtained with all pairs of primers in both the structural (Fig.2b, P11–28) and the vgrG (Fig.2b, P29–34) clusters. This indicates that these clusters are indeed arranged in single transcriptional units (Fig.2a). As expected, there was no amplification of a divergently transcribed region located outside the operons and used as a negative control (Fig.2b, P35–36). To map the promoters driving transcription of the structural (S) and the vgrG (V) operons, total RNA of P. putida cells in the stationary phase was analysed using RACE and primers complementary to tagB1 and PP3104 genes (Table S3). A single PCR band placed the transcriptional start site (+1) of the S operon in an adenine (A) located 121 bp upstream of the predicted translational start site of tagB1. This promoter was named P S (Fig.2c, left panel). Two PCR bands were obtained for the V operon, which mapped the+1 in a guanine (G) located 64 bp and a thymine (T) located 155 bp, respectively, upstream of the predicted translational start site of PP3104 (Fig.2c, right panel). These promoters were named PV1 and PV2. Unlike RACE, primer extension analysis of the Ps promoter also identified two different transcriptional start sites in this region (Fig.2d). The lower band positioned a +1 at the same adenine as the RACE analysis. The higher band mapped a second transcriptional start site in an adenine located 215 bp upstream of the tagB1 translational start site (Fig.2d). These two sequential promoters of the structural operon were named PS1 and PS2 (Fig.2). In silico analyses showed that the four identified promoter regions contain σ70 consensus sequences located at positions −35 and −10 from the+1 (Fig.2e). This suggests that transcription from these promoters depends on the σ70- loaded RNA polymerase (RNAP). 6 Bernal etal., Microbiology 2023;169:001295 Expression of K1-T6SS genes depends on the growth phase The K1- T6SS has previously been described as an active system under routine laboratory conditions and in the stationary phase of growth [5, 11], a condition in which other T6SS clusters have also been reported to be active [55, 56]. To analyse the effect of the growth phase on K1- T6SS gene expression, we performed quantitative reverse transcription PCR (qRT- PCR) using mRNA Fig. 2. Transcriptional organization of the K1- T6SS cluster. (a) The 44 kb chromosomal region containing the K1- T6SS genes is shown. Block arrows represent the different genes, their relative sizes and their transcriptional orientation, with the name of the gene or the PP number (http://www. pseudomonas.com/) indicated above the arrow. Genes encoding structural components are shown in black and genes encoding accessory proteins, effector and immunity proteins, and VgrG adaptors are represented in white. Genes encoding proteins of unknown function are depicted in grey. tssB1 and hcp1 genes from the S operon and PP3104 from the V operon are displayed in red, gold and green, respectively. The promoter regions identified or tested in this work are shown by coloured arrows above the genes. (b) RTPCR analysis to define the transcriptional units of the K1- T6SS cluster. mRNA from P. putida KT2440 wildtype strain grown until stationary phase was used. Primer pairs from P11 to P28 were used to amplify mRNA from the structural operon, P29 to P34 amplify the transcript of the vgrG operon. Primer pair P35–P36 was used as a negative control. Identification of the+1 site of the K1- T6SS structural and vgrR operons by 5′ RACE (c) and primer extension (d) analyses. P. putida KT2440 was grown in LB and samples were taken in stationary phase for total RNA isolation. (c) Pictures show cDNA products obtained by 5′-RACE analysis using P46 and P47 oligonucleotides (Table S3) in conjunction with adapter primers to amplify tagB1 and PP3104 5′cDNA ends, respectively. (d) The autoradiogram shows the cDNA products obtained after reverse transcription of 50 µg of total RNA with the 5′-end- labelled P45 oligonucleotide (Table S3) hybridizing with the tagB1 mRNA. (e) PS and PV promoter analyses. The P. putida KT2440 sequences corresponding to the regions upstream the tagB and PP3104 genes are shown. Nucleotides in bold, italic and underlined represent the+1 sites. The identified −10 and −35 σ70 binding sites are indicated. The sequences between the two promoters from the same region are in a contracted format as Nx where x represents the number of bases present between the two sequence segments. 7 Bernal etal., Microbiology 2023;169:001295 from P. putida cells grown until the exponential or stationary phase. Representative genes from both K1- T6SS operons were used in this analysis, i.e. hcp1 and tssB1 for the S operon, and PP3104 from the V operon. Expression of the three genes was induced considerably in the stationary phase compared to the exponential phase (12-, 5- and 2.6- fold, respectively) (Fig.3a). Interestingly, the qRT- PCR analysis showed that although hcp1 is located at the end of the S operon and tssB1 at the beginning, mRNA levels of hcp1 were higher than those of tssB1 (Fig.3a). This suggests that hcp1 could also be expressed from an internal promoter located upstream of hcp1 (Ph). To test this hypothesis and to analyse further the expression of the K1- T6SS promoters, we constructed transcriptional fusions of PS, PV and the putative Ph promoters to a promoterless lacZ gene using the pMP220 plasmid. The resulting plasmids, pMPS, pMPV and pMPh were introduced into the P. putida wildtype KT2440 strain and the β-galactosidase activity was determined Fig. 3. Expression of K1- T6SS genes. (a) mRNA levels of the indicated genes were obtained by qRT- PCR upon growth of the P. putida in LB until the exponential or stationary phase. The 2−ΔΔCT method was used to determine the fold change range in gene expression in the stationary versus the exponential phase. Data are the means±sd from three biological replicates (n=3), each one including three technical replicates. P- values were calculated by onesample ttest to a hypothetical value of 1 as described in the Methods. The horizontal dashed line denotes a fold change of 2. (b) β-galactosidase activity in the exponential and stationary phases of the P. putida wildtype strain bearing the pMP220- derived plasmids containing the indicated transcriptional fusion (PS, Ph and PV) to the lacZ gene. Strains were grown in LB until the exponential or stationary phase. Data are the means±sd from at least four replicates (n≥4), each one including two technical replicates. P- values were calculated by ttest analysis as described in the Methods. 8 Bernal etal., Microbiology 2023;169:001295 in the exponential and stationary phases of growth. pMPS contains the 345 bp upstream of the tagB1 translational start site that includes the previously identified PS1 and PS2 sequential promoters (Figs2 and 3). In agreement with the qRT- PCR results, transcription from this promoter was practically null in the exponential phase (71 Miller units, MU) and notably induced by more than 10- fold in the stationary phase (724 MU) (Fig.3b, left panel). Instead, the β-galactosidase activity of the pMP V construct that included the 458 bp upstream of the PP3104 translational start site containing the PV1 and PV2 promoters was high in both the exponential (~5000 MU) and the stationary phase (~9000 MU) (Fig.3b, right panel). This resulted in only 1.8- fold induction in the activity of this promoter in stationary compared to exponential phase, which was in accordance with the qRT- PCR results. The pMPh plasmid contains 509 bp upstream of the hcp1 translational start site, which includes 378 bp of the tagP1 3′ coding region. The β-galactosidase activity of this construct was negligible in both the exponential and the stationary phases (12 and 73 MU, respectively), and similar to that of the empty vector control (10.6 and 45 MU, respectively) (Fig.3b, middle panel). This indicates that there is not an internal promoter upstream of hcp1 and that the higher detection of the mRNA from this gene is likely due to the greater stability of this part of the transcript. Identification of global regulators involved in the control of K1-T6SS expression Several regulatory proteins have been involved in the control of the expression of T6SS clusters at different stages, including transcription, translational and posttranslational (assembling) levels. To identify the main transcriptional regulators involved in the control of the K1- T6SS cluster, we used a battery of P. putida globalregulator mutants (gacS, retS, turA, rpoS, rpoN and fleQ). We first focused on the Gac signalling cascade and the retS mutant. This mutant has been used extensively in P. aeruginosa to study the activity of the H1- T6SS because this system is not produced in the wildtype strain containing RetS [1, 57]. First, we compared the expression of hcp1, tssB1 and PP3104 genes in the P. putida wildtype versus the retS mutant in the exponential phase of growth by qRT- PCR. The three genes were significantly upregulated in the retS mutant, with foldchanges >3 for hcp1 and>1.5 for tssB1 and PP3104 (Fig.4a). Induction of the expression of the S operon in this mutant in the exponential phase of growth was confirmed by β-galactosidase using the PS::lacZ fusion (Fig.4b). This effect was not observed in the stationary phase (Fig.4b), which suggests that repression of the K1- T6SS by RetS occurs mainly during exponential growth. The activity of the PV promoter was not affected in the retS mutant in either the exponential or the stationary phase of growth (Fig.4c). However, activity of both PS and PV in the gacS response regulator mutant was considerably reduced (between two- and threefold) independently of the growth phase (Fig.4b, c). This indicates that GacS is a positive regulator of K1- T6SS expression, and this activity is repressed only by RetS in the exponential phase of growth. Further, transcription from PS and PV promoters stayed almost unchanged in the turA mutant (Fig.4b, c), indicating that this H- NS protein does not control K1- T6SS expression. Intriguingly, despite the fact that K1- T6SS expression is induced in the stationary phase of growth (Fig.3), the activity of both PS and PV promoters was not reduced but increased in absence of the stationary σ factor RpoS (Figs4b, c and S2). This is in accordance with the presence of σ70 but no RpoS consensus sequences in these promoter regions (Fig.2e). Strikingly, the transcriptional level from PS increased substantially in the absence of RpoN and FleQ (Fig.4b), especially in the exponential phase of growth, a condition in which the system is not expressed in the wildtype strain. This suggests that these regulators mostly repress expression from the PS promoter in exponential phase. The effect of RpoN is likely indirect since σ factors usually promote instead of inhibiting transcription. RpoN indirectly downregulates K1-T6SS expression We next aimed to analyse the effect of RpoN on the expression of the K1- T6SS S operon. This σ factor could either promote the transcription of an unknown T6SS repressor or directly repress the transcription from the PS promoter. RpoN recognizes and binds the 5′-TGGCAC- N5- TTGCW- 3′ (W is A/T) consensus sequence located at −12 and −24 bp from the transcriptional start site [39]. The minimal recognition sequence described to be RpoN- dependent is GG- N10- GC [39]. Interestingly, we identified a potential RpoN binding site in the PS2 promoter that overlaps with the σ70 binding site (Fig.5a). Because RpoN can stably interact with promoter sequences without initiating transcription [58], binding of this σ factor to PS2 could prevent the binding of the σ70 factor and therefore the σ70- mediated transcription from this promoter. It is important to note that in these cases, the socalled −12 and −24 boxes do not need to be located exactly in these positions from the transcriptional start site [58]. However, we have maintained this nomenclature of the boxes for clarity. To test whether RpoN binds to the PS2 promoter, we eliminated the conserved RpoN- binding GG and GC sequences in the PS2 fragment included in the pMP S plasmid. Because GG is part of both the −35 and −24 boxes, we produced two different constructs; M1 in which only the GC of the −12 box was changed to TT, and M2 in which both the GG and GC are changed to TT (Fig.5a). β-galactosidase assays showed that the activity of the M2 construct was completely null in both the wildtype and the rpoN mutant (Fig.5b), as expected from introducing changes in the σ70 -35 binding sequence. This indicates that this sequence in the PS2 promoter is crucial for expression of the structural operon. By contrast, the β-galactosidase activity of the M1 construct was similar to that of the wildtype promoter (Fig.5b). This indicates that RpoN does not bind and repress the expression from 9 Bernal etal., Microbiology 2023;169:001295 Fig. 4. Global regulators involved in K1- T6SS expression. (a) mRNA levels of the indicated genes were obtained by qRT- PCR upon growth of the P. putida and its isogenic retS mutant in LB until the exponential phase of growth. The 2−ΔΔCT method was used to determine the fold change range in gene expression in the retS versus the wildtype strain. Data are means±sd from three biological replicates (n=3), each one including three technical replicates. P- values were calculated by onesample ttest to a hypothetical value of 1 as described in the Methods. The horizontal dashed line denotes a fold change of 2. The activity of the PS (b) and PV (c) promoters in the indicated P. putida mutant strain was obtained by β-galactosidase assay. Strains bearing the indicated lacZ fusion were grown in LB until exponential or stationary phase. Data show the fold change of the promoter activity in the mutants relative to the KT2440 wildtype strain. Data are means±sd from at least four biological replicates (n≥4), each one including two technical replicates. P- values were calculated by onesample ttest to a hypothetical value of 1 as described in the Methods. The horizontal dashed lines denote a fold change of 1 in either direction represented in a log2 scale. 16 Bernal etal., Microbiology 2023;169:001295 63. Gonzalez D, Mavridou DAI. Making the best of aggression: the many dimensions of bacterial toxin regulation. Trends in Microbiology 2019;27:897–905. 64. Molina L, Bernal P, Udaondo Z, Segura A, Ramos J- L. Complete Genome Sequence of a Pseudomonas putida Clinical Isolate, Strain H8234. Genome Announc 2013;1:e00496- 13. 65. Molina L, Udaondo Z, Duque E, Fernández M, Bernal P, et al. Specific gene loci of clinical Pseudomonas putida isolates. PLoS One 2016;11:e0147478. 66. Borrero de Acuña JM, Bernal P. Plant holobiont interactions mediated by the type VI secretion system and the membrane vesicles: promising tools for a greener agriculture. Environ Microbiol 2021;23:1830–1836. 67. Weber BS, Ly PM, Irwin JN, Pukatzki S, Feldman MF. A multidrug resistance plasmid contains the molecular switch for type VI secretion in Acinetobacter baumannii. Proc Natl Acad Sci U S A 2015;112:9442–9447. 68. Ramos C, Mølbak L, Molin S. Bacterial activity in the rhizosphere analyzed at the singlecell level by monitoring ribosome contents and synthesis rates. Appl Environ Microbiol 2000;66:801–809. 69. Hug S, Liu Y, Heiniger B, Bailly A, Ahrens CH, etal. Differential expression of Paraburkholderia phymatum type VI secretion systems (T6SS) suggests a role of T6SS- b in early symbiotic interaction. Front Plant Sci 2021;12:69590. 70. Jiang F, Waterfield NR, Yang J, Yang G, Jin Q. A Pseudomonas aeruginosa type VI secretion phospholipase D effector targets both prokaryotic and eukaryotic cells. Cell Host Microbe 2014;15:600–610. 71. Pini C, Godoy P, Bernal P, Ramos J, Segura A. Regulation of the cyclopropane synthase cfaB gene in Pseudomonas putida KT2440. FEMS Microbiol Lett 2011;321:107–114. 72. Ramos- González MI, Molin S. Cloning, sequencing, and phenotypic characterization of the rpoS gene from Pseudomonas putida KT2440. J Bacteriol 1998;180:3421–3431. 73. Schuster M, Hawkins AC, Harwood CS, Greenberg EP. The Pseudomonas aeruginosa RpoS regulon and its relationship to quorum sensing. Mol Microbiol 2004;51:973–985. 74. Heurlier K, Dénervaud V, Pessi G, Reimmann C, Haas D. Negative control of quorum sensing by RpoN (sigma54) in Pseudomonas aeruginosa PAO1. J Bacteriol 2003;185:2227–2235. 75. Baraquet C, Harwood CS. Cyclic diguanosine monophosphate represses bacterial flagella synthesis by interacting with the walker a motif of the enhancerbinding protein FleQ. Proc Natl Acad Sci 2013;110:18478–18483. 76. Leal- Morales A, Pulido- Sánchez M, López- Sánchez A, Govantes F. Transcriptional organization and regulation of the Pseudomonas putida flagellar system. Environ Microbiol 2022;24:137–157. 77. Nie H, Xiao Y, Song M, Wu N, Peng Q, et al. Wsp system oppositely modulates antibacterial activity and biofilm formation via FleQ- FleN complex in Pseudomonas putida. Environ Microbiol 2022;24:1543–1559. 78. Ferreiro MD, Gallegos MT. Distinctive features of the Gac- Rsm pathway in plantassociated Pseudomonas. Environ Microbiol 2021;23:5670–5689. 79. Wang BX, Wheeler KM, Cady KC, Lehoux S, Cummings RD, etal. Mucin glycans signal through the sensor kinase rets to inhibit virulenceassociated traits in Pseudomonas aeruginosa. Current Biology 2021;31:90–102. 80. Records AR, Gross DC. Sensor kinases RetS and LadS regulate Pseudomonas syringae type VI secretion and virulence factors. J Bacteriol 2010;192:3584–3596. 81. Heeb S, Haas D. Regulatory roles of the GacS/GacA twocomponent system in plantassociated and other gramnegative bacteria. Mol Plant Microbe Interact 2001;14:1351–1363. Edited by: D. Mavridou Five reasons to publish your next article with a Microbiology Society journal 1. When you submit to our journals, you are supporting Society activities for your community. 2. Experience a fair, transparent process and critical, constructive review. 3. If you are at a Publish and Read institution, you’ll enjoy the benefits of Open Access across our journal portfolio. 4. Author feedback says our Editors are ‘thorough and fair’ and ‘patient and caring’. 5. Increase your reach and impact and share your research more widely. Find out more and submit your article at microbiologyresearch.org.