Computationally efficient assembly of a Pseudomonas aeruginosa gene expression compendium

Abstract

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. SignificancePseudomonas aeruginosa causes a wide range of infections including chronic infections associated with cystic fibrosis. P. aeruginosa infections are difficult to treat and people with CF-associated P. aeruginosa infections often have poor clinical outcomes. To aid the study of this important pathogen, we developed a methodology that facilitates analyses across experiments, strains, and conditions. We aligned, filtered for quality and normalized thousands of P. aeruginosa RNA-seq gene expression profiles that were publicly available via the National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA). The workflow that we present can be efficiently scaled to incorporate new data and applied to the analysis of other species.