Healthcare-associated infections (HAIs) cause mortality, morbidity, and waste of health care resources. HAIs are also an important driver of antimicrobial resistance, which is increasing around the world. Beginning in November 2016, we instituted an initiative to detect outbreaks of HAIs using prospective whole-genome sequencing-based surveillance of bacterial pathogens collected from hospitalized patients. Here, we describe the diversity of bacteria sampled from hospitalized patients at a single center, as revealed through systematic analysis of bacterial isolate genomes. We sequenced the genomes of 3,004 bacterial isolates from hospitalized patients collected over a 25-month period. We identified bacteria belonging to 97 distinct species, which were distributed among 14 groups of related species. Within these groups, isolates could be distinguished from one another by both average nucleotide identity (ANI) and principal-component analysis of accessory genes (PCA-A). Core genome genetic distances and rates of evolution varied among species, which has practical implications for defining shared ancestry during outbreaks and for our broader understanding of the origins of bacterial strains and species. Finally, antimicrobial resistance genes and putative mobile genetic elements were frequently observed, and our systematic analysis revealed patterns of occurrence across the different species sampled from our hospital. Overall, this study shows how understanding the population structure of diverse pathogens circulating in a single health care setting can improve the discriminatory power of genomic epidemiology studies and can help define the processes leading to strain and species differentiation. IMPORTANCE Hospitalized patients are at increased risk of becoming infected with antibiotic-resistant organisms. We used whole-genome sequencing to survey and compare over 3,000 clinical bacterial isolates collected from hospitalized patients at a large medical center over a 2-year period. We identified nearly 100 different bacterial species, which we divided into 14 different groups of related species. When we examined how genetic relatedness differed between species, we found that different species were likely evolving at different rates within our hospital. This is significant because the identification of bacterial outbreaks in the hospital currently relies on genetic similarity cutoffs, which are often applied uniformly across organisms. Finally, we found that antibiotic resistance genes and mobile genetic elements were abundant and were shared among the bacterial isolates we sampled. Overall, this study provides an in-depth view of the genomic diversity and evolutionary processes of bacteria sampled from hospitalized patients, as well as genetic similarity estimates that can inform hospital outbreak detection and prevention efforts.

We performed whole genome sequencing on SARS-CoV-2 from 59 vaccinated individuals from southwest Pennsylvania who tested positive between February and September, 2021. A comparison of mutations among vaccine breakthrough cases to a time-matched control group identified potential adaptive responses of SARS-CoV-2 to vaccination.

The emergence of the early COVID-19 epidemic in the United States (U.S.) went largely undetected, due to a lack of adequate testing and mitigation efforts. The city of New Orleans, Louisiana experienced one of the earliest and fastest accelerating outbreaks, coinciding with the annual Mardi Gras festival, which went ahead without precautions. To gain insight into the emergence of SARS-CoV-2 in the U.S. and how large, crowded events may have accelerated early transmission, we sequenced SARS-CoV-2 genomes during the first wave of the COVID-19 epidemic in Louisiana. We show that SARS-CoV-2 in Louisiana initially had limited sequence diversity compared to other U.S. states, and that one successful introduction of SARS-CoV-2 led to almost all of the early SARS-CoV-2 transmission in Louisiana. By analyzing mobility and genomic data, we show that SARS-CoV-2 was already present in New Orleans before Mardi Gras and that the festival dramatically accelerated transmission, eventually leading to secondary localized COVID-19 epidemics throughout the Southern U.S.. Our study provides an understanding of how superspreading during large-scale events played a key role during the early outbreak in the U.S. and can greatly accelerate COVID-19 epidemics on a local and regional scale.

