High-throughput protein characterization by complementation using DNA barcoded fragment libraries.

Our ability to predict, control, or design biological function is fundamentally limited by poorly annotated gene function. This can be particularly challenging in non-model systems. Accordingly, there is motivation for new high-throughput methods for accurate functional annotation. Here, we used complementation of auxotrophs and DNA barcode sequencing (Coaux-Seq) to enable high-throughput characterization of protein function. Fragment libraries from eleven genetically diverse bacteria were tested in twenty different auxotrophic strains of Escherichia coli to identify genes that complement missing biochemical activity. We recovered 41% of expected hits, with effectiveness ranging per source genome, and observed success even with distant E. coli relatives like Bacillus subtilis and Bacteroides thetaiotaomicron. Coaux-Seq provided the first experimental validation for 53 proteins, of which 11 are less than 40% identical to an experimentally characterized protein. Among the unexpected function identified was a sulfate uptake transporter, an O-succinylhomoserine sulfhydrylase for methionine synthesis, and an aminotransferase. We also identified instances of cross-feeding wherein protein overexpression and nearby non-auxotrophic strains enabled growth. Altogether, Coaux-Seq's utility is demonstrated, with future applications in ecology, health, and engineering.

Fitness factors impacting survival of a subsurface bacterium in contaminated groundwater.

Many factors contribute to the ability of a microbial species to persist when encountering complexly contaminated environments including time of exposure, the nature and concentration of contaminants, availability of nutritional resources, and possession of a combination of appropriate molecular mechanisms needed for survival. Herein we sought to identify genes that are most important for survival of Gram-negative Enterobacteriaceae in contaminated groundwater environments containing high concentrations of nitrate and metals using the metal-tolerant Oak Ridge Reservation (ORR) isolate, Pantoea sp. MT58 (MT58). Survival fitness experiments in which a randomly barcoded transposon insertion (RB-TnSeq) library of MT58 was exposed directly to contaminated ORR groundwater samples from across a nitrate and mixed metal contamination plume were used to identify genes important for survival with increasing exposure times and concentrations of contaminants, and availability of a carbon source. Genes involved in controlling and using carbon, encoding transcriptional regulators, and related to Gram-negative outer membrane processes were among those found to be important for survival in contaminated ORR groundwater. A comparative genomics analysis of 75 Pantoea genus strains allowed us to further separate the survival determinants into core and non-core genes in the Pantoea pangenome, revealing insights into the survival of subsurface microorganisms during contaminant plume intrusion.

Barcoded overexpression screens in gut Bacteroidales identify genes with roles in carbon utilization and stress resistance.

A mechanistic understanding of host-microbe interactions in the gut microbiome is hindered by poorly annotated bacterial genomes. While functional genomics can generate large gene-to-phenotype datasets to accelerate functional discovery, their applications to study gut anaerobes have been limited. For instance, most gain-of-function screens of gut-derived genes have been performed in Escherichia coli and assayed in a small number of conditions. To address these challenges, we develop Barcoded Overexpression BActerial shotgun library sequencing (Boba-seq). We demonstrate the power of this approach by assaying genes from diverse gut Bacteroidales overexpressed in Bacteroides thetaiotaomicron. From hundreds of experiments, we identify new functions and phenotypes for 29 genes important for carbohydrate metabolism or tolerance to antibiotics or bile salts. Highlights include the discovery of a D-glucosamine kinase, a raffinose transporter, and several routes that increase tolerance to ceftriaxone and bile salts through lipid biosynthesis. This approach can be readily applied to develop screens in other strains and additional phenotypic assays.

Molecular mechanisms and environmental adaptations of flagellar loss and biofilm growth of Rhodanobacter under environmental stress.

Biofilms aid bacterial adhesion to surfaces via direct and indirect mechanisms, and formation of biofilms is considered as an important strategy for adaptation and survival in suboptimal environmental conditions. However, the molecular underpinnings of biofilm formation in subsurface sediment/groundwater ecosystems where microorganisms often experience fluctuations in nutrient input, pH, and nitrate or metal concentrations are underexplored. We examined biofilm formation under different nutrient, pH, metal, and nitrate regimens of 16 Rhodanobacter strains isolated from subsurface groundwater wells spanning diverse levels of pH (3.5 to 5) and nitrates (13.7 to 146 mM). Eight Rhodanobacter strains demonstrated significant biofilm growth under low pH, suggesting adaptations for survival and growth at low pH. Biofilms were intensified under aluminum stress, particularly in strains possessing fewer genetic traits associated with biofilm formation, findings warranting further investigation. Through random barcode transposon-site sequencing (RB-TnSeq), proteomics, use of specific mutants, and transmission electron microscopy analysis, we discovered flagellar loss under aluminum stress, indicating a potential relationship between motility, metal tolerance, and biofilm growth. Comparative genomic analyses revealed the absence of flagella and chemotaxis genes and the presence of a putative type VI secretion system in the highly biofilm-forming strain FW021-MT20. In this study we identified genetic determinants associated with biofilm growth under metal stress in a predominant environmental genus, Rhodanobacter, and identified traits aiding survival and adaptation to contaminated subsurface environments.

