Impact of mutagenesis and lateral gene transfer processes in bacterial susceptibility to phage in food biocontrol and phage therapy.

Introduction: The emergence of resistance and interference mechanisms to phage infection can hinder the success of bacteriophage-based applications, but the significance of these mechanisms in phage therapy has not been determined. This work studies the emergence of Salmonella isolates with reduced susceptibility to a cocktail of three phages under three scenarios: i) Salmonella cultures (LAB), ii) biocontrol of cooked ham slices as a model of food safety (FOOD), and iii) oral phage therapy in broilers (PT). Methods: S. Typhimurium ATCC 14028 RifR variants with reduced phage susceptibility were isolated from the three scenarios and conventional and molecular microbiology techniques were applied to study them. Results and discussion: In LAB, 92% of Salmonella isolates lost susceptibility to all three phages 24 h after phage infection. This percentage was lower in FOOD, with 4.3% of isolates not susceptible to at least two of the three phages after seven days at 4°C following phage treatment. In PT, 9.7% and 3.3 % of isolates from untreated and treated broilers, respectively, displayed some mechanism of interference with the life cycle of some of the phages. In LAB and FOOD scenarios, resistant variants carrying mutations in rfc and rfaJ genes involved in lipopolysaccharide synthesis (phage receptor) were identified. However, in PT, the significant decrease of EOP, ECOI, and burst size observed in isolates was prompted by lateral gene transfer of large IncI1 plasmids, which may encode phage defense mechanisms. These data indicate that the acquisition of specific conjugative plasmids has a stronger impact than mutagenesis on the emergence of reduced phage-susceptibility bacteria in certain environments. In spite of this, neither mechanism seems to significantly impair the success of Salmonella biocontrol and oral phage therapy.

Systematic In Silico Assessment of Antimicrobial Resistance Dissemination across the Global Plasmidome.

The emergence of pathogenic strains resistant to multiple antimicrobials is a pressing problem in modern healthcare. Antimicrobial resistance is mediated primarily by dissemination of resistance determinants via horizontal gene transfer. The dissemination of some resistance genes has been well documented, but few studies have analyzed the patterns underpinning the dissemination of antimicrobial resistance genes. Analyzing the %GC content of plasmid-borne antimicrobial resistance genes relative to their host genome %GC content provides a means to efficiently detect and quantify dissemination of antimicrobial resistance genes. In this work we automate %GC content analysis to perform a comprehensive analysis of known antimicrobial resistance genes in publicly available plasmid sequences. We find that the degree to which antimicrobial resistance genes are disseminated depends primarily on the resistance mechanism. Our analysis identifies conjugative plasmids as primary dissemination vectors and indicates that most broadly disseminated genes have spread from single genomic backgrounds. We show that resistance dissemination profiles vary greatly among antimicrobials, oftentimes reflecting stewardship measures. Our findings establish %GC content analysis as a powerful, intuitive and scalable method to monitor the dissemination of resistance determinants using publicly available sequence data.

Non-canonical LexA proteins regulate the SOS response in the Bacteroidetes.

Lesions to DNA compromise chromosome integrity, posing a direct threat to cell survival. The bacterial SOS response is a widespread transcriptional regulatory mechanism to address DNA damage. This response is coordinated by the LexA transcriptional repressor, which controls genes involved in DNA repair, mutagenesis and cell-cycle control. To date, the SOS response has been characterized in most major bacterial groups, with the notable exception of the Bacteroidetes. No LexA homologs had been identified in this large, diverse and ecologically important phylum, suggesting that it lacked an inducible mechanism to address DNA damage. Here, we report the identification of a novel family of transcriptional repressors in the Bacteroidetes that orchestrate a canonical response to DNA damage in this phylum. These proteins belong to the S24 peptidase family, but are structurally different from LexA. Their N-terminal domain is most closely related to CI-type bacteriophage repressors, suggesting that they may have originated from phage lytic phase repressors. Given their role as SOS regulators, however, we propose to designate them as non-canonical LexA proteins. The identification of a new class of repressors orchestrating the SOS response illuminates long-standing questions regarding the origin and plasticity of this transcriptional network.

Importance of twitching and surface-associated motility in the virulence of Acinetobacter baumannii.

