# | Rank | Similarity | Title + Abs. | Year | PMID |
|---|---|---|---|---|---|
| 0 | 1 | 2 | 3 | 4 | 5 |
| 3809 | 0 | 1.0000 | High abundance of virulence gene homologues in marine bacteria. Marine bacteria can cause harm to single-celled and multicellular eukaryotes. However, relatively little is known about the underlying genetic basis for marine bacterial interactions with higher organisms. We examined whole-genome sequences from a large number of marine bacteria for the prevalence of homologues to virulence genes and pathogenicity islands known from bacteria that are pathogenic to terrestrial animals and plants. As many as 60 out of 119 genomes of marine bacteria, with no known association to infectious disease, harboured genes of virulence-associated types III, IV, V and VI protein secretion systems. Type III secretion was relatively uncommon, while type IV was widespread among alphaproteobacteria (particularly among roseobacters) and type VI was primarily found among gammaproteobacteria. Other examples included homologues of the Yersinia murine toxin and a phage-related 'antifeeding' island. Analysis of the Global Ocean Sampling metagenomic data indicated that virulence genes were present in up to 8% of the planktonic bacteria, with highest values in productive waters. From a marine ecology perspective, expression of these widely distributed genes would indicate that some bacteria infect or even consume live cells, that is, generate a previously unrecognized flow of organic matter and nutrients directly from eukaryotes to bacteria. | 2009 | 19207573 |
| 3830 | 1 | 0.9998 | Resistance Gene Carriage Predicts Growth of Natural and Clinical Escherichia coli Isolates in the Absence of Antibiotics. Bacterial pathogens that carry antibiotic resistance alleles sometimes pay a cost in the form of impaired growth in antibiotic-free conditions. This cost of resistance is expected to be a key parameter for understanding how resistance spreads and persists in pathogen populations. Analysis of individual resistance alleles from laboratory evolution and natural isolates has shown they are typically costly, but these costs are highly variable and influenced by genetic variation at other loci. It therefore remains unclear how strongly resistance is linked to impaired antibiotic-free growth in bacteria from natural and clinical scenarios, where resistance alleles are likely to coincide with other types of genetic variation. To investigate this, we measured the growth of 92 natural and clinical Escherichia coli isolates across three antibiotic-free environments. We then tested whether variation of antibiotic-free growth among isolates was predicted by their resistance to 10 antibiotics, while accounting for the phylogenetic structure of the data. We found that isolates with similar resistance profiles had similar antibiotic-free growth profiles, but it was not simply that higher average resistance was associated with impaired growth. Next, we used whole-genome sequences to identify antibiotic resistance genes and found that isolates carrying a greater number of resistance gene types grew relatively poorly in antibiotic-free conditions, even when the resistance genes they carried were different. This suggests that the resistance of bacterial pathogens is linked to growth costs in nature, but it is the total genetic burden and multivariate resistance phenotype that predict these costs, rather than individual alleles or mean resistance across antibiotics.IMPORTANCE Managing the spread of antibiotic resistance in bacterial pathogens is a major challenge for global public health. Central to this challenge is understanding whether resistance is linked to impaired bacterial growth in the absence of antibiotics, because this determines whether resistance declines when bacteria are no longer exposed to antibiotics. We studied 92 isolates of the key bacterial pathogen Escherichia coli; these isolates varied in both their antibiotic resistance genes and other parts of the genome. Taking this approach, rather than focusing on individual genetic changes associated with resistance as in much previous work, revealed that growth without antibiotics was linked to the number of specialized resistance genes carried and the combination of antibiotics to which isolates were resistant but was not linked to average antibiotic resistance. This approach provides new insights into the genetic factors driving the long-term persistence of antibiotic-resistant bacteria, which is important for future efforts to predict and manage resistance. | 2019 | 30530714 |
| 4375 | 2 | 0.9998 | Evidence 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. | 2005 | 16299586 |
| 9405 | 3 | 0.9998 | Functional 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. | 2021 | 34410638 |
| 3796 | 4 | 0.9998 | The presence of plasmids in bacterial hosts alters phage isolation and infectivity. Antibiotic resistance genes are often carried by plasmids, which spread intra- and inter genera bacterial populations, and also play a critical role in bacteria conferring phage resistance. However, it remains unknown about the influence of plasmids present in bacterial hosts on phage isolation and subsequent infectivity. In this study, using both Escherichia coli and Pseudomonas putida bacteria containing different plasmids, eight phages were isolated and characterized in terms of phage morphology and host range analysis, in conjunction with DNA and protein sequencing. We found that plasmids can influence both the phage isolation process and phage infectivity. In particular, the isolated phages exhibited different phage plaquing infectivity towards the same bacterial species containing different plasmids. Furthermore, the presence of plasmids was found to alter the expression of bacteria membrane protein, which correlates with bacterial cell surface receptors recognized by phages, thus affecting phage isolation and infectivity. Given the diverse and ubiquitous nature of plasmids, our findings highlight the need to consider plasmids as factors that can influence both phage isolation and infectivity. | 2022 | 37938681 |
| 9655 | 5 | 0.9998 | High genomic diversity of multi-drug resistant wastewater Escherichia coli. Wastewater treatment plants play an important role in the emergence of antibiotic resistance. They provide a hot spot for exchange of resistance within and between species. Here, we analyse and quantify the genomic diversity of the indicator Escherichia coli in a German wastewater treatment plant and we relate it to isolates' antibiotic resistance. Our results show a surprisingly large pan-genome, which mirrors how rich an environment a treatment plant is. We link the genomic analysis to a phenotypic resistance screen and pinpoint genomic hot spots, which correlate with a resistance phenotype. Besides well-known resistance genes, this forward genomics approach generates many novel genes, which correlated with resistance and which are partly completely unknown. A surprising overall finding of our analyses is that we do not see any difference in resistance and pan genome size between isolates taken from the inflow of the treatment plant and from the outflow. This means that while treatment plants reduce the amount of bacteria released into the environment, they do not reduce the potential for antibiotic resistance of these bacteria. | 2018 | 29895899 |
| 4381 | 6 | 0.9998 | Specific 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. | 2016 | 26820467 |
| 4050 | 7 | 0.9998 | Are Virulence and Antibiotic Resistance Genes Linked? A Comprehensive Analysis of Bacterial Chromosomes and Plasmids. Although pathogenic bacteria are the targets of antibiotics, these drugs also affect hundreds of commensal or mutualistic species. Moreover, the use of antibiotics is not only restricted to the treatment of infections but is also largely applied in agriculture and in prophylaxis. During this work, we tested the hypothesis that there is a correlation between the number and the genomic location of antibiotic resistance (AR) genes and virulence factor (VF) genes. We performed a comprehensive study of 16,632 reference bacterial genomes in which we identified and counted all orthologues of AR and VF genes in each of the locations: chromosomes, plasmids, or in both locations of the same genome. We found that, on a global scale, no correlation emerges. However, some categories of AR and VF genes co-occur preferentially, and in the mobilome, which supports the hypothesis that some bacterial pathogens are under selective pressure to be resistant to specific antibiotics, a fact that can jeopardize antimicrobial therapy for some human-threatening diseases. | 2022 | 35740113 |
| 8385 | 8 | 0.9998 | Function and Phylogeny of Bacterial Butyryl Coenzyme A:Acetate Transferases and Their Diversity in the Proximal Colon of Swine. Studying the host-associated butyrate-producing bacterial community is important, because butyrate is essential for colonic homeostasis and gut health. Previous research has identified the butyryl coenzyme A (CoA):acetate-CoA transferase (EC 2.3.8.3) as a gene of primary importance for butyrate production in intestinal ecosystems; however, this gene family (but) remains poorly defined. We developed tools for the analysis of butyrate-producing bacteria based on 12 putative but genes identified in the genomes of nine butyrate-producing bacteria obtained from the swine intestinal tract. Functional analyses revealed that eight of these genes had strong But enzyme activity. When but paralogues were found within a genome, only one gene per genome encoded strong activity, with the exception of one strain in which no gene encoded strong But activity. Degenerate primers were designed to amplify the functional but genes and were tested by amplifying environmental but sequences from DNA and RNA extracted from swine colonic contents. The results show diverse but sequences from swine-associated butyrate-producing bacteria, most of which clustered near functionally confirmed sequences. Here, we describe tools and a framework that allow the bacterial butyrate-producing community to be profiled in the context of animal health and disease. IMPORTANCE: Butyrate is a compound produced by the microbiota in the intestinal tracts of animals. This compound is of critical importance for intestinal health, and yet studying its production by diverse intestinal bacteria is technically challenging. Here, we present an additional way to study the butyrate-producing community of bacteria using one degenerate primer set that selectively targets genes experimentally demonstrated to encode butyrate production. This work will enable researchers to more easily study this very important bacterial function that has implications for host health and resistance to disease. | 2016 | 27613689 |
| 3811 | 9 | 0.9998 | Minor fitness costs in an experimental model of horizontal gene transfer in bacteria. Genes introduced by horizontal gene transfer (HGT) from other species constitute a significant portion of many bacterial genomes, and the evolutionary dynamics of HGTs are important for understanding the spread of antibiotic resistance and the emergence of new pathogenic strains of bacteria. The fitness effects of the transferred genes largely determine the fixation rates and the amount of neutral diversity of newly acquired genes in bacterial populations. Comparative analysis of bacterial genomes provides insight into what genes are commonly transferred, but direct experimental tests of the fitness constraints on HGT are scarce. Here, we address this paucity of experimental studies by introducing 98 random DNA fragments varying in size from 0.45 to 5 kb from Bacteroides, Proteus, and human intestinal phage into a defined position in the Salmonella chromosome and measuring the effects on fitness. Using highly sensitive competition assays, we found that eight inserts were deleterious with selection coefficients (s) ranging from ≈ -0.007 to -0.02 and 90 did not have significant fitness effects. When inducing transcription from a PBAD promoter located at one end of the insert, 16 transfers were deleterious and 82 were not significantly different from the control. In conclusion, a major fraction of the inserts had minor effects on fitness implying that extra DNA transferred by HGT, even though it does not confer an immediate selective advantage, could be maintained at selection-transfer balance and serve as raw material for the evolution of novel beneficial functions. | 2014 | 24536043 |
| 4374 | 10 | 0.9998 | Core 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. | 2022 | 35801696 |
| 3810 | 11 | 0.9998 | The Effect of the Presence and Absence of DNA Repair Genes on the Rate and Pattern of Mutation in Bacteria. Bacteria lose and gain repair genes as they evolve. Here, we investigate the consequences of gain and loss of 11 DNA repair genes across a broad range of bacteria. Using synonymous polymorphisms from bacteria and a set of 50 phylogenetically independent contrasts, we find no evidence that the presence or absence of these 11 genes affects either the overall level of diversity or the pattern of mutation. Using phylogenetic generalized linear squares yields a similar conclusion. It seems likely that the lack of an effect is due to variation in the genetic background and the environment which obscures any effects that the presence or absence of individual genes might have. | 2024 | 39376054 |
| 9658 | 12 | 0.9998 | Functional metagenomic libraries generated from anthropogenically impacted environments reveal importance of metabolic genes in biocide and antibiotic resistance. Anthropogenic activities result in the release of antimicrobial resistant bacteria and a cocktail of antimicrobial compounds into the environment that may directly select or indirectly co-select for antimicrobial resistance (AMR). Many studies use metagenome sequencing or qPCR-based approaches to study the environmental resistome but these methods are limited by a priori knowledge. In this study, a functional metagenomic approach was used to explore biocide resistance mechanisms in two contaminated environments and a pristine site, and to identify whether potentially novel genes conferring biocide resistance also conferred resistance or reduced susceptibility to antibiotics. Resistance was predominately mediated through novel mechanisms exclusive of the well-known qac efflux genes. UDP-galactose 4-epimerase (galE) -like genes were identified in both contaminated environments and were shown to confer cross-resistance to biocides and clinically important antibiotics for the first time (to our knowledge), compared to knockout mutants. GalE -like genes were also co-located with transposons, suggesting mobilisation potential. These results show that housekeeping genes may play a significant yet underappreciated role in AMR in environmental microbiomes. | 2023 | 36908773 |
| 9004 | 13 | 0.9998 | Shedding light on the bacterial resistance to toxic UV filters: a comparative genomic study. UV filters are toxic to marine bacteria that dominate the marine biomass. Ecotoxicology often studies the organism response but rarely integrates the toxicity mechanisms at the molecular level. In this study, in silico comparative genomics between UV filters sensitive and resistant bacteria were conducted in order to unravel the genes responsible for a resistance phenotype. The genomes of two environmentally relevant Bacteroidetes and three Firmicutes species were compared through pairwise comparison. Larger genomes were carried by bacteria exhibiting a resistant phenotype, favoring their ability to adapt to environmental stresses. While the antitoxin and CRISPR systems were the only distinctive features in resistant Bacteroidetes, Firmicutes displayed multiple unique genes that could support the difference between sensitive and resistant phenotypes. Several genes involved in ROS response, vitamin biosynthesis, xenobiotic degradation, multidrug resistance, and lipophilic compound permeability were shown to be exclusive to resistant species. Our investigation contributes to a better understanding of UV filters resistance phenotypes, by identifying pivotal genes involved in key pathways. | 2021 | 34760358 |
| 4665 | 14 | 0.9998 | A comprehensive survey of integron-associated genes present in metagenomes. BACKGROUND: Integrons are genomic elements that mediate horizontal gene transfer by inserting and removing genetic material using site-specific recombination. Integrons are commonly found in bacterial genomes, where they maintain a large and diverse set of genes that plays an important role in adaptation and evolution. Previous studies have started to characterize the wide range of biological functions present in integrons. However, the efforts have so far mainly been limited to genomes from cultivable bacteria and amplicons generated by PCR, thus targeting only a small part of the total integron diversity. Metagenomic data, generated by direct sequencing of environmental and clinical samples, provides a more holistic and unbiased analysis of integron-associated genes. However, the fragmented nature of metagenomic data has previously made such analysis highly challenging. RESULTS: Here, we present a systematic survey of integron-associated genes in metagenomic data. The analysis was based on a newly developed computational method where integron-associated genes were identified by detecting their associated recombination sites. By processing contiguous sequences assembled from more than 10 terabases of metagenomic data, we were able to identify 13,397 unique integron-associated genes. Metagenomes from marine microbial communities had the highest occurrence of integron-associated genes with levels more than 100-fold higher than in the human microbiome. The identified genes had a large functional diversity spanning over several functional classes. Genes associated with defense mechanisms and mobility facilitators were most overrepresented and more than five times as common in integrons compared to other bacterial genes. As many as two thirds of the genes were found to encode proteins of unknown function. Less than 1% of the genes were associated with antibiotic resistance, of which several were novel, previously undescribed, resistance gene variants. CONCLUSIONS: Our results highlight the large functional diversity maintained by integrons present in unculturable bacteria and significantly expands the number of described integron-associated genes. | 2020 | 32689930 |
| 3906 | 15 | 0.9998 | Survival in amoeba--a major selection pressure on the presence of bacterial copper and zinc resistance determinants? Identification of a "copper pathogenicity island". The presence of metal resistance determinants in bacteria usually is attributed to geological or anthropogenic metal contamination in different environments or associated with the use of antimicrobial metals in human healthcare or in agriculture. While this is certainly true, we hypothesize that protozoan predation and macrophage killing are also responsible for selection of copper/zinc resistance genes in bacteria. In this review, we outline evidence supporting this hypothesis, as well as highlight the correlation between metal resistance and pathogenicity in bacteria. In addition, we introduce and characterize the "copper pathogenicity island" identified in Escherichia coli and Salmonella strains isolated from copper- and zinc-fed Danish pigs. | 2015 | 26088177 |
| 4638 | 16 | 0.9998 | Comprehensive 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. | 2021 | 34060909 |
| 4151 | 17 | 0.9998 | Evolutionary relationships among genes for antibiotic resistance. The genes that determine resistance to antibiotics are commonly found encoded by extrachromosomal elements in bacteria. These were described first in Enterobacteriaceae and subsequently in a variety of other genera; their spread is associated with the increased use of antibiotics in human and animal medicine. Antibiotic-resistance genes that determine the production of enzymes which modify (detoxify) the antibiotics have been detected in antibiotic-producing organisms. It has been suggested that the producing strains provided the source of antibiotic-resistance genes that were then 'picked-up' by recombination. Recent studies of the nucleotide sequence of certain antibiotic-resistance genes indicate regions of strong homology in the encoded proteins. The implications of these similarities are discussed. | 1984 | 6559117 |
| 4153 | 18 | 0.9998 | Amino acid variation in the GyrA subunit of bacteria potentially associated with natural resistance to fluoroquinolone antibiotics. In studies of genetic diversity in natural microbial populations, we have analyzed nucleotide sequences in the quinolone resistance-determining region of the bacterial gyrA gene in ciprofloxacin-resistant and nonselected soil bacteria obtained from the environment. It is apparent that this sequence is highly variable, and resistance to fluoroquinolone antibiotics occurring in environmental populations of bacteria is due at least in part to natural sequence variation in this domain. We suggest that the development of new antimicrobial agents, including completely synthetic antimicrobials such as the fluoroquinolones, should incorporate the analysis of resistance mechanisms among microbes in natural environments; these studies could predict potential mechanisms of resistance to be encountered in subsequent clinical use of the agents and would guide chemical modification designed to evade resistance development. | 1997 | 9420056 |
| 9662 | 19 | 0.9998 | Species-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. | 2022 | 35040700 |