Development of an Amplicon Nanopore Sequencing Strategy for Detection of Mutations Conferring Intermediate Resistance to Vancomycin in Staphylococcus aureus Strains.

Staphylococcus aureus is a major cause of bacteremia and other hospital-acquired infections. The cell-wall active antibiotic vancomycin is commonly used to treat both methicillin-resistant (MRSA) and sensitive (MSSA) infections. Vancomycin intermediate S. aureus (VISA) variants can arise through de novo mutations. Here, we performed pilot experiments to develop a combined PCR/long-read sequencing-based method for detection of previously known VISA-causing mutations. Primers were designed to generate 10 amplicons covering 16 genes associated with the VISA phenotype. We sequenced amplicon pools as long reads with Oxford Nanopore adapter ligation on Flongle flow cells. We then detected mutations by mapping reads against a parental consensus or known reference sequence and comparing called variants against a database of known VISA mutations from laboratory selection. Each amplicon in the pool was sequenced to high (>1,000×) coverage, and no relationship was found between amplicon length and coverage. We also were able to detect the causative mutation (walK 646C>G) in a VISA mutant derived from the USA300 strain (N384-3 from parental strain N384). Mixing mutant (N384-3) and parental (N384) DNA at various ratios from 0 to 1 mutant suggested a mutation detection threshold of the average minor allele frequency (6.5%) at 95% confidence (two standard errors above mean mutation frequency). The study lays the groundwork for direct S. aureus antibiotic resistance genotype inference using rapid nanopore sequencing from clinical samples. IMPORTANCE Bacteremia mortality is known to increase rapidly with time after infection, making rapid diagnostics and treatment necessary. Successful treatment depends on correct administration of antibiotics based on knowledge of strain antibiotic susceptibility. Staphylococcus aureus is a major causative agent of bacteremia that is also commonly antibiotic resistant. In this work, we develop a method to accelerate detection of a complex, polygenic antibiotic resistance phenotype in S. aureus, vancomycin-intermediate resistance (VISA), through long-read genomic sequencing of amplicons representing genes most commonly mutated in VISA selection. This method both rapidly identifies VISA genotypes and incorporates the most comprehensive database of VISA genetic determinants known to date.

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.

Symbiont Genomic Features and Localization in the Bean Beetle Callosobruchus maculatus.

A pervasive pest of stored leguminous products, the bean beetle Callosobruchus maculatus (Coleoptera: Chrysomelidae) associates with a simple bacterial community during adulthood. Despite its economic importance, little is known about the compositional stability, heritability, localization, and metabolic potential of the bacterial symbionts of C. maculatus. In this study, we applied community profiling using 16S rRNA gene sequencing to reveal a highly conserved bacterial assembly shared between larvae and adults. Dominated by Firmicutes and Proteobacteria, this community is localized extracellularly along the epithelial lining of the bean beetle's digestive tract. Our analysis revealed that only one species, Staphylococcus gallinarum (phylum Firmicutes), is shared across all developmental stages. Isolation and whole-genome sequencing of S. gallinarum from the beetle gut yielded a circular chromosome (2.8 Mb) and one plasmid (45 kb). The strain encodes complete biosynthetic pathways for the production of B vitamins and amino acids, including tyrosine, which is increasingly recognized as an important symbiont-supplemented precursor for cuticle biosynthesis in beetles. A carbohydrate-active enzyme search revealed that the genome codes for a number of digestive enzymes, reflecting the nutritional ecology of C. maculatus. The ontogenic conservation of the gut microbiota in the bean beetle, featuring a "core" community composed of S. gallinarum, may be indicative of an adaptive role for the host. In clarifying symbiont localization and metabolic potential, we further our understanding and study of a costly pest of stored products. IMPORTANCE From supplementing essential nutrients to detoxifying plant secondary metabolites and insecticides, bacterial symbionts are a key source of adaptations for herbivorous insect pests. Despite the pervasiveness and geographical range of the bean beetle Callosobruchus maculatus, the role of microbial symbioses in its natural history remains understudied. Here, we demonstrate that the bean beetle harbors a simple gut bacterial community that is stable throughout development. This community localizes along the insect's digestive tract and is largely dominated by Staphylococcus gallinarum. In elucidating symbiont metabolic potential, we highlight its possible adaptive significance for a widespread agricultural pest.

