Seed microbiota revealed by a large-scale meta-analysis including 50 plant species.

Seed microbiota constitutes a primary inoculum for plants that is gaining attention owing to its role for plant health and productivity. Here, we performed a meta-analysis on 63 seed microbiota studies covering 50 plant species to synthesize knowledge on the diversity of this habitat. Seed microbiota are diverse and extremely variable, with taxa richness varying from one to thousands of taxa. Hence, seed microbiota presents a variable (i.e. flexible) microbial fraction but we also identified a stable (i.e. core) fraction across samples. Around 30 bacterial and fungal taxa are present in most plant species and in samples from all over the world. Core taxa, such as Pantoea agglomerans, Pseudomonas viridiflava, P. fluorescens, Cladosporium perangustum and Alternaria sp., are dominant seed taxa. The characterization of the core and flexible seed microbiota provided here will help uncover seed microbiota roles for plant health and design effective microbiome engineering.

Maternal effects shape the seed mycobiome in Quercus petraea.

The tree seed mycobiome has received little attention despite its potential role in forest regeneration and health. The aim of the present study was to analyze the processes shaping the composition of seed fungal communities in natural forests as seeds transition from the mother plant to the ground for establishment. We used metabarcoding approaches and confocal microscopy to analyze the fungal communities of seeds collected in the canopy and on the ground in four natural populations of sessile oak (Quercus petraea). Ecological processes shaping the seed mycobiome were inferred using joint species distribution models. Fungi were present in seed internal tissues, including the embryo. The seed mycobiome differed among oak populations and trees within the same population. Its composition was largely influenced by the mother, with weak significant environmental influences. The models also revealed several probable interactions among fungal pathogens and mycoparasites. Our results demonstrate that maternal effects, environmental filtering and biotic interactions all shape the seed mycobiome of sessile oak. They provide a starting point for future research aimed at understanding how maternal genes and environments interact to control the vertical transmission of fungal species that could then influence seed dispersal and germination, and seedling recruitment.

rpoB, a promising marker for analyzing the diversity of bacterial communities by amplicon sequencing.

BACKGROUND: Microbiome composition is frequently studied by the amplification and high-throughput sequencing of specific molecular markers (metabarcoding). Various hypervariable regions of the 16S rRNA gene are classically used to estimate bacterial diversity, but other universal bacterial markers with a finer taxonomic resolution could be employed. We compared specificity and sensitivity between a portion of the rpoB gene and the V3 V4 hypervariable region of the 16S rRNA gene. RESULTS: We first designed universal primers for rpoB suitable for use with Illumina sequencing-based technology and constructed a reference rpoB database of 45,000 sequences. The rpoB and V3 V4 markers were amplified and sequenced from (i) a mock community of 19 bacterial strains from both Gram-negative and Gram-positive lineages; (ii) bacterial assemblages associated with entomopathogenic nematodes. In metabarcoding analyses of mock communities with two analytical pipelines (FROGS and DADA2), the estimated diversity captured with the rpoB marker resembled the expected composition of these mock communities more closely than that captured with V3 V4. The rpoB marker had a higher level of taxonomic affiliation, a higher sensitivity (detection of all the species present in the mock communities), and a higher specificity (low rates of spurious OTU detection) than V3 V4. We compared the performance of the rpoB and V3 V4 markers in an animal ecosystem model, the infective juveniles of the entomopathogenic nematode Steinernema glaseri carrying the symbiotic bacteria Xenorhabdus poinarii. Both markers showed the bacterial community associated with this nematode to be of low diversity (< 50 OTUs), but only rpoB reliably detected the symbiotic bacterium X. poinarii. CONCLUSIONS: Our results confirm that different microbiota composition data may be obtained with different markers. We found that rpoB was a highly appropriate marker for assessing the taxonomic structure of mock communities and the nematode microbiota. Further studies on other ecosystems should be considered to evaluate the universal usefulness of the rpoB marker. Our data highlight two crucial elements that should be taken into account to ensure more reliable and accurate descriptions of microbial diversity in high-throughput amplicon sequencing analyses: i) the need to include mock communities as controls; ii) the advantages of using a multigenic approach including at least one housekeeping gene (rpoB is a good candidate) and one variable region of the 16S rRNA gene. This study will be useful to the growing scientific community describing bacterial communities by metabarcoding in diverse ecosystems.

