Computationally efficient assembly of a Pseudomonas aeruginosa gene expression compendium

Over the past two decades, thousands of RNA sequencing (RNA-seq) gene expression profiles of Pseudomonas aeruginosa have been made publicly available via the National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA). In the work we present here, we draw on over 2,300 P. aeruginosa transcriptomes from hundreds of studies performed by over seventy-five different research groups. We first developed a pipeline, using the Salmon pseudo-aligner and two different P. aeruginosa reference genomes (strains PAO1 and PA14), that transformed raw sequence data into a uniformly processed data in the form of sample-wise normalized counts. In this workflow, P. aeruginosa RNA-seq data are filtered using technically and biologically driven criteria with characteristics tailored to bacterial gene expression and that account for the effects of alignment to different reference genomes. The filtered data are then normalized to enable cross experiment comparisons. Finally, annotations are programmatically collected for those samples with sufficient meta-data and expression-based metrics are used to further enhance strain assignment for each sample. Our processing and quality control methods provide a scalable framework for taking full advantage of the troves of biological information hibernating in the depths of microbial gene expression data. The re-analysis of these data in aggregate is a powerful approach for hypothesis generation and testing, and this approach can be applied to transcriptome datasets in other species.