High-throughput genetics enables identification of nutrient utilization and accessory energy metabolism genes in a model methanogen.

Archaea are widespread in the environment and play fundamental roles in diverse ecosystems; however, characterization of their unique biology requires advanced tools. This is particularly challenging when characterizing gene function. Here, we generate randomly barcoded transposon libraries in the model methanogenic archaeon Methanococcus maripaludis and use high-throughput growth methods to conduct fitness assays (RB-TnSeq) across over 100 unique growth conditions. Using our approach, we identified new genes involved in nutrient utilization and response to oxidative stress. We identified novel genes for the usage of diverse nitrogen sources in M. maripaludis including a putative regulator of alanine deamination and molybdate transporters important for nitrogen fixation. Furthermore, leveraging the fitness data, we inferred that M. maripaludis can utilize additional nitrogen sources including ʟ-glutamine, ᴅ-glucuronamide, and adenosine. Under autotrophic growth conditions, we identified a gene encoding a domain of unknown function (DUF166) that is important for fitness and hypothesize that it has an accessory role in carbon dioxide assimilation. Finally, comparing fitness costs of oxygen versus sulfite stress, we identified a previously uncharacterized class of dissimilatory sulfite reductase-like proteins (Dsr-LP; group IIId) that is important during growth in the presence of sulfite. When overexpressed, Dsr-LP conferred sulfite resistance and enabled use of sulfite as the sole sulfur source. The high-throughput approach employed here allowed for generation of a large-scale data set that can be used as a resource to further understand gene function and metabolism in the archaeal domain.IMPORTANCEArchaea are widespread in the environment, yet basic aspects of their biology remain underexplored. To address this, we apply randomly barcoded transposon libraries (RB-TnSeq) to the model archaeon Methanococcus maripaludis. RB-TnSeq coupled with high-throughput growth assays across over 100 unique conditions identified roles for previously uncharacterized genes, including several encoding proteins with domains of unknown function (DUFs). We also expand on our understanding of carbon and nitrogen metabolism and characterize a group IIId dissimilatory sulfite reductase-like protein as a functional sulfite reductase. This data set encompasses a wide range of additional conditions including stress, nitrogen fixation, amino acid supplementation, and autotrophy, thus providing an extensive data set for the archaeal community to mine for characterizing additional genes of unknown function.

A phage tail-like bacteriocin suppresses competitors in metapopulations of pathogenic bacteria.

Bacteria can repurpose their own bacteriophage viruses (phage) to kill competing bacteria. Phage-derived elements are frequently strain specific in their killing activity, although there is limited evidence that this specificity drives bacterial population dynamics. Here, we identified intact phage and their derived elements in a metapopulation of wild plant-associated Pseudomonas genomes. We discovered that the most abundant viral cluster encodes a phage remnant resembling a phage tail called a tailocin, which bacteria have co-opted to kill bacterial competitors. Each pathogenic Pseudomonas strain carries one of a few distinct tailocin variants that target the variable polysaccharides in the outer membrane of co-occurring pathogenic Pseudomonas strains. Analysis of herbarium samples from the past 170 years revealed that the same tailocin and bacterial receptor variants have persisted in Pseudomonas populations. These results suggest that tailocin genetic diversity can be mined to develop targeted "tailocin cocktails" for microbial control.

Genomic and environmental controls on Castellaniella biogeography in an anthropogenically disturbed subsurface.

