Going through phages: a computational approach to revealing the role of prophage in Staphylococcus aureus. - Related Documents




#
Rank
Similarity
Title + Abs.
Year
PMID
012345
434701.0000Going through phages: a computational approach to revealing the role of prophage in Staphylococcus aureus. Prophages have important roles in virulence, antibiotic resistance, and genome evolution in Staphylococcus aureus . Rapid growth in the number of sequenced S. aureus genomes allows for an investigation of prophage sequences at an unprecedented scale. We developed a novel computational pipeline for phage discovery and annotation. We combined PhiSpy, a phage discovery tool, with VGAS and PROKKA, genome annotation tools to detect and analyse prophage sequences in nearly 10 011  S . aureus genomes, discovering thousands of putative prophage sequences with genes encoding virulence factors and antibiotic resistance. To our knowledge, this is the first large-scale application of PhiSpy on a large-scale set of genomes (10 011  S . aureus ). Determining the presence of virulence and resistance encoding genes in prophage has implications for the potential transfer of these genes/functions to other bacteria via transduction and thus can provide insight into the evolution and spread of these genes/functions between bacterial strains. While the phage we have identified may be known, these phages were not necessarily known or characterized in S. aureus and the clustering and comparison we did for phage based on their gene content is novel. Moreover, the reporting of these genes with the S. aureus genomes is novel.202337424556
966210.9998Species-Scale Genomic Analysis of Staphylococcus aureus Genes Influencing Phage Host Range and Their Relationships to Virulence and Antibiotic Resistance Genes. Phage therapy has been proposed as a possible alternative treatment for infections caused by the ubiquitous bacterial pathogen Staphylococcus aureus. However, successful therapy requires understanding the genetic basis of host range-the subset of strains in a species that could be killed by a particular phage. We searched diverse sets of S. aureus public genome sequences against a database of genes suggested from prior studies to influence host range to look for patterns of variation across the species. We found that genes encoding biosynthesis of molecules that were targets of S. aureus phage adsorption to the outer surface of the cell were the most conserved in the pangenome. Putative phage resistance genes that were core components of the pangenome genes had similar nucleotide diversity, ratio of nonsynonymous to synonymous substitutions, and functionality (measured by delta-bitscore) to other core genes. However, phage resistance genes that were not part of the core genome were significantly less consistent with the core genome phylogeny than all noncore genes in this set, suggesting more frequent movement between strains by horizontal gene transfer. Only superinfection immunity genes encoded by temperate phages inserted in the genome correlated with experimentally determined temperate phage resistance. Taken together, these results suggested that, while phage adsorption genes are heavily conserved in the S. aureus species, HGT may play a significant role in strain-specific evolution of host range patterns. IMPORTANCE Staphylococcus aureus is a widespread, hospital- and community-acquired pathogen that is commonly antibiotic resistant. It causes diverse diseases affecting both the skin and internal organs. Its ubiquity, antibiotic resistance, and disease burden make new therapies urgent, such as phage therapy, in which viruses specific to infecting bacteria clear infection. S. aureus phage host range not only determines whether phage therapy will be successful by killing bacteria but also horizontal gene transfer through transduction of host genetic material by phages. In this work, we comprehensively reviewed existing literature to build a list of S. aureus phage resistance genes and searched our database of almost 43,000 S. aureus genomes for these genes to understand their patterns of evolution, finding that prophages' superinfection immunity correlates best with phage resistance and HGT. These findings improved our understanding of the relationship between known phage resistance genes and phage host range in the species.202235040700
437520.9998Evidence of a large novel gene pool associated with prokaryotic genomic islands. Microbial genes that are "novel" (no detectable homologs in other species) have become of increasing interest as environmental sampling suggests that there are many more such novel genes in yet-to-be-cultured microorganisms. By analyzing known microbial genomic islands and prophages, we developed criteria for systematic identification of putative genomic islands (clusters of genes of probable horizontal origin in a prokaryotic genome) in 63 prokaryotic genomes, and then characterized the distribution of novel genes and other features. All but a few of the genomes examined contained significantly higher proportions of novel genes in their predicted genomic islands compared with the rest of their genome (Paired t test = 4.43E-14 to 1.27E-18, depending on method). Moreover, the reverse observation (i.e., higher proportions of novel genes outside of islands) never reached statistical significance in any organism examined. We show that this higher proportion of novel genes in predicted genomic islands is not due to less accurate gene prediction in genomic island regions, but likely reflects a genuine increase in novel genes in these regions for both bacteria and archaea. This represents the first comprehensive analysis of novel genes in prokaryotic genomic islands and provides clues regarding the origin of novel genes. Our collective results imply that there are different gene pools associated with recently horizontally transmitted genomic regions versus regions that are primarily vertically inherited. Moreover, there are more novel genes within the gene pool associated with genomic islands. Since genomic islands are frequently associated with a particular microbial adaptation, such as antibiotic resistance, pathogen virulence, or metal resistance, this suggests that microbes may have access to a larger "arsenal" of novel genes for adaptation than previously thought.200516299586
437430.9998Core genes can have higher recombination rates than accessory genes within global microbial populations. Recombination is essential to microbial evolution, and is involved in the spread of antibiotic resistance, antigenic variation, and adaptation to the host niche. However, assessing the impact of homologous recombination on accessory genes which are only present in a subset of strains of a given species remains challenging due to their complex phylogenetic relationships. Quantifying homologous recombination for accessory genes (which are important for niche-specific adaptations) in comparison to core genes (which are present in all strains and have essential functions) is critical to understanding how selection acts on variation to shape species diversity and genome structures of bacteria. Here, we apply a computationally efficient, non-phylogenetic approach to measure homologous recombination rates in the core and accessory genome using >100,000 whole genome sequences from Streptococcus pneumoniae and several additional species. By analyzing diverse sets of sequence clusters, we show that core genes often have higher recombination rates than accessory genes, and for some bacterial species the associated effect sizes for these differences are pronounced. In a subset of species, we find that gene frequency and homologous recombination rate are positively correlated. For S. pneumoniae and several additional species, we find that while the recombination rate is higher for the core genome, the mutational divergence is lower, indicating that divergence-based homologous recombination barriers could contribute to differences in recombination rates between the core and accessory genome. Homologous recombination may therefore play a key role in increasing the efficiency of selection in the most conserved parts of the genome.202235801696
940440.9998The Application of Transposon Insertion Sequencing in Identifying Essential Genes in B. fragilis. Essential genes are those that are indispensable for the survival of organism under specific growth conditions. Investigating essential genes in pathogenic bacteria not only helps to understand vital biological networks but also provides novel targets for drug development. Availability of genetic engineering tools and high-throughput sequencing methods has enabled essential genes identification in many pathogenic gram-positive and gram-negative bacteria. Bacteroides fragilis is one of the major bacteria specific of human gastrointestinal microbiota. When B. fragilis moves out of its niche, it turns into deadly pathogen. Here, we describe detailed method for the essential gene identification in B. fragilis. Generated transposon mutant pool can be used for other applications such as identification of genes responsible for drug resistance in B. fragilis.202234709623
940550.9998Functional Metagenomic Screening for Antimicrobial Resistance in the Oral Microbiome. A large proportion of bacteria, from a multitude of environments, are not yet able to be grown in the laboratory, and therefore microbiological and molecular biological investigations of these bacteria are challenging. A way to circumvent this challenge is to analyze the metagenome, the entire collection of DNA molecules that can be isolated from a particular environment or sample. This collection of DNA molecules can be sequenced and assembled to determine what is present and infer functional potential, or used as a PCR template to detect known target DNA and potentially unknown regions of DNA nearby those targets; however assigning functions to new or conserved hypothetical, functionally cryptic, genes is difficult. Functional metagenomics allows researchers to determine which genes are responsible for selectable phenotypes, such as resistance to antimicrobials and metabolic capabilities, without the prerequisite needs to grow the bacteria containing those genes or to already know which genes are of interest. It is estimated that a third of the resident species of the human oral cavity is not yet cultivable and, together with the ease of sample acquisition, makes this metagenome particularly suited to functional metagenomic studies. Here we describe the methodology related to the collection of saliva samples, extraction of metagenomic DNA, construction of metagenomic libraries, as well as the description of functional assays that have previously led to the identification of new genes conferring antimicrobial resistance.202134410638
437760.9997Pathogenicity and other genomic islands in plant pathogenic bacteria. SUMMARY Pathogenicity islands (PAIs) were first described in uropathogenic E. coli. They are now defined as regions of DNA that contain virulence genes and are present in the genome of pathogenic strains, but absent from or only rarely present in non-pathogenic variants of the same or related strains. Other features include a variable G+C content, distinct boundaries from the rest of the genome and the presence of genes related to mobile elements such as insertion sequences, integrases and transposases. Although PAIs have now been described in a wide range of both plant and animal pathogens it has become evident that the general features of PAIs are displayed by a number of regions of DNA with functions other than pathogenicity, such as symbiosis and antibiotic resistance, and the general term genomic islands has been adopted. This review will describe a range of genomic islands in plant pathogenic bacteria including those that carry effector genes, phytotoxins and the type III protein secretion cluster. The review will also consider some medically important bacteria in order to discuss the range, acquisition and stabilization of genomic islands.200320569400
966370.9997The structure of temperate phage-bacteria infection networks changes with the phylogenetic distance of the host bacteria. With their ability to integrate into the bacterial chromosome and thereby transfer virulence or drug-resistance genes across bacterial species, temperate phage play a key role in bacterial evolution. Thus, it is paramount to understand who infects whom to be able to predict the movement of DNA across the prokaryotic world and ultimately the emergence of novel (drug-resistant) pathogens. We empirically investigated lytic infection patterns among Vibrio spp. from distinct phylogenetic clades and their derived temperate phage. We found that across distantly related clades, infections occur preferentially within modules of the same clade. However, when the genetic distance of the host bacteria decreases, these clade-specific infections disappear. This indicates that the structure of temperate phage-bacteria infection networks changes with the phylogenetic distance of the host bacteria.201830429242
438680.9997Large-scale screening of a targeted Enterococcus faecalis mutant library identifies envelope fitness factors. Spread of antibiotic resistance among bacteria responsible for nosocomial and community-acquired infections urges for novel therapeutic or prophylactic targets and for innovative pathogen-specific antibacterial compounds. Major challenges are posed by opportunistic pathogens belonging to the low GC% gram-positive bacteria. Among those, Enterococcus faecalis is a leading cause of hospital-acquired infections associated with life-threatening issues and increased hospital costs. To better understand the molecular properties of enterococci that may be required for virulence, and that may explain the emergence of these bacteria in nosocomial infections, we performed the first large-scale functional analysis of E. faecalis V583, the first vancomycin-resistant isolate from a human bloodstream infection. E. faecalis V583 is within the high-risk clonal complex 2 group, which comprises mostly isolates derived from hospital infections worldwide. We conducted broad-range screenings of candidate genes likely involved in host adaptation (e.g., colonization and/or virulence). For this purpose, a library was constructed of targeted insertion mutations in 177 genes encoding putative surface or stress-response factors. Individual mutants were subsequently tested for their i) resistance to oxidative stress, ii) antibiotic resistance, iii) resistance to opsonophagocytosis, iv) adherence to the human colon carcinoma Caco-2 epithelial cells and v) virulence in a surrogate insect model. Our results identified a number of factors that are involved in the interaction between enterococci and their host environments. Their predicted functions highlight the importance of cell envelope glycopolymers in E. faecalis host adaptation. This study provides a valuable genetic database for understanding the steps leading E. faecalis to opportunistic virulence.201122194979
426190.9997Recovery and Characterization of Bacteria Resisting Infection by Lytic Bacteriophage. Bacteria and bacteriophages coexist and coevolve, bacteriophages being obligatory predators exerting an evolutionary pressure on their prey. Mechanisms in action vary depending on the bacterial genomic content and on the regulation of the bacteriophage cycle. To assess the multiplicity of bacterial genes involved in resistance as well as the changes in the bacteriophage interactions with the bacteria, it is necessary to isolate and investigate large numbers of independent resistant variants. Here we describe protocols that have been applied to the study of Pseudomonas aeruginosa and four of its virulent bacteriophages belonging to the Podoviridae and Myoviridae bacteriophage families. Mutations are identified using whole genome sequencing of resistant variants. Phenotypic analyses are performed to describe the changes conferred by the mutations.201829119434
4344100.9997Phenetic Comparison of Prokaryotic Genomes Using k-mers. Bacterial genomics studies are getting more extensive and complex, requiring new ways to envision analyses. Using the Ray Surveyor software, we demonstrate that comparison of genomes based on their k-mer content allows reconstruction of phenetic trees without the need of prior data curation, such as core genome alignment of a species. We validated the methodology using simulated genomes and previously published phylogenomic studies of Streptococcus pneumoniae and Pseudomonas aeruginosa. We also investigated the relationship of specific genetic determinants with bacterial population structures. By comparing clusters from the complete genomic content of a genome population with clusters from specific functional categories of genes, we can determine how the population structures are correlated. Indeed, the strain clustering based on a subset of k-mers allows determination of its similarity with the whole genome clusters. We also applied this methodology on 42 species of bacteria to determine the correlational significance of five important bacterial genomic characteristics. For example, intrinsic resistance is more important in P. aeruginosa than in S. pneumoniae, and the former has increased correlation of its population structure with antibiotic resistance genes. The global view of the pangenome of bacteria also demonstrated the taxa-dependent interaction of population structure with antibiotic resistance, bacteriophage, plasmid, and mobile element k-mer data sets.201728957508
8855110.9997Transposon Insertion Sequencing Elucidates Novel Gene Involvement in Susceptibility and Resistance to Phages T4 and T7 in Escherichia coli O157. Experiments using bacteriophage (phage) to infect bacterial strains have helped define some basic genetic concepts in microbiology, but our understanding of the complexity of bacterium-phage interactions is still limited. As the global threat of antibiotic resistance continues to increase, phage therapy has reemerged as an attractive alternative or supplement to treating antibiotic-resistant bacterial infections. Further, the long-used method of phage typing to classify bacterial strains is being replaced by molecular genetic techniques. Thus, there is a growing need for a complete understanding of the precise molecular mechanisms underpinning phage-bacterium interactions to optimize phage therapy for the clinic as well as for retrospectively interpreting phage typing data on the molecular level. In this study, a genomics-based fitness assay (TraDIS) was used to identify all host genes involved in phage susceptibility and resistance for a T4 phage infecting Shiga-toxigenic Escherichia coli O157. The TraDIS results identified both established and previously unidentified genes involved in phage infection, and a subset were confirmed by site-directed mutagenesis and phenotypic testing of 14 T4 and 2 T7 phages. For the first time, the entire sap operon was implicated in phage susceptibility and, conversely, the stringent starvation protein A gene (sspA) was shown to provide phage resistance. Identifying genes involved in phage infection and replication should facilitate the selection of bespoke phage combinations to target specific bacterial pathogens.IMPORTANCE Antibiotic resistance has diminished treatment options for many common bacterial infections. Phage therapy is an alternative option that was once popularly used across Europe to kill bacteria within humans. Phage therapy acts by using highly specific viruses (called phages) that infect and lyse certain bacterial species to treat the infection. Whole-genome sequencing has allowed modernization of the investigations into phage-bacterium interactions. Here, using E. coli O157 and T4 bacteriophage as a model, we have exploited a genome-wide fitness assay to investigate all genes involved in defining phage resistance or susceptibility. This knowledge of the genetic determinants of phage resistance and susceptibility can be used to design bespoke phage combinations targeted to specific bacterial infections for successful infection eradication.201830042196
4376120.9997Genetic exchanges are more frequent in bacteria encoding capsules. Capsules allow bacteria to colonize novel environments, to withstand numerous stresses, and to resist antibiotics. Yet, even though genetic exchanges with other cells should be adaptive under such circumstances, it has been suggested that capsules lower the rates of homologous recombination and horizontal gene transfer. We analysed over one hundred pan-genomes and thousands of bacterial genomes for the evidence of an association between genetic exchanges (or lack thereof) and the presence of a capsule system. We found that bacteria encoding capsules have larger pan-genomes, higher rates of horizontal gene transfer, and higher rates of homologous recombination in their core genomes. Accordingly, genomes encoding capsules have more plasmids, conjugative elements, transposases, prophages, and integrons. Furthermore, capsular loci are frequent in plasmids, and can be found in prophages. These results are valid for Bacteria, independently of their ability to be naturally transformable. Since we have shown previously that capsules are commonly present in nosocomial pathogens, we analysed their co-occurrence with antibiotic resistance genes. Genomes encoding capsules have more antibiotic resistance genes, especially those encoding efflux pumps, and they constitute the majority of the most worrisome nosocomial bacteria. We conclude that bacteria with capsule systems are more genetically diverse and have fast-evolving gene repertoires, which may further contribute to their success in colonizing novel niches such as humans under antibiotic therapy.201830576310
4637130.9997What Differentiates Probiotic from Pathogenic Bacteria? The Genetic Mobility of Enterococcus faecium Offers New Molecular Insights. Enterococcus faecium is a lactic acid bacterium with applications in food engineering and nutrigenomics, including as starter cultures in fermented foods. To differentiate the E. faecium probiotic from pathogenic bacteria, physiological analyses are often used but they do not guarantee that a bacterial strain is not pathogenic. We report here new findings and an approach based on comparison of the genetic mobility of (1) probiotic, (2) pathogenic, and (3) nonpathogenic and non-probiotic strains, so as to differentiate probiotics, and inform their safe use. The region of the 16S ribosomal DNA (rDNA) genes of different E. faecium strains native to Pernambuco-Brazil was used with the GenBank query sequence. Complete genomes were selected and divided into three groups as noted above to identify the mobile genetic elements (MGEs) (transposase, integrase, conjugative transposon protein and phage) and antibiotic resistance genes (ARGs), and to undertake pan-genome analysis and multiple genome alignment. Differences in the number of MGEs were found in ARGs, in the presence and absence of the genes that differentiate E. faecium probiotics and pathogenic bacteria genetically. Our data suggest that genetic mobility appears to be informative in differentiating between probiotic and pathogenic strains. While the present findings are not necessarily applicable to all probiotics, they offer novel molecular insights to guide future research in nutrigenomics, clinical medicine, and food engineering on new ways to differentiate pathogenic from probiotic bacteria.202032762606
4624140.9997Deciphering the distance to antibiotic resistance for the pneumococcus using genome sequencing data. Advances in genome sequencing technologies and genome-wide association studies (GWAS) have provided unprecedented insights into the molecular basis of microbial phenotypes and enabled the identification of the underlying genetic variants in real populations. However, utilization of genome sequencing in clinical phenotyping of bacteria is challenging due to the lack of reliable and accurate approaches. Here, we report a method for predicting microbial resistance patterns using genome sequencing data. We analyzed whole genome sequences of 1,680 Streptococcus pneumoniae isolates from four independent populations using GWAS and identified probable hotspots of genetic variation which correlate with phenotypes of resistance to essential classes of antibiotics. With the premise that accumulation of putative resistance-conferring SNPs, potentially in combination with specific resistance genes, precedes full resistance, we retrogressively surveyed the hotspot loci and quantified the number of SNPs and/or genes, which if accumulated would confer full resistance to an otherwise susceptible strain. We name this approach the 'distance to resistance'. It can be used to identify the creep towards complete antibiotics resistance in bacteria using genome sequencing. This approach serves as a basis for the development of future sequencing-based methods for predicting resistance profiles of bacterial strains in hospital microbiology and public health settings.201728205635
9406150.9997Proteomics as the final step in the functional metagenomics study of antimicrobial resistance. The majority of clinically applied antimicrobial agents are derived from natural products generated by soil microorganisms and therefore resistance is likely to be ubiquitous in such environments. This is supported by the fact that numerous clinically important resistance mechanisms are encoded within the genomes of such bacteria. Advances in genomic sequencing have enabled the in silico identification of putative resistance genes present in these microorganisms. However, it is not sufficient to rely on the identification of putative resistance genes, we must also determine if the resultant proteins confer a resistant phenotype. This will require an analysis pipeline that extends from the extraction of environmental DNA, to the identification and analysis of potential resistance genes and their resultant proteins and phenotypes. This review focuses on the application of functional metagenomics and proteomics to study antimicrobial resistance in diverse environments.201525784907
4638160.9997Comprehensive Scanning of Prophages in Lactobacillus: Distribution, Diversity, Antibiotic Resistance Genes, and Linkages with CRISPR-Cas Systems. Prophage integration, release, and dissemination exert various effects on host bacteria. In the genus Lactobacillus, they may cause bacteriophage contamination during fermentation and even regulate bacterial populations in the gut. However, little is known about their distribution, genetic architecture, and relationships with their hosts. Here, we conducted prophage prediction analysis on 1,472 genomes from 16 different Lactobacillus species and found prophage fragments in almost all lactobacilli (99.8%), with 1,459 predicted intact prophages identified in 64.1% of the strains. We present an uneven prophage distribution among Lactobacillus species; multihabitat species retained more prophages in their genomes than restricted-habitat species. Characterization of the genome features, average nucleotide identity, and landscape visualization presented a high genome diversity of Lactobacillus prophages. We detected antibiotic resistance genes in more than 10% of Lactobacillus prophages and validated that the occurrence of resistance genes conferred by prophage integration was possibly associated with phenotypic resistance in Lactobacillus plantarum. Furthermore, our broad and comprehensive examination of the distribution of CRISPR-Cas systems across the genomes predicted type I and type III systems as potential antagonistic elements of Lactobacillus prophage. IMPORTANCE Lactobacilli are inherent microorganisms in the human gut and are widely used in the food processing industries due to their probiotic properties. Prophages were reportedly hidden in numerous Lactobacillus genomes and can potentially contaminate entire batches of fermentation or modulate the intestinal microecology once they are released. Therefore, a comprehensive scanning of prophages in Lactobacillus is essential for the safety evaluation and application development of probiotic candidates. We show that prophages are widely distributed among lactobacilli; however, intact prophages are more common in multihabitat species and display wide variations in genome feature, integration site, and genomic organization. Our data of the prophage-mediated antibiotic resistance genes (ARGs) and the resistance phenotype of lactobacilli provide evidence for deciphering the putative role of prophages as vectors of the ARGs. Furthermore, understanding the association between prophages and CRISPR-Cas systems is crucial to appreciate the coevolution of phages and Lactobacillus.202134060909
4381170.9997Specific Gene Loci of Clinical Pseudomonas putida Isolates. Pseudomonas putida are ubiquitous inhabitants of soils and clinical isolates of this species have been seldom described. Clinical isolates show significant variability in their ability to cause damage to hosts because some of them are able to modulate the host's immune response. In the current study, comparisons between the genomes of different clinical and environmental strains of P. putida were done to identify genetic clusters shared by clinical isolates that are not present in environmental isolates. We show that in clinical strains specific genes are mostly present on transposons, and that this set of genes exhibit high identity with genes found in pathogens and opportunistic pathogens. The set of genes prevalent in P. putida clinical isolates, and absent in environmental isolates, are related with survival under oxidative stress conditions, resistance against biocides, amino acid metabolism and toxin/antitoxin (TA) systems. This set of functions have influence in colonization and survival within human tissues, since they avoid host immune response or enhance stress resistance. An in depth bioinformatic analysis was also carried out to identify genetic clusters that are exclusive to each of the clinical isolates and that correlate with phenotypical differences between them, a secretion system type III-like was found in one of these clinical strains, a determinant of pathogenicity in Gram-negative bacteria.201626820467
9645180.9997Horizontal Gene Transfers in prokaryotes show differential preferences for metabolic and translational genes. BACKGROUND: Horizontal gene transfer (HGT) is an important process, which contributes in bacterial pathogenesis and drug resistance. A number of methods have been proposed for detection of horizontal gene transfer. One successful approach to the detection of HGT events is due to Novichkov et al. (J. Bacteriology 186, 6575-85), who rely on comparing phylogenetic distances within a gene family with genomic distances of the source organisms. Building on their approach, we introduce outlier detection in the correlation between those two sets of distances. This approach is designed to detect horizontal transfers of core set of genes present in many bacteria. The principle behind method allows detection of xenologous gene displacements as well as acquisition of novel genes. RESULTS: Simulations indicated that our method performs better than Novichkov et al's original approach. The approach very efficiently identified HGT between distantly related bacteria and also a limited number of gene transfers between closely related bacteria. In combination with sequence similarity and likelihood tests, it yields a measure robust enough to derive a set of 171 genes deemed likely to have been horizontally transferred. Further analysis of these 171 established horizontal transfer events gave interesting insights in the direction of transfer. CONCLUSION: The majority of transfers between archaea and bacteria have occurred in the direction from bacteria to archaea rather than the other way round. Genes transferred between the archaea and bacteria are mostly metabolic genes. On the other hand, genes transferred within the bacterial phyla are mainly involved in translation.200919134215
3829190.9997Associations among Antibiotic and Phage Resistance Phenotypes in Natural and Clinical Escherichia coli Isolates. The spread of antibiotic resistance is driving interest in new approaches to control bacterial pathogens. This includes applying multiple antibiotics strategically, using bacteriophages against antibiotic-resistant bacteria, and combining both types of antibacterial agents. All these approaches rely on or are impacted by associations among resistance phenotypes (where bacteria resistant to one antibacterial agent are also relatively susceptible or resistant to others). Experiments with laboratory strains have shown strong associations between some resistance phenotypes, but we lack a quantitative understanding of associations among antibiotic and phage resistance phenotypes in natural and clinical populations. To address this, we measured resistance to various antibiotics and bacteriophages for 94 natural and clinical Escherichia coli isolates. We found several positive associations between resistance phenotypes across isolates. Associations were on average stronger for antibacterial agents of the same type (antibiotic-antibiotic or phage-phage) than different types (antibiotic-phage). Plasmid profiles and genetic knockouts suggested that such associations can result from both colocalization of resistance genes and pleiotropic effects of individual resistance mechanisms, including one case of antibiotic-phage cross-resistance. Antibiotic resistance was predicted by core genome phylogeny and plasmid profile, but phage resistance was predicted only by core genome phylogeny. Finally, we used observed associations to predict genes involved in a previously uncharacterized phage resistance mechanism, which we verified using experimental evolution. Our data suggest that susceptibility to phages and antibiotics are evolving largely independently, and unlike in experiments with lab strains, negative associations between antibiotic resistance phenotypes in nature are rare. This is relevant for treatment scenarios where bacteria encounter multiple antibacterial agents.IMPORTANCE Rising antibiotic resistance is making it harder to treat bacterial infections. Whether resistance to a given antibiotic spreads or declines is influenced by whether it is associated with altered susceptibility to other antibiotics or other stressors that bacteria encounter in nature, such as bacteriophages (viruses that infect bacteria). We used natural and clinical isolates of Escherichia coli, an abundant species and key pathogen, to characterize associations among resistance phenotypes to various antibiotics and bacteriophages. We found associations between some resistance phenotypes, and in contrast to past work with laboratory strains, they were exclusively positive. Analysis of bacterial genome sequences and horizontally transferred genetic elements (plasmids) helped to explain this, as well as our finding that there was no overall association between antibiotic resistance and bacteriophage resistance profiles across isolates. This improves our understanding of resistance evolution in nature, potentially informing new rational therapies that combine different antibacterials, including bacteriophages.201729089428