Paraburkholderia sabiae Uses One Type VI Secretion System (T6SS-1) as a Powerful Weapon against Notorious Plant Pathogens

Paraburkholderia sabiae LMG24235 is a nitrogen-fixing betaproteobacterium originally isolated from a root nodule of Mimosa caesalpiniifolia in Brazil. We show here that this strain effectively kills strains from several bacterial families (Burkholderiaceae, Pseudomonadaceae, Enterobacteriaceae) which include important plant pathogens in a contact-dependent manner. De novo assembly of the first complete genome of P. sabiae using long sequencing reads and subsequent annotation revealed two gene clusters predicted to encode type VI secretion systems (T6SS), which we named T6SS-1 and T6SS-3 according to previous classification methods (G. Shalom, J. G. Shaw, and M. S. Thomas, Microbiology, 153:2689-2699, 2007, https://doi.org/10.1099/mic.0.2007/006585-0). We created P. sabiae with mutations in each of the two T6SS gene clusters that abrogated their function, and the T6SS-1 mutant was no longer able to outcompete other strains in a contact-dependent manner. Notably, our analysis revealed that T6SS-1 is essential for competition against several important plant pathogens in vitro, including Burkholderia plantarii, Ralstonia solanacearum, Pseudomonas syringae, and Pectobacterium carotovorum. The 9-log reduction in P. syringae cells in the presence of P. sabiae was particularly remarkable. Importantly, in an in vivo assay, P. sabiae was able to protect potato tubers from bacterial soft rot disease caused by P. carotovorum, and this protection was partly dependent on T6SS-1. IMPORTANCE Rhizobia often display additional beneficial traits such as the production of plant hormones and the acquisition of limited essential nutrients that improve plant growth and enhance plant yields. Here, we show that the rhizobial strain P. sabiae antagonizes important phytopathogens such as P. carotovorum, P. syringae, and R. solanacearum and that this effect is due to contact-dependent killing mediated by one of two T6SS systems identified in the complete, de novo assembled genome sequence of P. sabiae. Importantly, co-inoculation of Solanum tuberosum tubers with P. sabiae also resulted in a drastic reduction of soft rot caused by P. carotovorum in an in vivo model system. This result highlights the protective potential of P. sabiae against important bacterial plant diseases, which makes it a valuable candidate for application as a biocontrol agent. It also emphasizes the particular potential of rhizobial inoculants that combine several beneficial effects such as plant growth promotion and biocontrol for sustainable agriculture.

Detection of Known and Novel Small Proteins in Pseudomonas stutzeri Using a Combination of Bottom-Up and Digest-Free Proteomics and Proteogenomics

Small proteins of around 50 aa in length have been largely overlooked in genetic and biochemical assays due to the inherent challenges with detecting and characterizing them. Recent discoveries of their critical roles in many biological processes have led to an increased recognition of the importance of small proteins for basic research and as potential new drug targets. One example is CcoM, a 36 aa subunit of the cbb3-type oxidase that plays an essential role in adaptation to oxygen-limited conditions in Pseudomonas stutzeri (P. stutzeri), a model for the clinically relevant, opportunistic pathogen Pseudomonas aeruginosa. However, as no comprehensive data were available in P. stutzeri, we devised an integrated, generic approach to study small proteins more systematically. Using the first complete genome as basis, we conducted bottom-up proteomics analyses and established a digest-free, direct-sequencing proteomics approach to study cells grown under aerobic and oxygen-limiting conditions. Finally, we also applied a proteogenomics pipeline to identify missed protein-coding genes. Overall, we identified 2921 known and 29 novel proteins, many of which were differentially regulated. Among 176 small proteins 16 were novel. Direct sequencing, featuring a specialized precursor acquisition scheme, exhibited advantages in the detection of small proteins with higher (up to 100%) sequence coverage and more spectral counts, including sequences with high proline content. Three novel small proteins, uniquely identified by direct sequencing and not conserved beyond P. stutzeri, were predicted to form an operon with a conserved protein and may represent de novo genes. These data demonstrate the power of this combined approach to study small proteins in P. stutzeri and show its potential for other prokaryotes.

Unraveling the small proteome of the plant symbiont Sinorhizobium meliloti by ribosome profiling and proteogenomics

The soil-dwelling plant symbiont Sinorhizobium meliloti is a major model organism of Alphaproteobacteria. Despite numerous detailed OMICS studies, information about small open reading frame (sORF)-encoded proteins (SEPs) is largely missing, because sORFs are poorly annotated and SEPs are hard to detect experimentally. However, given that SEPs can fulfill important functions, identification of translated sORFs is critical for analyzing their roles in bacterial physiology. Ribosome profiling (Ribo-seq) can detect translated sORFs with high sensitivity, but is not yet routinely applied to bacteria because it must be adapted for each species. Here, we established a Ribo-seq procedure for S. meliloti 2011 based on RNase I digestion and detected translation for 60% of the annotated coding sequences during growth in minimal medium. Using ORF prediction tools based on Ribo-seq data, subsequent filtering, and manual curation, the translation of 37 non-annotated sORFs with ≤ 70 amino acids was predicted with confidence. The Ribo-seq data were supplemented by mass spectrometry (MS) analyses from three sample preparation approaches and two integrated proteogenomic search database (iPtgxDB) types. Searches against standard and 20-fold smaller Ribo-seq data-informed custom iPtgxDBs confirmed 47 annotated SEPs and identified 11 additional novel SEPs. Epitope tagging and Western blot analysis confirmed the translation of 15 out of 20 SEPs selected from the translatome map. Overall, by combining MS and Ribo-seq approaches, the small proteome of S. meliloti was substantially expanded by 48 novel SEPs. Several of them are part of predicted operons and/or are conserved from Rhizobiaceae to Bacteria, suggesting important physiological functions.

A novel function of the key nitrogen-fixation activator NifA in beta-rhizobia: Repression of bacterial auxin synthesis during symbiosis

Rhizobia fix nitrogen within root nodules of host plants where nitrogenase expression is strictly controlled by its key regulator NifA. We recently discovered that in nodules infected by the beta-rhizobial strain Paraburkholderia phymatum STM815, NifA controls expression of two bacterial auxin synthesis genes. Both the iaaM and iaaH transcripts, as well as the metabolites indole-acetamide (IAM) and indole-3-acetic acid (IAA) showed increased abundance in nodules occupied by a nifA mutant compared to wild-type nodules. Here, we document the structural changes that a P. phymatum nifA mutant induces in common bean (Phaseolus vulgaris) nodules, eventually leading to hypernodulation. To investigate the role of the P. phymatum iaaMH genes during symbiosis, we monitored their expression in presence and absence of NifA over different stages of the symbiosis. The iaaMH genes were found to be under negative control of NifA in all symbiotic stages. While a P. phymatum iaaMH mutant produced the same number of nodules and nitrogenase activity as the wild-type strain, the nifA mutant produced more nodules than the wild-type that clustered into regularly-patterned root zones. Mutation of the iaaMH genes in a nifA mutant background reduced the presence of these nodule clusters on the root. We further show that the P. phymatum iaaMH genes are located in a region of the symbiotic plasmid with a significantly lower GC content and exhibit high similarity to two genes of the IAM pathway often used by bacterial phytopathogens to deploy IAA as a virulence factor. Overall, our data suggest that the increased abundance of rhizobial auxin in the non-fixing nifA mutant strain enables greater root infection rates and a role for bacterial auxin production in the control of early stage symbiotic interactions.

Differential Expression of Paraburkholderia phymatum Type VI Secretion Systems (T6SS) Suggests a Role of T6SS-b in Early Symbiotic Interaction

Paraburkholderia phymatum STM815, a rhizobial strain of the Burkholderiaceae family, is able to nodulate a broad range of legumes including the agriculturally important Phaseolus vulgaris (common bean). P. phymatum harbors two type VI Secretion Systems (T6SS-b and T6SS-3) in its genome that contribute to its high interbacterial competitiveness in vitro and in infecting the roots of several legumes. In this study, we show that P. phymatum T6SS-b is found in the genomes of several soil-dwelling plant symbionts and that its expression is induced by the presence of citrate and is higher at 20/28°C compared to 37°C. Conversely, T6SS-3 shows homologies to T6SS clusters found in several pathogenic Burkholderia strains, is more prominently expressed with succinate during stationary phase and at 37°C. In addition, T6SS-b expression was activated in the presence of germinated seeds as well as in P. vulgaris and Mimosa pudica root nodules. Phenotypic analysis of selected deletion mutant strains suggested a role of T6SS-b in motility but not at later stages of the interaction with legumes. In contrast, the T6SS-3 mutant was not affected in any of the free-living and symbiotic phenotypes examined. Thus, P. phymatum T6SS-b is potentially important for the early infection step in the symbiosis with legumes.

Genome-Wide Association Study Reveals Complex Genetic Architecture of Cadmium and Mercury Accumulation and Tolerance Traits in Medicago truncatula

Heavy metals are an increasing problem due to contamination from human sources that and can enter the food chain by being taken up by plants. Understanding the genetic basis of accumulation and tolerance in plants is important for reducing the uptake of toxic metals in crops and crop relatives, as well as for removing heavy metals from soils by means of phytoremediation. Following exposure of Medicago truncatula seedlings to cadmium (Cd) and mercury (Hg), we conducted a genome-wide association study using relative root growth (RRG) and leaf accumulation measurements. Cd and Hg accumulation and RRG had heritability ranging 0.44 - 0.72 indicating high genetic diversity for these traits. The Cd and Hg trait associations were broadly distributed throughout the genome, indicated the traits are polygenic and involve several quantitative loci. For all traits, candidate genes included several membrane associated ATP-binding cassette transporters, P-type ATPase transporters, oxidative stress response genes, and stress related UDP-glycosyltransferases. The P-type ATPase transporters and ATP-binding cassette protein-families have roles in vacuole transport of heavy metals, and our findings support their wide use in physiological plant responses to heavy metals and abiotic stresses. We also found associations between Cd RRG with the genes CAX3 and PDR3, two linked adjacent genes, and leaf accumulation of Hg associated with the genes NRAMP6 and CAX9. When plant genotypes with the most extreme phenotypes were compared, we found significant divergence in genomic regions using population genomics methods that contained metal transport and stress response gene ontologies. Several of these genomic regions show high linkage disequilibrium (LD) among candidate genes suggesting they have evolved together. Minor allele frequency (MAF) and effect size of the most significant SNPs was negatively correlated with large effect alleles being most rare. This is consistent with purifying selection against alleles that increase toxicity and abiotic stress. Conversely, the alleles with large affect that had higher frequencies that were associated with the exclusion of Cd and Hg. Overall, macroevolutionary conservation of heavy metal and stress response genes is important for improvement of forage crops by harnessing wild genetic variants in gene banks such as the Medicago HapMap collection.