Compendium-Wide Analysis of Pseudomonas aeruginosa Core and Accessory Genes Reveals Transcriptional Patterns across Strains PAO1 and PA14

Pseudomonas aeruginosa is a ubiquitous pathogen. There is much diversity among P. aeruginosa strains, including two divergent but well-studied strains, PAO1 and PA14. ABSTRACT Pseudomonas aeruginosa is an opportunistic pathogen that causes difficult-to-treat infections. Two well-studied divergent P. aeruginosa strain types, PAO1 and PA14, have significant genomic heterogeneity, including diverse accessory genes present in only some strains. Genome content comparisons find core genes that are conserved across both PAO1 and PA14 strains and accessory genes that are present in only a subset of PAO1 and PA14 strains. Here, we use recently assembled transcriptome compendia of publicly available P. aeruginosa RNA sequencing (RNA-seq) samples to create two smaller compendia consisting of only strain PAO1 or strain PA14 samples with each aligned to their cognate reference genome. We confirmed strain annotations and identified other samples for inclusion by assessing each sample’s median expression of PAO1-only or PA14-only accessory genes. We then compared the patterns of core gene expression in each strain. To do so, we developed a method by which we analyzed genes in terms of which genes showed similar expression patterns across strain types. We found that some core genes had consistent correlated expression patterns across both compendia, while others were less stable in an interstrain comparison. For each accessory gene, we also determined core genes with correlated expression patterns. We found that stable core genes had fewer coexpressed neighbors that were accessory genes. Overall, this approach for analyzing expression patterns across strain types can be extended to other groups of genes, like phage genes, or applied for analyzing patterns beyond groups of strains, such as samples with different traits, to reveal a deeper understanding of regulation. IMPORTANCE Pseudomonas aeruginosa is a ubiquitous pathogen. There is much diversity among P. aeruginosa strains, including two divergent but well-studied strains, PAO1 and PA14. Understanding how these different strain-level traits manifest is important for identifying targets that regulate different traits of interest. With the availability of thousands of PAO1 and PA14 samples, we created two strain-specific RNA-seq compendia where each one contains hundreds of samples from PAO1 or PA14 strains and used them to compare the expression patterns of core genes that are conserved in both strain types and to determine which core genes have expression patterns that are similar to those of accessory genes that are unique to one strain or the other using an approach that we developed. We found a subset of core genes with different transcriptional patterns across PAO1 and PA14 strains and identified those core genes with expression patterns similar to those of strain-specific accessory genes.

[1]  Anand V. Sastry,et al.  Machine learning from Pseudomonas aeruginosa transcriptomes identifies independently modulated sets of genes associated with known transcriptional regulators , 2022, Nucleic acids research.

[2]  Alexandra J. Lee,et al.  Computationally Efficient Assembly of Pseudomonas aeruginosa Gene Expression Compendia , 2022, bioRxiv.

[3]  Maxx H. Tessmer,et al.  Perspectives on the Pseudomonas aeruginosa Type III Secretion System Effector ExoU and Its Subversion of the Host Innate Immune Response to Infection , 2021, Toxins.

[4]  T. Vernet,et al.  Determination of the two-component systems regulatory network reveals core and accessory regulations across Pseudomonas aeruginosa lineages , 2021, bioRxiv.

[5]  C. Viti,et al.  Symbiotic and Nonsymbiotic Members of the Genus Ensifer (syn. Sinorhizobium) Are Separated into Two Clades Based on Comparative Genomics and High-Throughput Phenotyping , 2020, bioRxiv.

[6]  D. Hogan,et al.  Conditional antagonism in co-cultures of Pseudomonas aeruginosa and Candida albicans: An intersection of ethanol and phosphate signaling distilled from dual-seq transcriptomics , 2020, bioRxiv.

[7]  Kimberley A. Lewis,et al.  Social Cooperativity of Bacteria during Reversible Surface Attachment in Young Biofilms: a Quantitative Comparison of Pseudomonas aeruginosa PA14 and PAO1 , 2020, mBio.

[8]  M. Touchon,et al.  Phylogenetic background and habitat drive the genetic diversification of Escherichia coli , 2020, bioRxiv.

[9]  Pedro P. Vallejo Ramirez,et al.  Contextual Flexibility in Pseudomonas aeruginosa Central Carbon Metabolism during Growth in Single Carbon Sources , 2019, mBio.

[10]  A. Conesa,et al.  Differential Modulation of Quorum Sensing Signaling through QslA in Pseudomonas aeruginosa Strains PAO1 and PA14 , 2019, Journal of bacteriology.