Acinetobacter baumannii is a pathogen of increasing clinical importance worldwide, especially given its ability to readily acquire resistance determinants. Motile strains of this bacterium can move by either or both of two types of motility: (i) twitching, driven by type IV pili, and (ii) surface-associated motility, an appendage-independent form of movement. A. baumannii strain MAR002 possesses both twitching and surface-associated motility. In this study, we isolated spontaneous rifampin-resistant mutants of strain MAR002 in which point mutations in the rpoB gene were identified that resulted in an altered motility pattern. Transcriptomic analysis of mutants lacking twitching, surface-associated motility, or both led to the identification of deregulated genes within each motility phenotype, based on their level of expression and their biological function. Investigations of the corresponding knockout mutants revealed several genes involved in the motility of A. baumannii strain MAR002, including two involved in twitching (encoding a minor pilin subunit and an RND [resistance nodulation division] component), one in surface-associated motility (encoding an amino acid permease), and eight in both (encoding RND and ABC components, the energy transducer TonB, the porin OprD, the T6SS component TagF, an IclR transcriptional regulator, a PQQ-dependent sugar dehydrogenase, and a putative pectate lyase). Virulence assays showed the reduced pathogenicity of mutants with impairments in both types of motility or in surface-associated motility alone. By contrast, the virulence of twitching-affected mutants was not affected. These results shed light on the key role of surface-associated motility and the limited role of twitching in the pathogenicity of A. baumannii.

Flexible comparative genomics of prokaryotic transcriptional regulatory networks.

BACKGROUND: Comparative genomics methods enable the reconstruction of bacterial regulatory networks using available experimental data. In spite of their potential for accelerating research into the composition and evolution of bacterial regulons, few comparative genomics suites have been developed for the automated analysis of these regulatory systems. Available solutions typically rely on precomputed databases for operon and ortholog predictions, limiting the scope of analyses to processed complete genomes, and several key issues such as the transfer of experimental information or the integration of regulatory information in a probabilistic setting remain largely unaddressed. RESULTS: Here we introduce CGB, a flexible platform for comparative genomics of prokaryotic regulons. CGB has few external dependencies and enables fully customized analyses of newly available genome data. The platform automates the merging of experimental information and uses a gene-centered, Bayesian framework to generate and integrate easily interpretable results. We demonstrate its flexibility and power by analyzing the evolution of type III secretion system regulation in pathogenic Proteobacteria and by characterizing the SOS regulon of a new bacterial phylum, the Balneolaeota. CONCLUSIONS: Our results demonstrate the applicability of the CGB pipeline in multiple settings. CGB's ability to automatically integrate experimental information from multiple sources and use complete and draft genomic data, coupled with its non-reliance on precomputed databases and its easily interpretable display of gene-centered posterior probabilities of regulation provide users with an unprecedented level of flexibility in launching comparative genomics analyses of prokaryotic transcriptional regulatory networks. The analyses of type III secretion and SOS response regulatory networks illustrate instances of convergent and divergent evolution of these regulatory systems, showcasing the power of formal ancestral state reconstruction at inferring the evolutionary history of regulatory networks.

Exploration into the origins and mobilization of di-hydrofolate reductase genes and the emergence of clinical resistance to trimethoprim.

Trimethoprim is a synthetic antibacterial agent that targets folate biosynthesis by competitively binding to the di-hydrofolate reductase enzyme (DHFR). Trimethoprim is often administered synergistically with sulfonamide, another chemotherapeutic agent targeting the di-hydropteroate synthase (DHPS) enzyme in the same pathway. Clinical resistance to both drugs is widespread and mediated by enzyme variants capable of performing their biological function without binding to these drugs. These mutant enzymes were assumed to have arisen after the discovery of these synthetic drugs, but recent work has shown that genes conferring resistance to sulfonamide were present in the bacterial pangenome millions of years ago. Here, we apply phylogenetics and comparative genomics methods to study the largest family of mobile trimethoprim-resistance genes (dfrA). We show that most of the dfrA genes identified to date map to two large clades that likely arose from independent mobilization events. In contrast to sulfonamide resistance (sul) genes, we find evidence of recurrent mobilization in dfrA genes. Phylogenetic evidence allows us to identify novel dfrA genes in the emerging pathogen Acinetobacter baumannii, and we confirm their resistance phenotype in vitro. We also identify a cluster of dfrA homologues in cryptic plasmid and phage genomes, but we show that these enzymes do not confer resistance to trimethoprim. Our methods also allow us to pinpoint the chromosomal origin of previously reported dfrA genes, and we show that many of these ancient chromosomal genes also confer resistance to trimethoprim. Our work reveals that trimethoprim resistance predated the clinical use of this chemotherapeutic agent, but that novel mutations have likely also arisen and become mobilized following its widespread use within and outside the clinic. Hence, this work confirms that resistance to novel drugs may already be present in the bacterial pangenome, and stresses the importance of rapid mobilization as a fundamental element in the emergence and global spread of resistance determinants.

Origin of the Mobile Di-Hydro-Pteroate Synthase Gene Determining Sulfonamide Resistance in Clinical Isolates.