Genes Influencing Phage Host Range in Staphylococcus aureus on a Species-Wide Scale.

Staphylococcus aureus is a human pathogen that causes serious diseases, ranging from skin infections to septic shock. Bacteriophages (phages) are both natural killers of S. aureus, offering therapeutic possibilities, and important vectors of horizontal gene transfer (HGT) in the species. Here, we used high-throughput approaches to understand the genetic basis of strain-to-strain variation in sensitivity to phages, which defines the host range. We screened 259 diverse S. aureus strains covering more than 40 sequence types for sensitivity to eight phages, which were representatives of the three phage classes that infect the species. The phages were variable in host range, each infecting between 73 and 257 strains. Using genome-wide association approaches, we identified putative loci that affect host range and validated their function using USA300 transposon knockouts. In addition to rediscovering known host range determinants, we found several previously unreported genes affecting bacterial growth during phage infection, including trpA, phoR, isdB, sodM, fmtC, and relA We used the data from our host range matrix to develop predictive models that achieved between 40% and 95% accuracy. This work illustrates the complexity of the genetic basis for phage susceptibility in S. aureus but also shows that with more data, we may be able to understand much of the variation. With a knowledge of host range determination, we can rationally design phage therapy cocktails that target the broadest host range of S. aureus strains and address basic questions regarding phage-host interactions, such as the impact of phage on S. aureus evolution.IMPORTANCEStaphylococcus aureus is a widespread, hospital- and community-acquired pathogen, many strains of which are antibiotic resistant. It causes diverse diseases, ranging from local to systemic infection, and affects both the skin and many internal organs, including the heart, lungs, bones, and brain. Its ubiquity, antibiotic resistance, and disease burden make new therapies urgent. One alternative therapy to antibiotics is phage therapy, in which viruses specific to infecting bacteria clear infection. In this work, we identified and validated S. aureus genes that influence phage host range-the number of strains a phage can infect and kill-by testing strains representative of the diversity of the S. aureus species for phage host range and associating the genome sequences of strains with host range. These findings together improved our understanding of how phage therapy works in the bacterium and improve prediction of phage therapy efficacy based on the predicted host range of the infecting strain.

Gas gangrene-associated gliding motility is regulated by the Clostridium perfringens CpAL/VirSR system.

Clostridium perfringens strains cause a wide variety of human and animal disease, including gas gangrene or myonecrosis. Production of toxins required for myonecrosis, PFO and CPA, is regulated by the C. perfringens Agr-like (CpAL) system via the VirSR two-component system. Myonecrosis begins at the site of infection from where bacteria migrate deep into the host tissue likely using a previously described gliding motility phenotype. We therefore assessed whether gliding motility was under the control of the CpAL/VirSR regulon. The migration rate of myonecrosis-causing C. perfringens strain 13 (S13) was investigated during a 96 h period, including an adaptation phase with bacterial migration (∼1.4 mm/day) followed by a gliding phase allowing bacteria faster migration (∼8.6 mm/day). Gliding required both an intact CpAL system, and signaling through VirSR. Mutants lacking ΔagrB, or ΔvirR, were impaired for onward gliding while a complemented strain S13ΔagrB/pTS1303 had the gliding phenotype restored. Gene expression studies revealed upregulated transcription of pili genes (pilA1, pilA2 and pilT) whose encoded proteins were previously found to be required for gliding motility and CpAL/VirSR-regulated pfoA and cpa toxin genes. Compared to S13, transcription of cpa and pfoA significantly decreased in S13ΔagrB, or S13ΔvirR, strains but not that of pili genes. Further experiments demonstrated that mutants S13ΔpfoA and S13Δcpa migrated at the same rate as S13 wt. We demonstrated that CpAL/VirSR regulates C. perfringens gliding motility and that gliding bacteria have an increased transcription of toxin genes involved in myonecrosis.

Interaction between Streptococcus pneumoniae and Staphylococcus aureus Generates ·OH Radicals That Rapidly Kill Staphylococcus aureus Strains.