Castellaniella species have been isolated from a variety of mixed-waste environments including the nitrate and multiple metal-contaminated subsurface at the Oak Ridge Reservation (ORR). Previous studies examining microbial community composition and nitrate removal at ORR during biostimulation efforts reported increased abundances of members of the Castellaniella genus concurrent with increased denitrification rates. Thus, we asked how genomic and abiotic factors control the Castellaniella biogeography at the site to understand how these factors may influence nitrate transformation in an anthropogenically impacted setting. We report the isolation and characterization of several Castellaniella strains from the ORR subsurface. Five of these isolates match at 100% identity (at the 16S rRNA gene V4 region) to two Castellaniella amplicon sequence variants (ASVs), ASV1 and ASV2, that have persisted in the ORR subsurface for at least 2 decades. However, ASV2 has consistently higher relative abundance in samples taken from the site and was also the dominant blooming denitrifier population during a prior biostimulation effort. We found that the ASV2 representative strain has greater resistance to mixed metal stress than the ASV1 representative strains. We attribute this resistance, in part, to the large number of unique heavy metal resistance genes identified on a genomic island in the ASV2 representative genome. Additionally, we suggest that the relatively lower fitness of ASV1 may be connected to the loss of the nitrous oxide reductase (nos) operon (and associated nitrous oxide reductase activity) due to the insertion at this genomic locus of a mobile genetic element carrying copper resistance genes. This study demonstrates the value of integrating genomic, environmental, and phenotypic data to characterize the biogeography of key microorganisms in contaminated sites.

An expanded transcriptome atlas for Bacteroides thetaiotaomicron reveals a small RNA that modulates tetracycline sensitivity.

