Growth-Dependent Predation and Generalized Transduction of Antimicrobial Resistance by Bacteriophage

Bacteriophage (phage) are both predators and evolutionary drivers for bacteria, notably contributing to the spread of antimicrobial resistance (AMR) genes by generalized transduction. Our current understanding of this complex relationship is limited. We used an interdisciplinary approach to quantify how these interacting dynamics can lead to the evolution of multidrug-resistant bacteria. We cocultured two strains of methicillin-resistant Staphylococcus aureus, each harboring a different antibiotic resistance gene, with generalized transducing phage. After a growth phase of 8 h, bacteria and phage surprisingly coexisted at a stable equilibrium in our culture, the level of which was dependent on the starting concentration of phage. We detected double-resistant bacteria as early as 7 h, indicating that transduction of AMR genes had occurred. We developed multiple mathematical models of the bacteria and phage relationship and found that phage-bacteria dynamics were best captured by a model in which phage burst size decreases as the bacteria population reaches stationary phase and where phage predation is frequency-dependent. We estimated that one in every 108 new phage generated was a transducing phage carrying an AMR gene and that double-resistant bacteria were always predominantly generated by transduction rather than by growth. Our results suggest a shift in how we understand and model phage-bacteria dynamics. Although rates of generalized transduction could be interpreted as too rare to be significant, they are sufficient in our system to consistently lead to the evolution of multidrug-resistant bacteria. Currently, the potential of phage to contribute to the growing burden of AMR is likely underestimated

Genomic epidemiology reveals transmission patterns and dynamics of SARS-CoV-2 in Aotearoa New Zealand

New Zealand, a geographically remote Pacific island with easily sealable borders, implemented a nationwide 'lockdown' of all non-essential services to curb the spread of COVID-19. Here, we generate 649 SARS-CoV-2 genome sequences from infected patients in New Zealand with samples collected during the 'first wave', representing 56% of all confirmed cases in this time period. Despite its remoteness, the viruses imported into New Zealand represented nearly all of the genomic diversity sequenced from the global virus population. These data helped to quantify the effectiveness of public health interventions. For example, the effective reproductive number, Re of New Zealand's largest cluster decreased from 7 to 0.2 within the first week of lockdown. Similarly, only 19% of virus introductions into New Zealand resulted in ongoing transmission of more than one additional case. Overall, these results demonstrate the utility of genomic pathogen surveillance to inform public health and disease mitigation

Importance of patient bed pathways and length of stay differences in predicting COVID-19 hospital bed occupancy in England.

Background: Predicting bed occupancy for hospitalised patients with COVID-19 requires understanding of length of stay (LoS) in particular bed types. LoS can vary depending on the patient’s “bed pathway” - the sequence of transfers of individual patients between bed types during a hospital stay. In this study, we characterise these pathways, and their impact on predicted hospital bed occupancy. Methods: We obtained data from University College Hospital (UCH) and the ISARIC4C COVID-19 Clinical Information Network (CO-CIN) on hospitalised patients with COVID-19 who required care in general ward or critical care (CC) beds to determine possible bed pathways and LoS. We developed a discrete-time model to examine the implications of using either bed pathways or only average LoS by bed type to forecast bed occupancy. We compared model-predicted bed occupancy to publicly available bed occupancy data on COVID-19 in England between March and August 2020. Results: In both the UCH and CO-CIN datasets, 82% of hospitalised patients with COVID-19 only received care in general ward beds. We identified four other bed pathways, present in both datasets: “Ward, CC, Ward”, “Ward, CC”, “CC” and “CC, Ward”. Mean LoS varied by bed type, pathway, and dataset, between 1.78 and 13.53 days. For UCH, we found that using bed pathways improved the accuracy of bed occupancy predictions, while only using an average LoS for each bed type underestimated true bed occupancy. However, using the CO-CIN LoS dataset we were not able to replicate past data on bed occupancy in England, suggesting regional LoS heterogeneities. Conclusions: We identified five bed pathways, with substantial variation in LoS by bed type, pathway, and geography. This might be caused by local differences in patient characteristics, clinical care strategies, or resource availability, and suggests that national LoS averages may not be appropriate for local forecasts of bed occupancy for COVID-19. Trial registration: The ISARIC WHO CCP-UK study ISRCTN66726260 was retrospectively registered on 21/04/2020 and designated an Urgent Public Health Research Study by NIHR.</p