Streptococcus pneumoniae rapidly kills Staphylococcus aureus by producing membrane-permeable hydrogen peroxide (H2O2). The mechanism by which S. pneumoniae-produced H2O2 mediates S. aureus killing was investigated. An in vitro model that mimicked S. pneumoniae-S. aureus contact during colonization of the nasopharynx demonstrated that S. aureus killing required outcompeting densities of S. pneumoniae Compared to the wild-type strain, isogenic S. pneumoniae ΔlctO and S. pneumoniae ΔspxB, both deficient in production of H2O2, required increased density to kill S. aureus While residual H2O2 activity produced by single mutants was sufficient to eradicate S. aureus, an S. pneumoniae ΔspxB ΔlctO double mutant was unable to kill S. aureus A collection of 20 diverse methicillin-resistant S. aureus (MRSA) and methicillin-susceptible S. aureus (MSSA) strains showed linear sensitivity (R2 = 0.95) for S. pneumoniae killing, but the same strains had different susceptibilities when challenged with pure H2O2 (5 mM). There was no association between the S. aureus clonal complex and sensitivity to either S. pneumoniae or H2O2 To kill S. aureus, S. pneumoniae produced ∼180 µM H2O2 within 4 h of incubation, while the killing-defective S. pneumoniae ΔspxB and S. pneumoniae ΔspxB ΔlctO mutants produced undetectable levels. Remarkably, a sublethal dose (1 mM) of pure H2O2 incubated with S. pneumoniae ΔspxB eradicated diverse S. aureus strains, suggesting that S. pneumoniae bacteria may facilitate conversion of H2O2 to a hydroxyl radical (·OH). Accordingly, S. aureus killing was completely blocked by incubation with scavengers of ·OH radicals, dimethyl sulfoxide (Me2SO), thiourea, or sodium salicylate. The ·OH was detected in S. pneumoniae cells by spin trapping and electron paramagnetic resonance. Therefore, S. pneumoniae produces H2O2, which is rapidly converted to a more potent oxidant, hydroxyl radicals, to rapidly intoxicate S. aureus strains.IMPORTANCEStreptococcus pneumoniae strains produce hydrogen peroxide (H2O2) to kill bacteria in the upper airways, including pathogenic Staphylococcus aureus strains. The targets of S. pneumoniae-produced H2O2 have not been discovered, in part because of a lack of knowledge about the underlying molecular mechanism. We demonstrated that an increased density of S. pneumoniae kills S. aureus by means of H2O2 produced by two enzymes, SpxB and LctO. We discovered that SpxB/LctO-produced H2O2 is converted into a hydroxyl radical (·OH) that rapidly intoxicates and kills S. aureus We successfully inhibited the toxicity of ·OH with three different scavengers and detected ·OH in the supernatant. The target(s) of the hydroxyl radicals represents a new alternative for the development of antimicrobials against S. aureus infections.

Determinants of Phage Host Range in Staphylococcus Species.

Bacteria in the genus Staphylococcus are important targets for phage therapy due to their prevalence as pathogens and increasing antibiotic resistance. Here we review Staphylococcus outer surface features and specific phage resistance mechanisms that define the host range, the set of strains that an individual phage can potentially infect. Phage infection goes through five distinct phases: attachment, uptake, biosynthesis, assembly, and lysis. Adsorption inhibition, encompassing outer surface teichoic acid receptor alteration, elimination, or occlusion, limits successful phage attachment and entry. Restriction-modification systems (in particular, type I and IV systems), which target phage DNA inside the cell, serve as the major barriers to biosynthesis as well as transduction and horizontal gene transfer between clonal complexes and species. Resistance to late stages of infection occurs through mechanisms such as assembly interference, in which staphylococcal pathogenicity islands siphon away superinfecting phage proteins to package their own DNA. While genes responsible for teichoic acid biosynthesis, capsule, and restriction-modification are found in most Staphylococcus strains, a variety of other host range determinants (e.g., clustered regularly interspaced short palindromic repeats, abortive infection, and superinfection immunity) are sporadic. The fitness costs of phage resistance through teichoic acid structure alteration could make staphylococcal phage therapies promising, but host range prediction is complex because of the large number of genes involved, and the roles of many of these are unknown. In addition, little is known about the genetic determinants that contribute to host range expansion in the phages themselves. Future research must identify host range determinants, characterize resistance development during infection and treatment, and examine population-wide genetic background effects on resistance selection.

Whole-Genome Sequences of Staphylococcus aureus Isolates from Cystic Fibrosis Lung Infections.