[11]  Ying Xu,et al.  DOOR: a prokaryotic operon database for genome analyses and functional inference , 2019, Briefings Bioinform..

[12]  X. Deng,et al.  An integrated genomic regulatory network of virulence-related transcriptional factors in Pseudomonas aeruginosa , 2019, Nature Communications.

[13]  Tapio Salakoski,et al.  The CAFA challenge reports improved protein function prediction and new functional annotations for hundreds of genes through experimental screens , 2019, Genome Biology.

[14]  J. Jeukens,et al.  The Pseudomonas aeruginosa Pan-Genome Provides New Insights on Its Population Structure, Horizontal Gene Transfer, and Pathogenicity , 2018, Genome biology and evolution.

[15]  Klaus Hornischer,et al.  BACTOME—a reference database to explore the sequence- and gene expression-variation landscape of Pseudomonas aeruginosa clinical isolates , 2018, Nucleic Acids Res..

[16]  E. Lander,et al.  Defining the core essential genome of Pseudomonas aeruginosa , 2018, Proceedings of the National Academy of Sciences.

[17]  Edoardo Pasolli,et al.  Whole-genome epidemiology, characterisation, and phylogenetic reconstruction of Staphylococcus aureus strains in a paediatric hospital , 2018, Genome Medicine.

[18]  S. Lory,et al.  Insertion sequences drive the emergence of a highly adapted human pathogen , 2018, bioRxiv.

[19]  S. Rice,et al.  Association between possession of ExoU and antibiotic resistance in Pseudomonas aeruginosa , 2018, PloS one.

[20]  A. Bush,et al.  Metabolic Phenotyping and Strain Characterisation of Pseudomonas aeruginosa Isolates from Cystic Fibrosis Patients Using Rapid Evaporative Ionisation Mass Spectrometry , 2018, Scientific Reports.

[21]  Ernest Y Lee,et al.  Multigenerational memory and adaptive adhesion in early bacterial biofilm communities , 2018, Proceedings of the National Academy of Sciences.

[22]  H. Zeighami,et al.  Diversity of virulence genes in multidrug resistant Pseudomonas aeruginosa isolated from burn wound infections. , 2018, Microbial pathogenesis.

[23]  Jie Tan,et al.  Unsupervised Extraction of Stable Expression Signatures from Public Compendia with an Ensemble of Neural Networks. , 2017, Cell systems.

[24]  F. Chang,et al.  Outbreak of imipenem‐resistant Acinetobacter baumannii in different wards at a regional hospital related to untrained bedside caregivers , 2017, American journal of infection control.

[25]  E. Reboud,et al.  Pseudomonas aeruginosa renews its virulence factors. , 2016, Environmental microbiology reports.

[26]  L. Rahme,et al.  Evidence for Direct Control of Virulence and Defense Gene Circuits by the Pseudomonas aeruginosa Quorum Sensing Regulator, MvfR , 2016, Scientific Reports.

[27]  Z. Karami,et al.  Production and Characterization of a Nitrilase from Pseudomonas aeruginosa RZ44 and its Potential for Nitrile Biotransformation , 2016, Iranian journal of biotechnology.

[28]  T. van Opijnen,et al.  Strain Dependent Genetic Networks for Antibiotic-Sensitivity in a Bacterial Pathogen with a Large Pan-Genome , 2016, PLoS pathogens.

[29]  J. Rello,et al.  Pseudomonas aeruginosa ventilator-associated pneumonia management , 2016, Infection and drug resistance.

[30]  C. Greene,et al.  ADAGE-Based Integration of Publicly Available Pseudomonas aeruginosa Gene Expression Data with Denoising Autoencoders Illuminates Microbe-Host Interactions , 2016, mSystems.

[31]  Xinghua Lu,et al.  Learning a hierarchical representation of the yeast transcriptomic machinery using an autoencoder model , 2016, BMC Bioinformatics.

[32]  M. Spies,et al.  Inhibition of Pseudomonas aeruginosa ExsA DNA-Binding Activity by N-Hydroxybenzimidazoles , 2015, Antimicrobial Agents and Chemotherapy.

[33]  Geoffrey L. Winsor,et al.  Clinical utilization of genomics data produced by the international Pseudomonas aeruginosa consortium , 2015, Front. Microbiol..

[34]  Ilias Tagkopoulos,et al.  Microbial Forensics: Predicting Phenotypic Characteristics and Environmental Conditions from Large-Scale Gene Expression Profiles , 2015, PLoS Comput. Biol..

