Total transcript amplification (TTA) from single eukaryotic cells for transcriptome analysis is established, but TTA from a single prokaryotic cell presents additional challenges with much less starting material, the lack of poly(A)-tails, and the fact that the messages can be polycistronic. Here, we describe a novel method for single-bacterium TTA using a model organism, Burkholderia thailandensis, exposed to a subinhibitory concentration of the antibacterial agent, glyphosate. Utilizing a B. thailandensis microarray to assess the TTA method showed low fold-change bias (less than twofold difference and Pearson correlation coefficient R % 0.87-0.89) and drop-outs (4%-6% of 2842 detectable genes), compared with data obtained from the larger-scale nonamplified RNA samples. Further analysis of the microarray data suggests that B. thailandensis, when exposed to the aromatic amino acid biosynthesis inhibitor glyphosate, induces (or represses) genes to possibly recuperate and balance the intracellular amino acid pool. We validated our single-cell microarray data at the multi-cell and single-cell levels with lacZ and gfp reporter-gene fusions, respectively. Sanger sequencing of 192 clones generated from the TTA product of a single cell, with and without enrichment by elimination of rRNA and tRNA, detected only B. thailandensis sequences with no contamination. These data indicate that RNA-seq of TTA from a single cell is possible using this novel method.

A close interconnection between nutrient metabolism and virulence factor expression contributes to the pathophysiology of Pseudomonas aeruginosa as a successful pathogen. P. aeruginosa fatty acid (FA) degradation is complicated with multiple acyl-CoA synthetase homologs (FadDs) expressed in vivo in lung tissue during cystic fibrosis infections. The promoters of two genetically linked P. aeruginosa fadD genes (fadD1 and fadD2) were mapped and northern blot analysis indicated they could exist on two different transcripts. These FadDs contain ATP/AMP signature and FA-binding motifs highly homologous to those of the Escherichia coli FadD. Upon introduction into an E. coli fadD -/fadR - double mutant, both P. aeruginosa fadDs functionally complemented the E. coli fadD -/fadR - mutant, allowing degradation of different chain-length FAs. Chromosomal mutagenesis, growth analysis, induction studies, and determination of kinetic parameters suggested that FadD1 has a substrate preference for long-chain FAs while FadD2 prefers shorter-chain FAs. When compared to the wild type strain, the fadD2 mutant exhibited decreased production of lipase, protease, rhamnolipid and phospholipase, and retardation of both swimming and swarming motilities. Interestingly, fadD1 mutant showed only increased swarming motility. Growth analysis of the fadD mutants showed noticeable deficiencies in utilizing FAs and phosphatidylcholine (major components of lung surfactant) as the sole carbon source. This defect translated into decreased in vivo fitness of P. aeruginosa in a BALB/c mouse lung infection model, supporting the role of lipids as a significant nutrient source for this bacterium in vivo.

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.

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.