Crop Residues in Wheat-Oilseed Rape Rotation System: a Pivotal, Shifting Platform for Microbial Meetings.

Crop residues are a crucial ecological niche with a major biological impact on agricultural ecosystems. In this study, we used a combined diachronic and synchronic field experiment based on wheat-oilseed rape rotations to test the hypothesis that plant is a structuring factor of microbial communities in crop residues, and that this effect decreases over time with their likely progressive degradation and colonisation by other microorganisms. We characterised an entire fungal and bacterial community associated with 150 wheat and oilseed rape residue samples at a plurennial scale by metabarcoding. The impact of plant species on the residue microbiota decreased over time and our data revealed turnover, with the replacement of oligotrophs, often plant-specific genera (such as pathogens) by copiotrophs, belonging to more generalist genera. Within a single cropping season, the plant-specific genera and species were gradually replaced by taxa that are likely to originate from the soil. These changes occurred more rapidly for bacteria than for fungi, known to degrade complex compounds. Overall, our findings suggest that crop residues constitute a key fully-fledged microbial ecosystem. Taking into account this ecosystem, that has been neglected for too long, is essential, not only to improve the quantitative management of residues, the presence of which can be detrimental to crop health, but also to identify groups of beneficial microorganisms. Our findings are of particular importance, because the wheat-oilseed rape rotation, in which no-till practices are frequent, is particularly widespread in the European arable cropping systems.

Complete genomic sequence of Raphanus sativus cryptic virus 4 (RsCV4), a novel alphapartitivirus from radish.

The present work reports the discovery and complete genome sequencing of a virus from symptomless radish seedlings, classifiable as a novel member of the genus Alphapartitivirus, family Partitiviridae. Total RNA extracted from germinating seedlings was sequenced using Illumina technology. Bioinformatic analysis of the RNA-seq data revealed two contigs representing the near full-length genomic sequences of two genomic RNAs representing a new virus. Analysis of the genome sequence (excluding the polyA tail, RNA1: 1976 nt and RNA2: 1751 nt, respectively) showed a genomic organization typical of viruses classed within the Partitiviridae, with each genomic RNA encoding a single open reading frame (ORF). Phylogenetic analysis of the RNA dependent RNA polymerase (RNA1 ORF) and of the capsid protein (RNA2 ORF) clearly showed the new virus can be classified within the genus Alphapartitivirus, but sequence divergence establishes it as a new species, for which the name "Raphanus sativus cryptic virus 4" is proposed.

The severity of nonalcoholic fatty liver disease is associated with gut dysbiosis and shift in the metabolic function of the gut microbiota.

UNLABELLED: Several animal studies have emphasized the role of gut microbiota in nonalcoholic fatty liver disease (NAFLD). However, data about gut dysbiosis in human NAFLD remain scarce in the literature, especially studies including the whole spectrum of NAFLD lesions. We aimed to evaluate the association between gut dysbiosis and severe NAFLD lesions, that is, nonalcoholic steatohepatitis (NASH) and fibrosis, in a well-characterized population of adult NAFLD. Fifty-seven patients with biopsy-proven NAFLD were enrolled. Taxonomic composition of gut microbiota was determined using 16S ribosomal RNA gene sequencing of stool samples. Thirty patients had F0/F1 fibrosis stage at liver biopsy (10 with NASH), and 27 patients had significant F≥2 fibrosis (25 with NASH). Bacteroides abundance was significantly increased in NASH and F≥2 patients, whereas Prevotella abundance was decreased. Ruminococcus abundance was significantly higher in F≥2 patients. By multivariate analysis, Bacteroides abundance was independently associated with NASH and Ruminococcus with F≥2 fibrosis. Stratification according to the abundance of these two bacteria generated three patient subgroups with increasing severity of NAFLD lesions. Based on imputed metagenomic profiles, Kyoto Encyclopedia of Genes and Genomes pathways significantly related to NASH and fibrosis F≥2 were mostly related to carbohydrate, lipid, and amino acid metabolism. CONCLUSION: NAFLD severity associates with gut dysbiosis and a shift in metabolic function of the gut microbiota. We identified Bacteroides as independently associated with NASH and Ruminococcus with significant fibrosis. Thus, gut microbiota analysis adds information to classical predictors of NAFLD severity and suggests novel metabolic targets for pre-/probiotics therapies.

Ancestral acquisitions, gene flow and multiple evolutionary trajectories of the type three secretion system and effectors in Xanthomonas plant pathogens.