[35]  G. Raman,et al.  Clinical and economic consequences of hospital-acquired resistant and multidrug-resistant Pseudomonas aeruginosa infections: a systematic review and meta-analysis , 2014, Antimicrobial Resistance and Infection Control.

[36]  Jasmine Lee,et al.  The hierarchy quorum sensing network in Pseudomonas aeruginosa , 2014, Protein & Cell.

[37]  Egon A Ozer,et al.  Characterization of the core and accessory genomes of Pseudomonas aeruginosa using bioinformatic tools Spine and AGEnt , 2014, BMC Genomics.

[38]  T. Blackwell,et al.  Structure and function of OprD protein in Pseudomonas aeruginosa: from antibiotic resistance to novel therapies. , 2012, International journal of medical microbiology : IJMM.

[39]  A. Dash,et al.  Plasmodium vivax lineages: geographical distribution, tandem repeat polymorphism, and phylogenetic relationship , 2011, Malaria Journal.

[40]  Raymond Lo,et al.  Pseudomonas Genome Database: improved comparative analysis and population genomics capability for Pseudomonas genomes , 2010, Nucleic Acids Res..

[41]  A. Martínez-Antonio,et al.  The Regulatory Network of Pseudomonas aeruginosa , 2011, Microbial Informatics and Experimentation.

[42]  A. Filloux,et al.  Protein secretion systems in Pseudomonas aeruginosa: A wealth of pathogenic weapons. , 2010, International journal of medical microbiology : IJMM.

[43]  Egon A. Ozer,et al.  The Accessory Genome of Pseudomonas aeruginosa , 2010, Microbiology and Molecular Biology Reviews.

[44]  Zixin Deng,et al.  Pathogenicity Islands PAPI-1 and PAPI-2 Contribute Individually and Synergistically to the Virulence of Pseudomonas aeruginosa Strain PA14 , 2010, Infection and Immunity.

[45]  A. Mengoni,et al.  Comparative genomics of VirR regulons in Clostridium perfringens strains , 2010, BMC Microbiology.

[46]  M. Schuster,et al.  Global position analysis of the Pseudomonas aeruginosa quorum‐sensing transcription factor LasR , 2009, Molecular microbiology.

[47]  J. Engel,et al.  Role of Pseudomonas aeruginosa type III effectors in disease. , 2009, Current opinion in microbiology.

[48]  S. Abedon Phage evolution and ecology. , 2009, Advances in applied microbiology.

[49]  E. Denamur,et al.  Escherichia , 2020, Definitions.

[50]  Daniel J. Kvitek,et al.  Variations in Stress Sensitivity and Genomic Expression in Diverse S. cerevisiae Isolates , 2008, PLoS genetics.

[51]  S. Brady,et al.  Paerucumarin, a New Metabolite Produced by the pvc Gene Cluster from Pseudomonas aeruginosa , 2008, Journal of bacteriology.

[52]  Kai Li,et al.  Exploring the functional landscape of gene expression: directed search of large microarray compendia , 2007, Bioinform..

[53]  M. Whiteley,et al.  Nutritional Cues Control Pseudomonas aeruginosa Multicellular Behavior in Cystic Fibrosis Sputum , 2007, Journal of bacteriology.

[54]  D. Hogan,et al.  Examination of Pseudomonas aeruginosa lasI regulation and 3-oxo-C12-homoserine lactone production using a heterologous Escherichia coli system. , 2007, FEMS microbiology letters.

[55]  Sean R. Davis,et al.  GEOquery: a bridge between the Gene Expression Omnibus (GEO) and BioConductor , 2007, Bioinform..

[56]  Wiley Interscience,et al.  Predictors of mortality in adults with cystic fibrosis , 2007, Pediatric pulmonology.

[57]  David A. D'Argenio,et al.  Growth phenotypes of Pseudomonas aeruginosa lasR mutants adapted to the airways of cystic fibrosis patients , 2007, Molecular microbiology.

[58]  Delbert Dueck,et al.  Clustering by Passing Messages Between Data Points , 2007, Science.

[59]  O. Kuipers,et al.  Regulation of Gene Expression in Streptococcus pneumoniae by Response Regulator 09 Is Strain Dependent , 2006, Journal of bacteriology.

[60]  Li Li,et al.  Genomic analysis reveals that Pseudomonas aeruginosa virulence is combinatorial , 2006, Genome Biology.

[61]  D. Newman,et al.  The phenazine pyocyanin is a terminal signalling factor in the quorum sensing network of Pseudomonas aeruginosa , 2006, Molecular microbiology.

[62]  Matthew A. Hibbs,et al.  Finding function: evaluation methods for functional genomic data , 2006, BMC Genomics.