Background Most hospitals use traditional infection prevention (IP) methods for outbreak detection. We developed the Enhanced Detection System for Healthcare-Associated Transmission (EDS-HAT), which combines whole-genome sequencing (WGS) surveillance and machine learning (ML) of the electronic health record (EHR) to identify undetected outbreaks and the responsible transmission routes, respectively. Methods We performed WGS surveillance of healthcare-associated bacterial pathogens from November 2016 to November 2018. EHR ML was used to identify the transmission routes for WGS-detected outbreaks, which were investigated by an IP expert. Potential infections prevented were estimated and compared with traditional IP practice during the same period. Results Of 3165 isolates, there were 2752 unique patient isolates in 99 clusters involving 297 (10.8%) patient isolates identified by WGS; clusters ranged from 2-14 patients. At least 1 transmission route was detected for 65.7% of clusters. During the same time, traditional IP investigation prompted WGS for 15 suspected outbreaks involving 133 patients, for which transmission events were identified for 5 (3.8%). If EDS-HAT had been running in real time, 25-63 transmissions could have been prevented. EDS-HAT was found to be cost-saving and more effective than traditional IP practice, with overall savings of $192 408-$692 532. Conclusions EDS-HAT detected multiple outbreaks not identified using traditional IP methods, correctly identified the transmission routes for most outbreaks, and would save the hospital substantial costs. Traditional IP practice misidentified outbreaks for which transmission did not occur. WGS surveillance combined with EHR ML has the potential to save costs and enhance patient safety. Whole-genome sequencing surveillance of bacterial pathogens and machine learning of the electronic health record finds previously undetected outbreaks and their transmission routes, which can increase patient safety and save costs.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike protein (S) plays critical roles in host cell entry. Non-synonymous substitutions affecting S are not uncommon and have become fixed in a number of SARS-CoV-2 lineages. A subset of such mutations enable escape from neutralizing antibodies or are thought to enhance transmission through mechanisms such as increased affinity for the cell entry receptor, angiotensin-converting enzyme 2 (ACE2). Independent genomic surveillance programs based in New Mexico and Louisiana contemporaneously detected the rapid rise of numerous clade 20G (lineage B.1.2) infections carrying a Q677P substitution in S. The variant was first detected in the US on October 23, yet between 01 Dec 2020 and 19 Jan 2021 it rose to represent 27.8% and 11.3% of all SARS-CoV-2 genomes sequenced from Louisiana and New Mexico, respectively. Q677P cases have been detected predominantly in the south central and southwest United States; as of 03 Feb 2021, GISAID data show 499 viral sequences of this variant from the USA. Phylogenetic analyses revealed the independent evolution and spread of at least six distinct Q677H sub-lineages, with first collection dates ranging from mid-August to late November 2020. Four 677H clades from clade 20G (B.1.2), 20A (B.1.234), and 20B (B.1.1.220, and B.1.1.222) each contain roughly 100 or fewer sequenced cases, while a distinct pair of clade 20G clusters are represented by 754 and 298 cases, respectively. Although sampling bias and founder effects may have contributed to the rise of S:677 polymorphic variants, the proximity of this position to the polybasic cleavage site at the S1/S2 boundary are consistent with its potential functional relevance during cell entry, suggesting parallel evolution of a trait that may confer an advantage in spread or transmission. Taken together, our findings demonstrate simultaneous convergent evolution, thus providing an impetus to further evaluate S:677 polymorphisms for effects on proteolytic processing, cell tropism, and transmissibility.

As selection frequently favors noncooperating defectors in mixed populations with cooperators, mechanisms that promote cooperation stability clearly exist. One potential mechanism is bacterial cell-to-cell communication, quorum sensing (QS), which can allow cooperators to prevent invasion by defectors. However, the impact of QS on widespread maintenance of cooperation in well-mixed conditions has not been experimentally demonstrated over extended evolutionary timescales. Here, we use wild-type (WT) Vibrio campbellii that regulates cooperation with QS and an unconditional cooperating (UC) mutant to examine the evolutionary origins and dynamics of novel defectors during a long-term evolution experiment. We found that UC lineages were completely outcompeted by defectors, whereas functioning QS enabled the maintenance of cooperative variants in most WT populations. Sequencing evolved populations revealed multiple luxR mutations that swept the UC lineages. However, the evolution of mutant lineages with reduced levels of bioluminescence (dims) occurred in many WT lineages. These dim variants also decreased other cooperative phenotypes regulated by QS, including protease production, indicating they result from changes to QS regulation. This diminished investment phenotype optimizes a tradeoff between cooperative input and growth output and suggests that decreasing the cost of QS could be a favorable strategy for maintaining the cooperative behaviors it regulates.