Staphylococcus aureus is an early colonizer in the lungs of individuals with cystic fibrosis (CF), but surprisingly, only a limited number of genomes from CF-associated S. aureus isolates have been sequenced. Here, we present the whole-genome sequences of 65 S. aureus isolates obtained from 50 individuals with CF.

MetaCRAST: reference-guided extraction of CRISPR spacers from unassembled metagenomes.

Clustered regularly interspaced short palindromic repeat (CRISPR) systems are the adaptive immune systems of bacteria and archaea against viral infection. While CRISPRs have been exploited as a tool for genetic engineering, their spacer sequences can also provide valuable insights into microbial ecology by linking environmental viruses to their microbial hosts. Despite this importance, metagenomic CRISPR detection remains a major challenge. Here we present a reference-guided CRISPR spacer detection tool (Metagenomic CRISPR Reference-Aided Search Tool-MetaCRAST) that constrains searches based on user-specified direct repeats (DRs). These DRs could be expected from assembly or taxonomic profiles of metagenomes. We compared the performance of MetaCRAST to those of two existing metagenomic CRISPR detection tools-Crass and MinCED-using both real and simulated acid mine drainage (AMD) and enhanced biological phosphorus removal (EBPR) metagenomes. Our evaluation shows MetaCRAST improves CRISPR spacer detection in real metagenomes compared to the de novo CRISPR detection methods Crass and MinCED. Evaluation on simulated metagenomes show it performs better than de novo tools for Illumina metagenomes and comparably for 454 metagenomes. It also has comparable performance dependence on read length and community composition, run time, and accuracy to these tools. MetaCRAST is implemented in Perl, parallelizable through the Many Core Engine (MCE), and takes metagenomic sequence reads and direct repeat queries (FASTA or FASTQ) as input. It is freely available for download at https://github.com/molleraj/MetaCRAST.

Determining virus-host interactions and glycerol metabolism profiles in geographically diverse solar salterns with metagenomics.

Solar salterns are excellent model ecosystems for studying virus-microbial interactions because of their low microbial diversity, environmental stability, and high viral density. By using the power of CRISPR spacers to link viruses to their prokaryotic hosts, we explored virus-host interactions in geographically diverse salterns. Using taxonomic profiling, we identified hosts such as archaeal Haloquadratum, Halorubrum, and Haloarcula and bacterial Salinibacter, and we found that community composition related to not only salinity but also local environmental dynamics. Characterizing glycerol metabolism genes in these metagenomes suggested Halorubrum and Haloquadratum possess most dihydroxyacetone kinase genes while Salinibacter possesses most glycerol-3-phosphate dehydrogenase genes. Using two different methods, we detected fewer CRISPR spacers in Haloquadratum-dominated compared with Halobacteriaceae-dominated saltern metagenomes. After CRISPR detection, spacers were aligned against haloviral genomes to map virus to host. While most alignments for each saltern metagenome linked viruses to Haloquadratum walsbyi, there were also alignments indicating interactions with the low abundance taxa Haloarcula and Haloferax. Further examination of the dinucleotide and trinucleotide usage differences between paired viruses and their hosts confirmed viruses and hosts had similar nucleotide usage signatures. Detection of cas genes in the salterns supported the possibility of CRISPR activity. Taken together, our studies suggest similar virus-host interactions exist in different solar salterns and that the glycerol metabolism gene dihydroxyacetone kinase is associated with Haloquadratum and Halorubrum.

Biochemical characterization and zinc binding group (ZBGs) inhibition studies on the catalytic domain of MMP7 (cdMMP7).

Matrix metalloproteinase 7 (MMP7/matrilysin-1) has been implicated in many pathological conditions, such as in cancer and inflammatory diseases; therefore, MMP7 has been targeted for drugs. Success in developing a clinical inhibitor, which exhibits suitable specificity and selectivity, will likely require structural and/or kinetic evaluation of enzyme/inhibitor interactions. To enable these future studies we herein describe the over-expression, purification, and characterization of the catalytic domain of MMP7 (cdMMP7). cdMMP7 was over-expressed in an E. coli over-expression system, and the resulting enzyme was processed into inclusion bodies, which were subsequently solubilized, enabling the enzyme to be re-folded into a catalytically-active form. cdMMP7 was shown to bind 1.8eq of Zn(II), exhibit steady-state kinetic constants of 0.4s-1 for kcat and 23µM for Km, and yield CD and fluorescence spectra that are consistent with a properly-folded enzyme. Pre-steady state kinetic studies yielded kinetic mechanisms of cdMMP7, and these mechanisms are similar to those of other MMPs. Inhibition studies on cdMMP7 with four zinc binding group (ZBG) inhibitors showed that maltol, thiomaltol, and allothiomaltol are better inhibitors with lower IC50 values and lower Kd values against cdMMP7 and cdMMP16 than the commonly-used ZBG inhibitor acetohydroxamic acid. Docking studies suggest that improved inhibitory character may be due to interactions with the S1' substrate binding pocket. Finally, a ZnCo-heterobimetallic analog of cdMMP7 with Co(II) bound in the catalytic site was prepared and characterized. This study describes a well-characterized analog of MMP7 that is available for future inhibitor design efforts.