[63]  M. Whiteley,et al.  Microarray Analysis of the Osmotic Stress Response in Pseudomonas aeruginosa , 2006, Journal of bacteriology.

[64]  D. Tifrea,et al.  A chemosensory system that regulates biofilm formation through modulation of cyclic diguanylate levels. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[65]  S. Molin,et al.  Heterogeneity of Biofilms Formed by Nonmucoid Pseudomonas aeruginosa Isolates from Patients with Cystic Fibrosis , 2005, Journal of Clinical Microbiology.

[66]  Jonathan R Edwards,et al.  Overview of nosocomial infections caused by gram-negative bacilli. , 2005, Clinical infectious diseases : an official publication of the Infectious Diseases Society of America.

[67]  Pradeep K. Singh,et al.  Cystic Fibrosis Sputum Supports Growth and Cues Key Aspects of Pseudomonas aeruginosa Physiology , 2005, Journal of bacteriology.

[68]  D. Bilton,et al.  Phenotypic variability of Pseudomonas aeruginosa in sputa from patients with acute infective exacerbation of cystic fibrosis and its impact on the validity of antimicrobial susceptibility testing. , 2005, The Journal of antimicrobial chemotherapy.

[69]  Dean Cheng,et al.  Pseudomonas aeruginosa Genome Database and PseudoCAP: facilitating community-based, continually updated, genome annotation , 2004, Nucleic Acids Res..

[70]  A. Hauser,et al.  Relative Contributions of Pseudomonas aeruginosa ExoU, ExoS, and ExoT to Virulence in the Lung , 2004, Infection and Immunity.

[71]  P. Giffard,et al.  Genetic Diversity among Community Methicillin-Resistant Staphylococcus aureus Strains Causing Outpatient Infections in Australia , 2004, Journal of Clinical Microbiology.

[72]  Hongwei D. Yu,et al.  Cross-Sectional Analysis of Clinical and Environmental Isolates of Pseudomonas aeruginosa: Biofilm Formation, Virulence, and Genome Diversity , 2004, Infection and Immunity.

[73]  S. Worgall,et al.  Transcriptome Analysis of the Response of Pseudomonas aeruginosa to Hydrogen Peroxide , 2004, Journal of bacteriology.

[74]  R. Drijber,et al.  TNT Biotransformation and Detoxification by a Pseudomonas Aeruginosa Strain , 2003, Biodegradation.

[75]  J. Ramos Genomics, life style and molecular architecture , 2004 .

[76]  S. de Bentzmann,et al.  Genotypic and phenotypic analysis of type III secretion system in a cohort of Pseudomonas aeruginosa bacteremia isolates: evidence for a possible association between O serotypes and exo genes. , 2003, The Journal of infectious diseases.

[77]  J. Engel Molecular Pathogenesis of Acute Pseudomonas Aeruginosa Infections , 2003 .

[78]  M. Vasil,et al.  GeneChip® expression analysis of the iron starvation response in Pseudomonas aeruginosa: identification of novel pyoverdine biosynthesis genes , 2002, Molecular microbiology.

[79]  Frederick M. Ausubel,et al.  Identification of Virulence Genes in a Pathogenic Strain of Pseudomonas aeruginosa by Representational Difference Analysis , 2002, Journal of bacteriology.

[80]  Roger E. Bumgarner,et al.  Gene expression in Pseudomonas aeruginosa biofilms , 2001, Nature.

[81]  S. Lory,et al.  Complete genome sequence of Pseudomonas aeruginosa PAO1, an opportunistic pathogen , 2000, Nature.

[82]  F. Ausubel,et al.  Killing of Caenorhabditis elegans by Pseudomonas aeruginosa used to model mammalian bacterial pathogenesis. , 1999, Proceedings of the National Academy of Sciences of the United States of America.

[83]  Lei Zhu,et al.  ExoU expression by Pseudomonas aeruginosa correlates with acute cytotoxicity and epithelial injury , 1997, Molecular microbiology.

[84]  F. Ausubel,et al.  Common virulence factors for bacterial pathogenicity in plants and animals. , 1995, Science.

[85]  S. Lory,et al.  Products of three accessory genes, pilB, pilC, and pilD, are required for biogenesis of Pseudomonas aeruginosa pili , 1990, Journal of bacteriology.

[86]  D. Galas,et al.  Functional promoters created by the insertion of transposable element IS1. , 1986, Journal of molecular biology.

[87]  N. Glansdorff,et al.  Activation of gene expression by IS2 and IS3. , 1981, Cold Spring Harbor symposia on quantitative biology.