Many bacteria cycle between sessile and motile forms in which they must sense and respond to internal and external signals to coordinate appropriate physiology. Maintaining fitness requires genetic networks that have been honed in variable environments to integrate these signals. The identity of the major regulators and how their control mechanisms evolved remain largely unknown in most organisms. During four different evolution experiments with the opportunist betaproteobacterium Burkholderia cenocepacia in a biofilm model, mutations were most frequently selected in the conserved gene rpfR. RpfR uniquely integrates two major signaling systems-quorum sensing and the motile-sessile switch mediated by cyclic-di-GMP-by two domains that sense, respond to, and control the synthesis of the auto-inducer cis-2-dodecenoic acid (BDSF). The BDSF response in turn regulates the activity of diguanylate cyclase and phosphodiesterase domains acting on cyclic-di-GMP. Parallel adaptive substitutions evolved in each of these domains to produce unique life history strategies by regulating cyclic-di-GMP levels, global transcriptional responses, biofilm production, and polysaccharide composition. These phenotypes translated into distinct ecology and biofilm structures that enabled mutants to coexist and produce more biomass than expected from their constituents grown alone. This study shows that when bacterial populations are selected in environments challenging the limits of their plasticity, the evolved mutations not only alter genes at the nexus of signaling networks but also reveal the scope of their regulatory functions.

Significance The evolution of new metabolic pathways has been a driver of diversification from the last universal common ancestor 3.8 billion y ago to the present. Bioinformatic evidence suggests that many pathways were assembled by recruiting promiscuous enzymes to serve new functions. However, the processes by which new pathways have emerged are lost in time. We have little information about the environmental conditions that fostered emergence of new pathways, the genome context in which new pathways emerged, and the types of mutations that elevated flux through inefficient new pathways. Experimental laboratory evolution has allowed us to evolve a new pathway and identify mechanisms by which mutations increase fitness when an inefficient new pathway becomes important for survival.

New enzymes often evolve by gene amplification and divergence. Previous experimental studies have followed the evolutionary trajectory of an amplified gene, but have not considered mutations elsewhere in the genome when fitness is limited by an evolving gene. We have evolved a strain of Escherichia coli in which a secondary promiscuous activity has been recruited to serve an essential function. The gene encoding the ‘weak-link’ enzyme amplified in all eight populations, but mutations improving the newly needed activity occurred in only one. Most adaptive mutations occurred elsewhere in the genome. Some mutations increase expression of the enzyme upstream of the weak-link enzyme, pushing material through the dysfunctional metabolic pathway. Others enhance production of a co-substrate for a downstream enzyme, thereby pulling material through the pathway. Most of these latter mutations are detrimental in wild-type E. coli, and thus would require reversion or compensation once a sufficient new activity has evolved.

Bacterial populations vary in their stress tolerance and population structure depending upon whether growth occurs in well-mixed or structured environments. We hypothesized that evolution in biofilms would generate greater genetic diversity than well-mixed environments and lead to different pathways of antibiotic resistance. We used experimental evolution and whole genome sequencing to test how the biofilm lifestyle influenced the rate, genetic mechanisms, and pleiotropic effects of resistance to ciprofloxacin in Acinetobacter baumannii populations. Both evolutionary dynamics and the identities of mutations differed between lifestyle. Planktonic populations experienced selective sweeps of mutations including the primary topoisomerase drug targets, whereas biofilm-adapted populations acquired mutations in regulators of efflux pumps. An overall trade-off between fitness and resistance level emerged, wherein biofilm-adapted clones were less resistant than planktonic but more fit in the absence of drug. However, biofilm populations developed collateral sensitivity to cephalosporins, demonstrating the clinical relevance of lifestyle on the evolution of resistance.