Deciphering the evolutionary history and transmission patterns of virulence determinants is necessary to understand the emergence of novel pathogens. The main virulence determinant of most pathogenic proteobacteria is the type three secretion system (T3SS). The Xanthomonas genus includes bacteria responsible for numerous epidemics in agroecosystems worldwide and represents a major threat to plant health. The main virulence factor of Xanthomonas is the Hrp2 family T3SS; however, this system is not conserved in all strains and it has not been previously determined whether the distribution of T3SS in this bacterial genus has resulted from losses or independent acquisitions. Based on comparative genomics of 82 genome sequences representing the diversity of the genus, we have inferred three ancestral acquisitions of the Hrp2 cluster during Xanthomonas evolution followed by subsequent losses in some commensal strains and re-acquisition in some species. While mutation was the main force driving polymorphism at the gene level, interspecies homologous recombination of large fragments expanding through several genes shaped Hrp2 cluster polymorphism. Horizontal gene transfer of the entire Hrp2 cluster also occurred. A reduced core effectome composed of xopF1, xopM, avrBs2 and xopR was identified that may allow commensal strains overcoming plant basal immunity. In contrast, stepwise accumulation of numerous type 3 effector genes was shown in successful pathogens responsible for epidemics. Our data suggest that capacity to intimately interact with plants through T3SS would be an ancestral trait of xanthomonads. Since its acquisition, T3SS has experienced a highly dynamic evolutionary history characterized by intense gene flux between species that may reflect its role in host adaptation.

Terroir is a key driver of seed-associated microbial assemblages.

Seeds have evolved in association with diverse microbial assemblages that may influence plant growth and health. However, little is known about the composition of seed-associated microbial assemblages and the ecological processes shaping their structures. In this work, we monitored the relative influence of the host genotypes and terroir on the structure of the seed microbiota through metabarcoding analysis of different microbial assemblages associated to five different bean cultivars harvested in two distinct farms. Overall, few bacterial and fungal operational taxonomic units (OTUs) were conserved across all seed samples. The lack of shared OTUs between samples is explained by a significant effect of the farm site on the structure of microbial assemblage, which explained 12.2% and 39.7% of variance in bacterial and fungal diversity across samples. This site-specific effect is reflected by the significant enrichment of 70 OTUs in Brittany and 88 OTUs in Luxembourg that lead to differences in co-occurrence patterns. In contrast, variance in microbial assemblage structure was not explained by host genotype. Altogether, these results suggest that seed-associated microbial assemblage is determined by niche-based processes and that the terroir is a key driver of these selective forces.

Emergence shapes the structure of the seed microbiota.

Seeds carry complex microbial communities, which may exert beneficial or deleterious effects on plant growth and plant health. To date, the composition of microbial communities associated with seeds has been explored mainly through culture-based diversity studies and therefore remains largely unknown. In this work, we analyzed the structures of the seed microbiotas of different plants from the family Brassicaceae and their dynamics during germination and emergence through sequencing of three molecular markers: the ITS1 region of the fungal internal transcribed spacer, the V4 region of 16S rRNA gene, and a species-specific bacterial marker based on a fragment of gyrB. Sequence analyses revealed important variations in microbial community composition between seed samples. Moreover, we found that emergence strongly influences the structure of the microbiota, with a marked reduction of bacterial and fungal diversity. This shift in the microbial community composition is mostly due to an increase in the relative abundance of some bacterial and fungal taxa possessing fast-growing abilities. Altogether, our results provide an estimation of the role of the seed as a source of inoculum for the seedling, which is crucial for practical applications in developing new strategies of inoculation for disease prevention.

Pairwise transcriptomic analysis of the interactions between the ectomycorrhizal fungus Laccaria bicolor S238N and three beneficial, neutral and antagonistic soil bacteria.

Ectomycorrhizal fungi are surrounded by bacterial communities with which they interact physically and metabolically during their life cycle. These bacteria can have positive or negative effects on the formation and the functioning of ectomycorrhizae. However, relatively little is known about the mechanisms by which ectomycorrhizal fungi and associated bacteria interact. To understand how ectomycorrhizal fungi perceive their biotic environment and the mechanisms supporting interactions between ectomycorrhizal fungi and soil bacteria, we analysed the pairwise transcriptomic responses of the ectomycorrhizal fungus Laccaria bicolor (Basidiomycota: Agaricales) when confronted with beneficial, neutral or detrimental soil bacteria. Comparative analyses of the three transcriptomes indicated that the fungus reacted differently to each bacterial strain. Similarly, each bacterial strain produced a specific and distinct response to the presence of the fungus. Despite these differences in responses observed at the gene level, we found common classes of genes linked to cell-cell interaction, stress response and metabolic processes to be involved in the interaction of the four microorganisms.

