Mapping genetic and phenotypic diversity of Pseudomonas aeruginosa across clinical and environmental isolation sites.
Pseudomonas aeruginosa is a clinically significant opportunistic pathogen adept at thriving in both host-associated and environmental settings. To define the extent to which P. aeruginosa isolates specialize across niches and identify genotype-phenotype correlates, we performed whole genome sequencing and comprehensive phenotypic characterization of 125 P. aeruginosa isolates from diverse clinical and environmental sites, evaluating virulence-associated traits, including motility, cytotoxicity, biofilm formation, pyocyanin production, and antimicrobial susceptibility. We identify that genomic diversity does not correlate with isolation source or most virulence phenotypes. Instead, we find that the two major P. aeruginosa clades (Groups A and B) delineate phylogeny and cytotoxicity, with Group B strains showing significantly higher cytotoxicity than Group A. Sequence analysis revealed previously uncharacterized alleles of genes encoding type III secretion effector proteins. We observed high variability amongst strains and isolation sources in all four assayed virulence phenotypes. Antimicrobial resistance (AMR) is exclusively observed in clinical isolates, not environmental, reflecting antibiotic exposure-driven selection. Bacterial GWAS revealed a statistically significant association between cytotoxicity and exoU presence, and we identified a novel exoU allelic variant with decreased cytotoxicity, demonstrating that functional diversity within well-characterized virulence factors may still influence pathogenic outcomes. In summary, our analyses of 125 diverse isolates suggest that the ability of P. aeruginosa to thrive across diverse niches is driven by broadly conserved genetic repertoire rather than niche-specific accessory genes.