The diffusible signal factors (DSFs) are a family of quorum-sensing autoinducers (AIs) produced and detected by numerous gram-negative bacteria. The DSF family AIs are fatty acids, differing in their acyl chain length, branching, and substitution but having in common a cis-2 double bond that is required for their activity. In both human and plant pathogens, DSFs regulate diverse phenotypes, including virulence factor expression, antibiotic resistance, and biofilm dispersal. Despite their widespread relevance to both human health and agriculture, the molecular basis of DSF recognition by their cellular receptors remained a mystery. Here, we report the first structure-function studies of the DSF receptor regulation of pathogenicity factor R (RpfR). We present the X-ray crystal structure of the RpfR DSF-binding domain in complex with the Burkholderia DSF (BDSF), which to our knowledge is the first structure of a DSF receptor in complex with its AI. To begin to understand the mechanistic role of the BDSF-RpfR contacts observed in the biologically important complex, we have also determined the X-ray crystal structure of the RpfR DSF-binding domain in complex with the inactive, saturated isomer of BDSF, dodecanoic acid (C12:0). In addition to these ligand-receptor complex structures, we report the discovery of a previously overlooked RpfR domain and show that it binds to and negatively regulates the DSF synthase regulation of pathogenicity factor F (RpfF). We have named this RpfR region the RpfF interaction (FI) domain, and we have determined its X-ray crystal structure alone and in complex with RpfF. These X-ray crystal structures, together with extensive complementary in vivo and in vitro functional studies, reveal the molecular basis of DSF recognition and the importance of the cis-2 double bond to DSF function. Finally, we show that throughout cellular growth, the production of BDSF by RpfF is post-translationally controlled by the RpfR N-terminal FI domain, affecting the cellular concentration of the bacterial second messenger bis-(3'-5')-cyclic dimeric guanosine monophosphate (c-di-GMP). Thus, in addition to describing the molecular basis for the binding and specificity of a DSF for its receptor, we describe a receptor-synthase interaction regulating bacterial quorum-sensing signaling and second messenger signal transduction.

Unregulated protein-tyrosine kinase signaling is a common feature of AML, often involving mutations in Flt3 and overexpression of myeloid Src-family kinases (Hck, Fgr, Lyn). Here we show that high-level expression of these Src kinases predicts poor survival in a large cohort of AML patients. To test the therapeutic benefit of Flt3 and Src-family kinase inhibition, we used the pyrrolopyrimidine kinase inhibitor A-419259. This compound potently inhibits Hck, Fgr, and Lyn as well as Flt3 bearing an activating internal tandem duplication (ITD). Flt3-ITD expression sensitized human TF-1 myeloid cells to growth arrest by A-419259, supporting direct action on the Flt3-ITD kinase domain. Cells transformed with the Flt3-ITD mutants D835Y and F691L were resistant to A-419259, while co-expression of Hck or Fgr restored inhibitor sensitivity to Flt3-ITD D835Y. Conversely, Hck and Fgr mutants with engineered A-419259 resistance mutations decreased sensitivity of TF-1/Flt3-ITD cells. To investigate de novo resistance mechanisms, A-419259-resistant Flt3-ITD+ AML cell populations were derived via long-term dose escalation. Whole exome sequencing identified a distinct Flt3-ITD kinase domain mutation (N676S/T) among all A-419259 target kinases in each of six independent resistant cell populations. These studies show that Hck and Fgr expression influences inhibitor sensitivity and the pathway to acquired resistance in Flt3-ITD+ AML.