Seedling microbiota engineering using bacterial synthetic community inoculation on seeds.

Synthetic Communities (SynComs) are being developed and tested to manipulate plant microbiota and improve plant health. To date, only few studies proposed the use of SynCom on seed despite its potential for plant microbiota engineering. We developed and presented a simple and effective seedling microbiota engineering method using SynCom inoculation on seeds. The method was successful using a wide diversity of SynCom compositions and bacterial strains that are representative of the common bean seed microbiota. First, this method enables the modulation of seed microbiota composition and community size. Then, SynComs strongly outcompeted native seed and potting soil microbiota and contributed on average to 80% of the seedling microbiota. We showed that strain abundance on seed was a main driver of an effective seedling microbiota colonization. Also, selection was partly involved in seed and seedling colonization capacities since strains affiliated to Enterobacteriaceae and Erwiniaceae were good colonizers while Bacillaceae and Microbacteriaceae were poor colonizers. Additionally, the engineered seed microbiota modified the recruitment and assembly of seedling and rhizosphere microbiota through priority effects. This study shows that SynCom inoculation on seeds represents a promising approach to study plant microbiota assembly and its consequence on plant fitness.

The type VI secretion system of Stenotrophomonas rhizophila CFBP13503 limits the transmission of Xanthomonas campestris pv. campestris 8004 from radish seeds to seedlings.

Stenotrophomonas rhizophila CFBP13503 is a seedborne commensal bacterial strain, which is efficiently transmitted to seedlings and can outcompete the phytopathogenic bacterium Xanthomonas campestris pv. campestris (Xcc8004). The type VI secretion system (T6SS), an interference contact-dependent mechanism, is a critical component of interbacterial competition. The involvement of the T6SS of S. rhizophila CFBP13503 in the inhibition of Xcc8004 growth and seed-to-seedling transmission was assessed. The T6SS cluster of S. rhizophila CFBP13503 and nine putative effectors were identified. Deletion of two T6SS structural genes, hcp and tssB, abolished the competitive advantage of S. rhizophila against Xcc8004 in vitro. The population sizes of these two bacterial species were monitored in seedlings after inoculation of radish seeds with mixtures of Xcc8004 and either S. rhizophila wild-type (wt) strain or isogenic hcp mutant. A significant decrease in the population size of Xcc8004 was observed during confrontation with the S. rhizophila wt in comparison with T6SS-deletion mutants in germinated seeds and seedlings. We found that the T6SS distribution among 835 genomes of the Stenotrophomonas genus is scarce. In contrast, in all available S. rhizophila genomes, T6SS clusters are widespread and mainly belong to the T6SS group i4. In conclusion, the T6SS of S. rhizophila CFBP13503 is involved in the antibiosis against Xcc8004 and reduces seedling transmission of Xcc8004 in radish. The distribution of this T6SS cluster in the S. rhizophila complex could make it possible to exploit these strains as biocontrol agents against X. campestris pv. campestris.

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

BACKGROUND: Pseudomonas fluorescens F113 is a plant growth-promoting rhizobacterium (PGPR) isolated from the sugar-beet rhizosphere. This bacterium has been extensively studied as a model strain for genetic regulation of secondary metabolite production in P. fluorescens, as a candidate biocontrol agent against phytopathogens, and as a heterologous host for expression of genes with biotechnological application. The F113 genome sequence and annotation has been recently reported. RESULTS: Comparative analysis of 50 genome sequences of strains belonging to the P. fluorescens group has revealed the existence of five distinct subgroups. F113 belongs to subgroup I, which is mostly composed of strains classified as P. brassicacearum. The core genome of these five strains is highly conserved and represents approximately 76% of the protein-coding genes in any given genome. Despite this strong conservation, F113 also contains a large number of unique protein-coding genes that encode traits potentially involved in the rhizocompetence of this strain. These features include protein coding genes required for denitrification, diterpenoids catabolism, motility and chemotaxis, protein secretion and production of antimicrobial compounds and insect toxins. CONCLUSIONS: The genome of P. fluorescens F113 is composed of numerous protein-coding genes, not usually found together in previously sequenced genomes, which are potentially decisive during the colonisation of the rhizosphere and/or interaction with other soil organisms. This includes genes encoding proteins involved in the production of a second flagellar apparatus, the use of abietic acid as a growth substrate, the complete denitrification pathway, the possible production of a macrolide antibiotic and the assembly of multiple protein secretion systems.