Investigating the position of the hairpin loop in New Delhi metallo-ß-lactamase, NDM-1, during catalysis and inhibitor binding.

In an effort to examine the relative position of a hairpin loop in New Delhi metallo-ß-lactamase, NDM-1, during catalysis, rapid freeze quench double electron electron resonance (RFQ-DEER) spectroscopy was used. A doubly-labeled mutant of NDM-1, which had one spin label on the invariant loop at position 69 and another label at position 235, was prepared and characterized. The reaction of the doubly spin labeled mutant with chromacef was freeze quenched at 500µs and 10ms. DEER results showed that the average distance between labels decreased by 4Å in the 500µs quenched sample and by 2Å in the 10ms quenched sample, as compared to the distance in the unreacted enzyme, although the peaks corresponding to distance distributions were very broad. DEER spectra with the doubly spin labeled enzyme with two inhibitors showed that the distance between the loop residue at position 69 and the spin label at position 235 does not change upon inhibitor binding. This study suggests that the hairpin loop in NDM-1 moves over the metal ion during the catalysis and then moves back to its original position after hydrolysis, which is consistent with a previous hypothesis based on NMR solution studies on a related metallo-ß-lactamase. This study also demonstrates that this loop motion occurs in the millisecond time domain.

Biochemical, mechanistic, and spectroscopic characterization of metallo-ß-lactamase VIM-2.

This study examines metal binding to metallo-ß-lactamase VIM-2, demonstrating the first successful preparation of a Co(II)-substituted VIM-2 analogue. Spectroscopic studies of the half- and fully metal loaded enzymes show that both Zn(II) and Co(II) bind cooperatively, where the major species present, regardless of stoichiometry, are apo- and di-Zn (or di-Co) enzymes. We determined the di-Zn VIM-2 structure to a resolution of 1.55 Å, and this structure supports results from spectroscopic studies. Kinetics, both steady-state and pre-steady-state, show that VIM-2 utilizes a mechanism that proceeds through a very short-lived anionic intermediate when chromacef is used as the substrate. Comparison with other B1 enzymes shows that those that bind Zn(II) cooperatively are better poised to protonate the intermediate on its formation, compared to those that bind Zn(II) non-cooperatively, which uniformly build up substantial amounts of the intermediate.

Structural and kinetic studies on metallo-ß-lactamase IMP-1.

In an effort to probe for metal binding to metallo-ß-lactamase (MßL) IMP-1, the enzyme was overexpressed, purified, and characterized. The resulting enzyme was shown to bind 2 equiv of Zn(II), exhibit significant catalytic activity, and yield EXAFS results similar to crystallographic data previously reported. Rapid kinetic studies showed that IMP-1 does not stabilize a nitrocefin-derived reaction intermediate; rather, the enzyme follows a simple Michaelis mechanism to hydrolyze nitrocefin. Metal-substituted and metal-reconstituted analogues of IMP-1 were prepared by directly adding metal ion stocks to metal-free enzyme, which was generated by dialysis versus EDTA. UV-vis studies on IMP-1 containing 1 equiv of Co(II) showed a strong ligand-to-metal charge transition at 340 nm, and the intensity of this feature increased when the second equivalent of Co(II) was added to the enzyme. EXAFS fits on IMP-1 containing 1 equiv of Co(II) strongly suggest the presence of a metal-metal interaction, and EPR spectra of the IMP-1 containing 1 and 2 equiv of Co(II) are very similar. Taken together, steady-state kinetic and spectroscopic studies suggest that metal binding to metal-free IMP-1 follows a positive-cooperative mode.