Bacteria adapt to infections by evolving variants that are more fit and persistent. These recalcitrant variants are typically observed in chronic infections. However, it is unclear when and why these variants evolve. To address these questions, we used a porcine chronic wound model to study the evolutionary dynamics of Pseudomonas aeruginosa in a mixed-strain infection. We isolated hyperbiofilm variants that persisted early in the infection. Interstrain interactions were also observed, where adapted variants acquired CRISPR-mediated immunity to phages. We show that when initiating infection, P. aeruginosa experiences strong positive selection for hyperbiofilm phenotypes produced by mutants of a single chemosensory system, the Wsp pathway. We predict that hyperbiofilm variants are early adaptations to infection and that interstrain interactions may influence bacterial burden and infection outcomes.

Synonymous mutations do not alter the specified amino acid but may alter the structure or function of an mRNA in ways that impact fitness. There are few examples in the literature, however, in which the effects of synonymous mutations on microbial growth rates have been measured, and even fewer for which the underlying mechanism is understood. We evolved four populations of a strain of Salmonella enterica in which a promiscuous enzyme has been recruited to replace an essential enzyme. A previously identified point mutation increases the enzyme’s ability to catalyze the newly needed reaction (required for arginine biosynthesis) but decreases its ability to catalyze its native reaction (required for proline biosynthesis). The poor performance of this enzyme limits growth rate on glucose. After 260 generations, we identified two synonymous mutations in the first six codons of the gene encoding the weak-link enzyme that increase growth rate by 41 and 67%. We introduced all possible synonymous mutations into the first six codons and found substantial effects on growth rate; one doubles growth rate, and another completely abolishes growth. Computational analyses suggest that these mutations affect either the stability of a stem-loop structure that sequesters the start codon or the accessibility of the region between the Shine-Dalgarno sequence and the start codon. Thus, these mutations would be predicted to affect translational efficiency and thereby indirectly affect mRNA stability because translating ribosomes protect mRNA from degradation. Experimental data support these hypotheses. We conclude that the effects of the synonymous mutations are due to a combination of effects on mRNA stability and translation efficiency that alter levels of the weak-link enzyme. These findings suggest that synonymous mutations can have profound effects on fitness under strong selection and that their importance in evolution may be under-appreciated.

ESG and CWM performed the experimental work. JSH performed the colony PCR. SBC infected and sampled the porcine burn wounds. DS generated the sequence library. MW and GRL quantified the strain frequency in the wounds. ESG, CWM, MW, VSC, and DJW conceptualized the project and wrote the manuscript.

Background. Data on the use of ceftolozane-tazobactam and emergence of ceftolozane-tazobactam resistance during multidrug resistant (MDR)-Pseudomonas aeruginosa infections are limited.Methods. We performed a retrospective study of 21 patients treated with ceftolozane-tazobactam for MDR-P. aeruginosa infections. Whole genome sequencing and quantitative real-time polymerase chain reaction were performed on longitudinal isolates.Results. Median age was 58 years; 9 patients (43%) were transplant recipients. Median simplified acute physiology score-II (SAPS-II) was 26. Eighteen (86%) patients were treated for respiratory tract infections; others were treated for bloodstream, complicated intraabdominal infections, or complicated urinary tract infections. Ceftolozane-tazobactam was discontinued in 1 patient (rash). Thirty-day all-cause and attributable mortality rates were 10% (2/21) and 5% (1/21), respectively; corresponding 90-day mortality rates were 48% (10/21) and 19% (4/21). The ceftolozane-tazobactam failure rate was 29% (6/21). SAPS-II score was the sole predictor of failure. Ceftolozane-tazobactam resistance emerged in 3 (14%) patients. Resistance was associated with de novo mutations, rather than acquisition of resistant nosocomial isolates. ampC overexpression and mutations were identified as potential resistance determinants.Conclusions. In this small study, ceftolozane-tazobactam was successful in treating 71% of patients with MDR-P. aeruginosa infections, most of whom had pneumonia. The emergence of ceftolozane-tazobactam resistance in 3 patients is worrisome and may be mediated in part by AmpC-related mechanisms. More research on treatment responses and resistance during various types of MDR-P. aeruginosa infections is needed to define ceftolozane-tazobactam's place in the armamentarium.