Molecular evolution of LysR-type transcriptional regulation in Pseudomonas aeruginosa.

Signal perception and transduction through tightly coordinated circuits is integral to the survival and persistence of microbes in diverse ecological niches. The capacity to adapt to changes in the environment is central to their ability to thrive under adverse circumstances. Signal dependent transcriptional regulators are a key mechanism through which microbes assimilate environmental cues and mediate the appropriate adaptive response. By far the largest class of transcriptional regulator is the LysR-class, which is universally distributed among bacteria, archaea, and even eukaryotic organisms. The number of LysR-Type Transcriptional Regulators (LTTRs) varies among species with one of the largest repertoires encoded in the genome of the nosocomial pathogen Pseudomonas aeruginosa. To understand the evolutionary basis for this, we undertook to analyse the relationship between the LTTRs, both at the species and genus level. Phylogenetic analysis of the complete Pseudomonas LTTR dataset revealed significant cluster patterns based on full length and domain analysis. Interestingly, evidence of acquisition through horizontal gene transfer was rare, with divergent evolution apparently favoured. Furthermore, genes that appear to have been acquired, as well as those with a non-classical topological arrangement were clustered in distinct groups in the phylogenetic trees, indicating some ancestral association. The conservation within clusters identified in this study will provide a useful platform for future molecular analyses.

Subtilomycin: a new lantibiotic from Bacillus subtilis strain MMA7 isolated from the marine sponge Haliclona simulans.

Bacteriocins are attracting increased attention as an alternative to classic antibiotics in the fight against infectious disease and multidrug resistant pathogens. Bacillus subtilis strain MMA7 isolated from the marine sponge Haliclona simulans displays a broad spectrum antimicrobial activity, which includes Gram-positive and Gram-negative pathogens, as well as several pathogenic Candida species. This activity is in part associated with a newly identified lantibiotic, herein named as subtilomycin. The proposed biosynthetic cluster is composed of six genes, including protein-coding genes for LanB-like dehydratase and LanC-like cyclase modification enzymes, characteristic of the class I lantibiotics. The subtilomycin biosynthetic cluster in B. subtilis strain MMA7 is found in place of the sporulation killing factor (skf) operon, reported in many B. subtilis isolates and involved in a bacterial cannibalistic behaviour intended to delay sporulation. The presence of the subtilomycin biosynthetic cluster appears to be widespread amongst B. subtilis strains isolated from different shallow and deep water marine sponges. Subtilomycin possesses several desirable industrial and pharmaceutical physicochemical properties, including activity over a wide pH range, thermal resistance and water solubility. Additionally, the production of the lantibiotic subtilomycin could be a desirable property should B. subtilis strain MMA7 be employed as a probiotic in aquaculture applications.

Genome sequence of the biocontrol strain Pseudomonas fluorescens F113.

Pseudomonas fluorescens F113 is a plant growth-promoting rhizobacterium (PGPR) that has biocontrol activity against fungal plant pathogens and is a model for rhizosphere colonization. Here, we present its complete genome sequence, which shows that besides a core genome very similar to those of other strains sequenced within this species, F113 possesses a wide array of genes encoding specialized functions for thriving in the rhizosphere and interacting with eukaryotic organisms.

Single Seed Microbiota: Assembly and Transmission from Parent Plant to Seedling.

