Characterization and diversity of defense systems in Providencia pathogen

Introduction Providencia species are emerging opportunistic pathogens associated with multidrug-resistant infections, yet their molecular defense mechanisms against phage or mobile genetic elements remain poorly characterized. Methods We present a comprehensive pan-genomic analysis of antiviral defense systems across 73 complete genomes (or chromosomes) of Providencia stuartii (n = 31) and Providencia rettgeri (n = 42), using DefenseFinder and CRISPRCasFinder. We further expanded analysis of contig/scaffold assemblies to confirm conservation of core defense profiles across assembly types. BacMGEnet was employed to derive spacer-MGE interaction networks. Phylogenetic reconstruction and gene gain and loss modeling were performed to assess evolutionary patterns. To validate functionality, we experimentally tested the anti-phage activity of Gabija and Septu in heterologous E. coli assays, including point mutation analysis of conserved residues. Results We reveal a diverse and complex defense repertoire dominated by restriction-modification systems and CRISPR-Cas Class 1 Type I-F, with significant contributions from toxin-antitoxin, GAPS2, PsyrTA, and Mokosh systems. Notably, defense genes are non-randomly distributed, often clustering into genomic islands suggestive of horizontal acquisition. Expanded analysis confirms conservation of core defense profiles across assembly types, supporting the utility of lower-quality data when complete genomes are scarce. Comparative analysis uncovers species-specific differences, with P. rettgeri harboring a higher abundance of non-CRISPR systems. BacMGEnet-derived spacer-MGE interaction networks further highlight species-specific dynamics, dense, hub-driven networks in P. stuartii versus sparser networks in P. rettgeri . Correlation analysis indicates potential associations between specific defense systems and virulence or antibiotic resistance genes. Phylogenetic reconstruction and gene gain and loss modeling further highlight dynamic evolutionary patterns. Both Gabija and Septu systems conferred robust, phage-specific protection; point mutations in conserved residues (GajA E465K and PtuB H53K) abolished defense. Discussion Our findings unveil a multi-layered, modular immune architecture in Providencia, providing crucial insights into its genome plasticity, phage resistance, and adaptation in clinical environments. This work establishes a foundation for understanding the role of defense systems in the evolution and pathogenicity of the Providencia genus.