Sulfonamides are synthetic chemotherapeutic agents that work as competitive inhibitors of the di-hydro-pteroate synthase (DHPS) enzyme, encoded by the folP gene. Resistance to sulfonamides is widespread in the clinical setting and predominantly mediated by plasmid- and integron-borne sul1-3 genes encoding mutant DHPS enzymes that do not bind sulfonamides. In spite of their clinical importance, the genetic origin of sul1-3 genes remains unknown. Here we analyze sul genes and their genetic neighborhoods to uncover sul signature elements that enable the elucidation of their genetic origin. We identify a protein sequence Sul motif associated with sul-encoded proteins, as well as consistent association of a phosphoglucosamine mutase gene (glmM) with the sul2 gene. We identify chromosomal folP genes bearing these genetic markers in two bacterial families: the Rhodobiaceae and the Leptospiraceae. Bayesian phylogenetic inference of FolP/Sul and GlmM protein sequences clearly establishes that sul1-2 and sul3 genes originated as a mobilization of folP genes present in, respectively, the Rhodobiaceae and the Leptospiraceae, and indicate that the Rhodobiaceae folP gene was transferred from the Leptospiraceae. Analysis of %GC content in folP/sul gene sequences supports the phylogenetic inference results and indicates that the emergence of the Sul motif in chromosomally encoded FolP proteins is ancient and considerably predates the clinical introduction of sulfonamides. In vitro assays reveal that both the Rhodobiaceae and the Leptospiraceae, but not other related chromosomally encoded FolP proteins confer resistance in a sulfonamide-sensitive Escherichia coli background, indicating that the Sul motif is associated with sulfonamide resistance. Given the absence of any known natural sulfonamides targeting DHPS, these results provide a novel perspective on the emergence of resistance to synthetic chemotherapeutic agents, whereby preexisting resistant variants in the vast bacterial pangenome may be rapidly selected for and disseminated upon the clinical introduction of novel chemotherapeuticals.

Comparative genomics of the DNA damage-inducible network in the Patescibacteria.

Metagenomics provide unprecedented insights into the genetic diversity of uncultivated bacteria inhabiting natural environments. Recent surveys have uncovered a major radiation of candidate phyla encompassing the Patescibacteria superphylum. Patescibacteria have small genomes and a presumed symbiotic or parasitic lifestyle, but the difficulty in culturing representative members constrains the characterization of behavioural and adaptive traits. Here we combine in silico and in vitro approaches to characterize the SOS transcriptional response to DNA damage in the Patescibacteria superphylum. Leveraging comparative genomics methods, we identify and experimentally define a novel binding motif for the SOS transcriptional repressor LexA, and we use this motif to characterize the conserved elements of the SOS regulatory network in Patescibacteria. The Patescibacteria LexA-binding motif has unusual direct-repeat structure, and comparative analyses reveal sequence and structural similarities with the distant Acidobacteria LexA protein. Our results reveal a shared core SOS network, complemented by varying degrees of LexA regulation of other core SOS functions. This work illustrates how the combination of computational and experimental methods can leverage metagenomic data to characterize transcriptional responses in uncultivated bacteria. The report of an operational SOS response in presumed symbiotic and parasitic bacteria hints at an intermediate step in the process of genome reduction.

Genomics of Three New Bacteriophages Useful in the Biocontrol of Salmonella.

Non-typhoid Salmonella is the principal pathogen related to food-borne diseases throughout the world. Widespread antibiotic resistance has adversely affected human health and has encouraged the search for alternative antimicrobial agents. The advances in bacteriophage therapy highlight their use in controlling a broad spectrum of food-borne pathogens. One requirement for the use of bacteriophages as antibacterials is the characterization of their genomes. In this work, complete genome sequencing and molecular analyses were carried out for three new virulent Salmonella-specific bacteriophages (UAB_Phi20, UAB_Phi78, and UAB_Phi87) able to infect a broad range of Salmonella strains. Sequence analysis of the genomes of UAB_Phi20, UAB_Phi78, and UAB_Phi87 bacteriophages did not evidence the presence of known virulence-associated and antibiotic resistance genes, and potential immunoreactive food allergens. The UAB_Phi20 genome comprised 41,809 base pairs with 80 open reading frames (ORFs); 24 of them with assigned function. Genome sequence showed a high homology of UAB_Phi20 with Salmonella bacteriophage P22 and other P22likeviruses genus of the Podoviridae family, including ST64T and ST104. The DNA of UAB_Phi78 contained 44,110 bp including direct terminal repeats (DTR) of 179 bp and 58 putative ORFs were predicted and 20 were assigned function. This bacteriophage was assigned to the SP6likeviruses genus of the Podoviridae family based on its high similarity not only with SP6 but also with the K1-5, K1E, and K1F bacteriophages, all of which infect Escherichia coli. The UAB_Phi87 genome sequence consisted of 87,669 bp with terminal direct repeats of 608 bp; although 148 ORFs were identified, putative functions could be assigned to only 29 of them. Sequence comparisons revealed the mosaic structure of UAB_Phi87 and its high similarity with bacteriophages Felix O1 and wV8 of E. coli with respect to genetic content and functional organization. Phylogenetic analysis of large terminase subunits confirms their packaging strategies and grouping to the different phage genus type. All these studies are necessary for the development and the use of an efficient cocktail with commercial applications in bacteriophage therapy against Salmonella.