The seed acts as the primary inoculum source for the plant microbiota. Understanding the processes involved in its assembly and dynamics during germination and seedling emergence has the potential to allow for the improvement of crop establishment. Changes in the bacterial community structure were tracked in 1,000 individual seeds that were collected throughout seed developments of beans and radishes. Seeds were associated with a dominant bacterial taxon that represented more than 75% of all reads. The identity of this taxon was highly variable between the plants and within the seeds of the same plant. We identified selection as the main ecological process governing the succession of dominant taxa during seed filling and maturation. In a second step, we evaluated the seedling transmission of seed-borne taxa in 160 individual plants. While the initial bacterial abundance on seeds was not a good predictor of seedling transmission, the identities of the seed-borne taxa modified the phenotypes of seedlings. Overall, this work revealed that individual seeds are colonized by a few bacterial taxa of highly variable identity, which appears to be important for the early stages of plant development. IMPORTANCE Seeds are key components of plant fitness and are central to the sustainability of the agri-food system. Both the seed quality for food consumption and the seed vigor in agricultural settings can be influenced by the seed microbiota. Understanding the ecological processes involved in seed microbiota assembly will inform future practices for promoting the presence of important seed microorganisms for plant health and productivity. Our results highlighted that seeds were associated with one dominant bacterial taxon of variable taxonomic identity. This variety of dominant taxa was due to (i) spatial heterogeneity between and within plants and (ii) primary succession during seed development. According to neutral models, selection was the main driver of microbial community assembly for both plant species.

Defining Composition and Function of the Rhizosphere Microbiota of Barley Genotypes Exposed to Growth-Limiting Nitrogen Supplies.

The microbiota populating the rhizosphere, the interface between roots and soil, can modulate plant growth, development, and health. These microbial communities are not stochastically assembled from the surrounding soil, but their composition and putative function are controlled, at least partially, by the host plant. Here, we use the staple cereal barley as a model to gain novel insights into the impact of differential applications of nitrogen, a rate-limiting step for global crop production, on the host genetic control of the rhizosphere microbiota. Using a high-throughput amplicon sequencing survey, we determined that nitrogen availability for plant uptake is a factor promoting the selective enrichment of individual taxa in the rhizosphere of wild and domesticated barley genotypes. Shotgun sequencing and metagenome-assembled genomes revealed that this taxonomic diversification is mirrored by a functional specialization, manifested by the differential enrichment of multiple Gene Ontology terms, of the microbiota of plants exposed to nitrogen conditions limiting barley growth. Finally, a plant soil feedback experiment revealed that host control of the barley microbiota underpins the assembly of a phylogenetically diverse group of bacteria putatively required to sustain plant performance under nitrogen-limiting supplies. Taken together, our observations indicate that under nitrogen conditions limiting plant growth, host-microbe and microbe-microbe interactions fine-tune the host genetic selection of the barley microbiota at both taxonomic and functional levels. The disruption of these recruitment cues negatively impacts plant growth. IMPORTANCE The microbiota inhabiting the rhizosphere, the thin layer of soil surrounding plant roots, can promote the growth, development, and health of their host plants. Previous research indicated that differences in the genetic composition of the host plant coincide with variations in the composition of the rhizosphere microbiota. This is particularly evident when looking at the microbiota associated with input-demanding modern cultivated varieties and their wild relatives, which have evolved under marginal conditions. However, the functional significance of these differences remains to be fully elucidated. We investigated the rhizosphere microbiota of wild and cultivated genotypes of the global crop barley and determined that nutrient conditions limiting plant growth amplify the host control on microbes at the root-soil interface. This is reflected in a plant- and genotype-dependent functional specialization of the rhizosphere microbiota, which appears to be required for optimal plant growth. These findings provide novel insights into the significance of the rhizosphere microbiota for plant growth and sustainable agriculture.

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

Bacteria encode multiple protein secretion systems that are crucial for interaction with the environment and with hosts. In recent years, attention has focused on type VI secretion systems (T6SSs), which are specialized transporters widely encoded in Proteobacteria. The myriad of processes associated with these secretion systems could be explained by subclasses of T6SS, each involved in specialized functions. To assess diversity and predict function associated with different T6SSs, comparative genomic analysis of 34 Pseudomonas genomes was performed. This identified 70 T6SSs, with at least one locus in every strain, except for Pseudomonas stutzeri A1501. By comparing 11 core genes of the T6SS, it was possible to identify five main Pseudomonas phylogenetic clusters, with strains typically carrying T6SSs from more than one clade. In addition, most strains encode additional vgrG and hcp genes, which encode extracellular structural components of the secretion apparatus. Using a combination of phylogenetic and meta-analysis of transcriptome datasets it was possible to associate specific subsets of VgrG and Hcp proteins with each Pseudomonas T6SS clade. Moreover, a closer examination of the genomic context of vgrG genes in multiple strains highlights a number of additional genes associated with these regions. It is proposed that these genes may play a role in secretion or alternatively could be new T6S effectors.