Bacterial cooperative associations and dynamics in the biofilm microenvironments are of special interest in recent years. Knowledge of localized gene-expression and corresponding bacterial behaviors within the biofilm architecture at a global scale has been limited, due to a lack of robust technology to study limited number of cells in stratified layers of biofilms. With our recent pioneering developments in single bacterial cell transcriptomic analysis technology, we generated herein an unprecedented spatial transcriptome map of the mature in vitro Pseudomonas aeruginosa biofilm model, revealing contemporaneous yet altered bacterial behaviors at different layers within the biofilm architecture (i.e., surface, middle, and interior of the biofilm). Many genes encoding unknown functions were highly expressed at the biofilm-solid interphase, exposing a critical gap in the knowledge of their activities that maybe unique to this interior niche. Several genes of unknown functions are critical for biofilm formation. The in vivo importance of these unknown proteins was validated in invertebrate (fruit fly) and vertebrate (mouse) models. We envisage the future value of this report to the community, in aiding the further pathophysiological understanding of P. aeruginosa biofilms. Our approach will open doors to the study of bacterial functional genomics of different species in numerous settings.

Pseudomonas aeruginosa can grow to very high-cell-density (HCD) during infection of the cystic fibrosis (CF) lung. Phosphatidylcholine (PC), the major component of lung surfactant, has been hypothesized to support HCD growth of P. aeruginosa in vivo. The phosphorylcholine headgroup, a glycerol molecule, and two long-chain fatty acids (FAs) are released by enzymatic cleavage of PC by bacterial phospholipase C and lipases. Three different bacterial pathways, the choline, glycerol, and fatty acid degradation pathways, are then involved in the degradation of these PC components. Here, we identified five potential FA degradation (Fad) related fadBA-operons (fadBA1-5, each encoding 3-hydroxyacyl-CoA dehydrogenase and acyl-CoA thiolase). Through mutagenesis and growth analyses, we showed that three (fadBA145) of the five fadBA-operons are dominant in medium-chain and long-chain Fad. The triple fadBA145 mutant also showed reduced ability to degrade PC in vitro. We have previously shown that by partially blocking Fad, via mutagenesis of fadBA5 and fadDs, we could significantly reduce the ability of P. aeruginosa to replicate on FA and PC in vitro, as well as in the mouse lung. However, no studies have assessed the ability of mutants, defective in choline and/or glycerol degradation in conjunction with Fad, to grow on PC or in vivo. Hence, we constructed additional mutants (ΔfadBA145ΔglpD, ΔfadBA145ΔbetAB, and ΔfadBA145ΔbetABΔglpD) significantly defective in the ability to degrade FA, choline, and glycerol and, therefore, PC. The analysis of these mutants in the BALB/c mouse lung infection model showed significant inability to utilize PC in vitro, resulted in decreased replication fitness and competitiveness in vivo compared to the complement strain, although there was little to no variation in typical virulence factor production (e.g., hemolysin, lipase, and protease levels). This further supports the hypothesis that lung surfactant PC serves as an important nutrient for P. aeruginosa during CF lung infection.

The fatty acid (FA) degradation pathway of Pseudomonas aeruginosa, an opportunistic pathogen, was recently shown to be involved in nutrient acquisition during BALB/c mouse lung infection model. The source of FA in the lung is believed to be phosphatidylcholine, the major component of lung surfactant. Previous research indicated that P. aeruginosa has more than two fatty acyl-CoA synthetase genes (fadD; PA3299 and PA3300), which are responsible for activation of FAs using ATP and coenzyme A. Through a bioinformatics approach, 11 candidate genes were identified by their homology to the Escherichia coli FadD in the present study. Four new homologues of fadD (PA1617, PA2893, PA3860, and PA3924) were functionally confirmed by their ability to complement the E. coli fadD mutant on FA-containing media. Growth phenotypes of 17 combinatorial fadD mutants on different FAs, as sole carbon sources, indicated that the four new fadD homologues are involved in FA degradation, bringing the total number of P. aeruginosa fadD genes to six. Of the four new homologues, fadD4 (PA1617) contributed the most to the degradation of different chain length FAs. Growth patterns of various fadD mutants on plant-based perfumery substances, citronellic and geranic acids, as sole carbon and energy sources indicated that fadD4 is also involved in the degradation of these plant-derived compounds. A decrease in fitness of the sextuple fadD mutant, relative to the ΔfadD1D2 mutant, was only observed during BALB/c mouse lung infection at 24 h.

Prokaryotic cell transcriptomics has been limited to mixed or sub-population dynamics and individual cells within heterogeneous populations, which has hampered further understanding of spatiotemporal and stage-specific processes of prokaryotic cells within complex environments. Here we develop a ‘TRANSITomic’ approach to profile transcriptomes of single Burkholderia pseudomallei cells as they transit through host cell infection at defined stages, yielding pathophysiological insights. We find that B. pseudomallei transits through host cells during infection in three observable stages: vacuole entry; cytoplasmic escape and replication; and membrane protrusion, promoting cell-to-cell spread. The B. pseudomallei ‘TRANSITome’ reveals dynamic gene-expression flux during transit in host cells and identifies genes that are required for pathogenesis. We find several hypothetical proteins and assign them to virulence mechanisms, including attachment, cytoskeletal modulation, and autophagy evasion. The B. pseudomallei ‘TRANSITome’ provides prokaryotic single-cell transcriptomics information enabling high-resolution understanding of host-pathogen interactions.

It is generally believed that the Pseudomonas aeruginosa biofilm matrix itself acts as a molecular sieve or sink that contributes to significant levels of drug resistance, but it is becoming more apparent that multidrug efflux pumps induced during biofilm growth significantly enhance resistance levels. We present here a novel transcriptional regulator, PA3898, which controls biofilm formation and multidrug efflux pumps in P. aeruginosa. A mutant of this regulator significantly reduced the ability of P. aeruginosa to produce biofilm in vitro and affected its in vivo fitness and pathogenesis in Drosophila melanogaster and BALB/c mouse lung infection models. Transcriptome analysis revealed that PA3898 modulates essential virulence genes/pathways, including multidrug efflux pumps and phenazine biosynthesis. Chromatin immunoprecipitation sequencing (ChIP-seq) identified its DNA binding sequences and confirmed that PA3898 directly interacts with promoter regions of four genes/operons, two of which are mexAB-oprM and phz2. Coimmunoprecipitation revealed a regulatory partner of PA3898 as PA2100, and both are required for binding to DNA in electrophoretic mobility shift assays. PA3898 and PA2100 were given the names MdrR1 and MdrR2, respectively, as novel repressors of the mexAB-oprM multidrug efflux operon and activators for another multidrug efflux pump, EmrAB. The interaction between MdrR1 and MdrR2 at the promoter regions of their regulons was further characterized via localized surface plasmon resonance and DNA footprinting. These regulators directly repress the mexAB-oprM operon, independent of its well-established MexR regulator. Mutants of mdrR1 and mdrR2 caused increased resistance to multiple antibiotics in P. aeruginosa, validating the significance of these newly discovered regulators.