Plasticity in gene expression allows bacteria to adapt to diverse environments. This is particularly relevant in the dynamic niche of the human intestinal tract; however, transcriptional networks remain largely unknown for gut-resident bacteria. Here we apply differential RNA sequencing (RNA-seq) and conventional RNA-seq to the model gut bacterium Bacteroides thetaiotaomicron to map transcriptional units and profile their expression levels across 15 in vivo-relevant growth conditions. We infer stress- and carbon source-specific transcriptional regulons and expand the annotation of small RNAs (sRNAs). Integrating this expression atlas with published transposon mutant fitness data, we predict conditionally important sRNAs. These include MasB, which downregulates tetracycline tolerance. Using MS2 affinity purification and RNA-seq, we identify a putative MasB target and assess its role in the context of the MasB-associated phenotype. These data-publicly available through the Theta-Base web browser ( http://micromix.helmholtz-hiri.de/bacteroides/ )-constitute a valuable resource for the microbiome community.

A weaponized phage suppresses competitors in historical and modern metapopulations of pathogenic bacteria.

Bacteriophages, the viruses of bacteria, are proposed to drive bacterial population dynamics, yet direct evidence of their impact on natural populations is limited. Here we identified viral sequences in a metapopulation of wild plant-associated Pseudomonas spp. genomes. We discovered that the most abundant viral cluster does not encode an intact phage but instead encodes a tailocin - a phage-derived element that bacteria use to kill competitors for interbacterial warfare. Each pathogenic Pseudomonas sp. strain carries one of a few distinct tailocin variants, which target variable polysaccharides in the outer membrane of co-occurring pathogenic strains. Analysis of historic herbarium samples from the last 170 years revealed that the same tailocin and receptor variants have persisted in the Pseudomonas populations for at least two centuries, suggesting the continued use of a defined set of tailocin haplotypes and receptors. These results indicate that tailocin genetic diversity can be mined to develop targeted "tailocin cocktails" for microbial control. One-Sentence Summary: Bacterial pathogens in a host-associated metapopulation use a repurposed prophage to kill their competitors.

Randomly barcoded transposon mutant libraries for gut commensals I: Strategies for efficient library construction.

Randomly barcoded transposon mutant libraries are powerful tools for studying gene function and organization, assessing gene essentiality and pathways, discovering potential therapeutic targets, and understanding the physiology of gut bacteria and their interactions with the host. However, construction of high-quality libraries with uniform representation can be challenging. In this review, we survey various strategies for barcoded library construction, including transposition systems, methods of transposon delivery, optimal library size, and transconjugant selection schemes. We discuss the advantages and limitations of each approach, as well as factors to consider when selecting a strategy. In addition, we highlight experimental and computational advances in arraying condensed libraries from mutant pools. We focus on examples of successful library construction in gut bacteria and their application to gene function studies and drug discovery. Given the need for understanding gene function and organization in gut bacteria, we provide a comprehensive guide for researchers to construct randomly barcoded transposon mutant libraries.

Randomly barcoded transposon mutant libraries for gut commensals II: Applying libraries for functional genetics.

The critical role of the intestinal microbiota in human health and disease is well recognized. Nevertheless, there are still large gaps in our understanding of the functions and mechanisms encoded in the genomes of most members of the gut microbiota. Genome-scale libraries of transposon mutants are a powerful tool to help us address this gap. Recent advances in barcoded transposon mutagenesis have dramatically lowered the cost of mutant fitness determination in hundreds of in vitro and in vivo experimental conditions. In an accompanying review, we discuss recent advances and caveats for the construction of pooled and arrayed barcoded transposon mutant libraries in human gut commensals. In this review, we discuss how these libraries can be used across a wide range of applications, the technical aspects involved, and expectations for such screens.

A chemically-defined growth medium to support Lactobacillus-Acetobacter sp. community analysis.

Lactobacilli and Acetobacter sp. are commercially important bacteria that often form communities in natural fermentations, including food preparations, spoilage, and in the digestive tract of the fruit fly Drosophila melanogaster. Communities of these bacteria are widespread and prolific, despite numerous strain-specific auxotrophies, suggesting they have evolved nutrient interdependencies that regulate their growth. The use of a chemically-defined medium (CDM) supporting the growth of both groups of bacteria would facilitate the identification of the molecular mechanisms for the metabolic interactions between them. While numerous CDMs have been developed that support specific strains of lactobacilli or Acetobacter, there has not been a medium formulated to support both genera. We developed such a medium, based on a previous CDM designed for growth of lactobacilli, by modifying the nutrient abundances to improve growth yield. We further simplified the medium by substituting casamino acids in place of individual amino acids and the standard Wolfe's vitamins and mineral stocks in place of individual vitamins and minerals, resulting in a reduction from 40 to 8 stock solutions. These stock solutions can be used to prepare several CDM formulations that support robust growth of numerous lactobacilli and Acetobacters. Here, we provide the composition and several examples of its use, which is important for tractability in dissecting the genetic and metabolic basis of natural bacterial species interactions.

A mutant fitness compendium in Bifidobacteria reveals molecular determinants of colonization and host-microbe interactions.

Bifidobacteria commonly represent a dominant constituent of human gut microbiomes during infancy, influencing nutrition, immune development, and resistance to infection. Despite interest as a probiotic therapy, predicting the nutritional requirements and health-promoting effects of Bifidobacteria is challenging due to major knowledge gaps. To overcome these deficiencies, we used large-scale genetics to create a compendium of mutant fitness in Bifidobacterium breve (Bb). We generated a high density, randomly barcoded transposon insertion pool in Bb, and used this pool to determine Bb fitness requirements during colonization of germ-free mice and chickens with multiple diets and in response to hundreds of in vitro perturbations. To enable mechanistic investigation, we constructed an ordered collection of insertion strains covering 1462 genes. We leveraged these tools to improve models of metabolic pathways, reveal unexpected host- and diet-specific requirements for colonization, and connect the production of immunomodulatory molecules to growth benefits. These resources will greatly reduce the barrier to future investigations of this important beneficial microbe.

Simvastatin induces human gut bacterial cell surface genes.

Drugs intended to target mammalian cells can have broad off-target effects on the human gut microbiota with potential downstream consequences for drug efficacy and side effect profiles. Yet, despite a rich literature on antibiotic resistance, we still know very little about the mechanisms through which commensal bacteria evade non-antibiotic drugs. Here, we focus on statins, one of the most prescribed drug types in the world and an essential tool in the prevention and treatment of high circulating cholesterol levels. Prior work in humans, mice, and cell culture support an off-target effect of statins on human gut bacteria; however, the genetic determinants of statin sensitivity remain unknown. We confirmed that simvastatin inhibits the growth of diverse human gut bacterial strains grown in communities and in pure cultures. Drug sensitivity varied between phyla and was dose-dependent. We selected two representative simvastatin-sensitive species for more in-depth analysis: Eggerthella lenta (phylum: Actinobacteriota) and Bacteroides thetaiotaomicron (phylum: Bacteroidota). Transcriptomics revealed that both bacterial species upregulate genes in response to simvastatin that alter the cell membrane, including fatty acid biogenesis (E. lenta) and drug efflux systems (B. thetaiotaomicron). Transposon mutagenesis identified a key efflux system in B. thetaiotaomicron that enables growth in the presence of statins. Taken together, these results emphasize the importance of the bacterial cell membrane in countering the off-target effects of host-targeted drugs. Continued mechanistic dissection of the various mechanisms through which the human gut microbiota evades drugs will be essential to understand and predict the effects of drug administration in human cohorts and the potential downstream consequences for health and disease.

Genome-wide fitness profiling reveals molecular mechanisms that bacteria use to interact with Trichoderma atroviride exometabolites.

Trichoderma spp. are ubiquitous rhizosphere fungi capable of producing several classes of secondary metabolites that can modify the dynamics of the plant-associated microbiome. However, the bacterial-fungal mechanisms that mediate these interactions have not been fully characterized. Here, a random barcode transposon-site sequencing (RB-TnSeq) approach was employed to identify bacterial genes important for fitness in the presence of Trichoderma atroviride exudates. We selected three rhizosphere bacteria with RB-TnSeq mutant libraries that can promote plant growth: the nitrogen fixers Klebsiella michiganensis M5aI and Herbaspirillum seropedicae SmR1, and Pseudomonas simiae WCS417. As a non-rhizosphere species, Pseudomonas putida KT2440 was also included. From the RB-TnSeq data, nitrogen-fixing bacteria competed mainly for iron and required the siderophore transport system TonB/ExbB for optimal fitness in the presence of T. atroviride exudates. In contrast, P. simiae and P. putida were highly dependent on mechanisms associated with membrane lipid modification that are required for resistance to cationic antimicrobial peptides (CAMPs). A mutant in the Hog1-MAP kinase (Δtmk3) gene of T. atroviride showed altered expression patterns of many nonribosomal peptide synthetase (NRPS) biosynthetic gene clusters with potential antibiotic activity. In contrast to exudates from wild-type T. atroviride, bacterial mutants containing lesions in genes associated with resistance to antibiotics did not show fitness defects when RB-TnSeq libraries were exposed to exudates from the Δtmk3 mutant. Unexpectedly, exudates from wild-type T. atroviride and the Δtmk3 mutant rescued purine auxotrophic mutants of H. seropedicae, K. michiganensis and P. simiae. Metabolomic analysis on exudates from wild-type T. atroviride and the Δtmk3 mutant showed that both strains excrete purines and complex metabolites; functional Tmk3 is required to produce some of these metabolites. This study highlights the complex interplay between Trichoderma-metabolites and soil bacteria, revealing both beneficial and antagonistic effects, and underscoring the intricate and multifaceted nature of this relationship.

Geochemical constraints on bacteriophage infectivity in terrestrial environments.

Lytic phages can be potent and selective inhibitors of microbial growth and can have profound impacts on microbiome composition and function. However, there is uncertainty about the biogeochemical conditions under which phage predation modulates microbial ecosystem function, particularly in terrestrial systems. Ionic strength is critical for infection of bacteria by many phages, but quantitative data is limited on the ion thresholds for phage infection that can be compared with environmental ion concentrations. Similarly, while carbon composition varies in the environment, we do not know how this variability influences the impact of phage predation on microbiome function. Here, we measured the half-maximal effective concentrations (EC50) of 80 different inorganic ions for the infection of E. coli with two canonical dsDNA and ssRNA phages, T4 and MS2, respectively. Many alkaline earth metals and alkali metals enabled lytic infection but the ionic strength thresholds varied for different ions between phages. Additionally, using a freshwater nitrate-reducing microbiome, we found that the ability of lytic phages to influence nitrate reduction end-products depended upon the carbon source as well as ionic strength. For all phage:host pairs, the ion EC50s for phage infection exceeded the ion concentrations found in many terrestrial freshwater systems. Thus, our findings support a model where phages most influence terrestrial microbial functional ecology in hot spots and hot moments such as metazoan guts, drought influenced soils, or biofilms where ion concentration is locally or transiently elevated and nutrients are available to support the growth of specific phage hosts.

Large-scale genetic characterization of the model sulfate-reducing bacterium, Desulfovibrio vulgaris Hildenborough.

Sulfate-reducing bacteria (SRB) are obligate anaerobes that can couple their growth to the reduction of sulfate. Despite the importance of SRB to global nutrient cycles and their damage to the petroleum industry, our molecular understanding of their physiology remains limited. To systematically provide new insights into SRB biology, we generated a randomly barcoded transposon mutant library in the model SRB Desulfovibrio vulgaris Hildenborough (DvH) and used this genome-wide resource to assay the importance of its genes under a range of metabolic and stress conditions. In addition to defining the essential gene set of DvH, we identified a conditional phenotype for 1,137 non-essential genes. Through examination of these conditional phenotypes, we were able to make a number of novel insights into our molecular understanding of DvH, including how this bacterium synthesizes vitamins. For example, we identified DVU0867 as an atypical L-aspartate decarboxylase required for the synthesis of pantothenic acid, provided the first experimental evidence that biotin synthesis in DvH occurs via a specialized acyl carrier protein and without methyl esters, and demonstrated that the uncharacterized dehydrogenase DVU0826:DVU0827 is necessary for the synthesis of pyridoxal phosphate. In addition, we used the mutant fitness data to identify genes involved in the assimilation of diverse nitrogen sources and gained insights into the mechanism of inhibition of chlorate and molybdate. Our large-scale fitness dataset and RB-TnSeq mutant library are community-wide resources that can be used to generate further testable hypotheses into the gene functions of this environmentally and industrially important group of bacteria.

An integrated transcriptomics-functional genomics approach reveals a small RNA that modulates Bacteroides thetaiotaomicron sensitivity to tetracyclines.

Gene expression plasticity allows bacteria to adapt to diverse environments, tie their metabolism to available nutrients, and cope with stress. This is particularly relevant in a niche as dynamic and hostile as the human intestinal tract, yet transcriptional networks remain largely unknown in gut Bacteroides spp. Here, we map transcriptional units and profile their expression levels in Bacteroides thetaiotaomicron over a suite of 15 defined experimental conditions that are relevant in vivo , such as variation of temperature, pH, and oxygen tension, exposure to antibiotic stress, and growth on simple carbohydrates or on host mucin-derived glycans. Thereby, we infer stress- and carbon source-specific transcriptional regulons, including conditional expression of capsular polysaccharides and polysaccharide utilization loci, and expand the annotation of small regulatory RNAs (sRNAs) in this organism. Integrating this comprehensive expression atlas with transposon mutant fitness data, we identify conditionally important sRNAs. One example is MasB, whose inactivation led to increased bacterial tolerance of tetracyclines. Using MS2 affinity purification coupled with RNA sequencing, we predict targets of this sRNA and discuss their potential role in the context of the MasB-associated phenotype. Together, this transcriptomic compendium in combination with functional sRNA genomics-publicly available through a new iteration of the 'Theta-Base' web browser (www.helmholtz-hiri.de/en/datasets/bacteroides-v2)-constitutes a valuable resource for the microbiome and sRNA research communities alike.

Multicopy suppressor screens reveal convergent evolution of single-gene lysis proteins.

Single-strand RNA (ssRNA) Fiersviridae phages cause host lysis with a product of single gene (sgl for single-gene lysis; product Sgl) that induces autolysis. Many different Sgls have been discovered, but the molecular targets of only a few have been identified. In this study, we used a high-throughput genetic screen to uncover genome-wide host suppressors of diverse Sgls. In addition to validating known molecular mechanisms, we discovered that the Sgl of PP7, an ssRNA phage of Pseudomonas aeruginosa, targets MurJ, the flippase responsible for lipid II export, previously shown to be the target of the Sgl of coliphage M. These two Sgls, which are unrelated and predicted to have opposite membrane topology, thus represent a case of convergent evolution. We extended the genetic screens to other uncharacterized Sgls and uncovered a common set of multicopy suppressors, suggesting that these Sgls act by the same or similar mechanism.

Genome-wide identification of fitness determinants in the Xanthomonas campestris bacterial pathogen during early stages of plant infection.

Plant diseases are an important threat to food production. While major pathogenicity determinants required for disease have been extensively studied, less is known on how pathogens thrive during host colonization, especially at early infection stages. Here, we used randomly barcoded-transposon insertion site sequencing (RB-TnSeq) to perform a genome-wide screen and identify key bacterial fitness determinants of the vascular pathogen Xanthomonas campestris pv campestris (Xcc) during infection of the cauliflower host plant (Brassica oleracea). This high-throughput analysis was conducted in hydathodes, the natural entry site of Xcc, in xylem sap and in synthetic media. Xcc did not face a strong bottleneck during hydathode infection. In total, 181 genes important for fitness were identified in plant-associated environments with functional enrichment in genes involved in metabolism but only few genes previously known to be involved in virulence. The biological relevance of 12 genes was independently confirmed by phenotyping single mutants. Notably, we show that XC_3388, a protein with no known function (DUF1631), plays a key role in the adaptation and virulence of Xcc possibly through c-di-GMP-mediated regulation. This study revealed yet unsuspected social behaviors adopted by Xcc individuals when confined inside hydathodes at early infection stages.