Threshold-free genomic cluster detection to track transmission pathways in health-care settings: a genomic epidemiology analysis

Machine learning reveals the dynamic importance of accessory sequences for Salmonella outbreak clustering

Bacterial typing at whole-genome scales is now feasible owing to decreasing costs in high-throughput sequencing and the recent advances in computation. The unprecedented resolution of whole-genome typing is achieved by genotyping the variable segments of bacterial genomes that can fluctuate significantly in gene content. However, due to the transient and hypervariable nature of many accessory elements, the value of the added resolution in outbreak investigations remains disputed. To assess the analytical value of bacterial accessory genomes in clustering epidemiologically related cases, we trained classifiers on a set of genomes collected from 24 Salmonella enterica outbreaks of food, animal, or environmental origin. The models demonstrated high precision and recall on unseen test data with near-perfect accuracy in classifying clonal and short-term outbreaks. Annotating the genomic features important for cluster classification revealed functional enrichment of molecular fingerprints in genes involved in membrane transportation, trafficking, and carbohydrate metabolism. Importantly, we discovered polymorphisms in mobile genetic elements (MGEs) and gain/loss of MGEs to be informative in defining outbreak clusters. To quantify the ability of MGE variations to cluster outbreak clones, we devised a reference-free tree-building algorithm inspired by colored de Bruijn graphs, which enabled topological comparisons between MGE and standard typing methods. Systematic evaluation of clustering MGEs on an unseen dataset of 34 Salmonella outbreaks yielded mixed results that exemplified the power of accessory sequence variations when core genomes of unrelated cases are insufficiently discriminatory, as well as the distortion of outbreak signals by microevolution events or the incomplete assembly of MGEs.

IMPORTANCE

Gene-by-gene typing is widely used to detect clusters of foodborne illnesses that share a common origin. It remains actively debated whether the inclusion of accessory sequences in bacterial typing schema is informative or deleterious for cluster definitions in outbreak investigations due to the potential confounding effects of horizontal gene transfer. By training machine learning models on a curated set of historical Salmonella outbreaks, we revealed an enriched presence of outbreak distinguishing features in a wide range of mobile genetic elements. Systematic comparison of the efficacy of clustering different accessory elements against standard sequence typing methods led to our cataloging of scenarios where accessory sequence variations were beneficial and uninformative to resolving outbreak clusters. The presented work underscores the complexity of the molecular trends in enteric outbreaks and seeks to inspire novel computational ways to exploit whole-genome sequencing data in enteric disease surveillance and management.

Evidence of Lactobacillus strains shared between the female urinary and vaginal microbiota

Lactobacillus species are common inhabitants of the 'healthy' female urinary and vaginal communities, often associated with a lack of symptoms in both anatomical sites. Given identification by prior studies of similar bacterial species in both communities, it has been hypothesized that the two microbiotas are in fact connected. Here, we carried out whole-genome sequencing of 49 Lactobacillus strains, including 16 paired urogenital samples from the same participant. These strains represent five different Lactobacillus species: L. crispatus, L. gasseri, L. iners, L. jensenii, and L. paragasseri. Average nucleotide identity (ANI), alignment, single-nucleotide polymorphism (SNP), and CRISPR comparisons between strains from the same participant were performed. We conducted simulations of genome assemblies and ANI comparisons and present a statistical method to distinguish between unrelated, related, and identical strains. We found that 50 % of the paired samples have identical strains, evidence that the urinary and vaginal communities are connected. Additionally, we found evidence of strains sharing a common ancestor. These results establish that microbial sharing between the urinary tract and vagina is not limited to uropathogens. Knowledge that these two anatomical sites can share lactobacilli in females can inform future clinical approaches.

Diversity and Persistence of MRSA and VRE in Skilled Nursing Facilities: Environmental Screening, Whole Genome Sequencing, Development of a Dispersion Index

BACKGROUND

Environmental contamination with methicillin-resistant Staphylococcus aureus (MRSA) and vancomycin-resistant Enterococcus (VRE) in skilled nursing facilities (SNFs) may contribute to patient acquisition. We assessed diversity and association of MRSA and VRE isolates in a SNF wing and developed a mathematical index to define each strain's tendency to persist in rooms and spread horizontally.

METHODS

Longitudinal study of MRSA and VRE colonization and contamination among successive patient occupancies in a cluster of nine SNF private rooms during eight months characterized by microbiological testing and whole genome isolate typing. 'Dispersion index" of a strain is defined as the number of rooms it was found in (including the patient), divided by the average of times it was found consecutively in the same room.

FINDINGS

MRSA (ten strain types) and VRE (seven types) were recovered from room or patient in 16.4% and 35.6% of the occupancies, respectively. MRSA showed moderate horizontal spread and several episodes of same-room persistence (three distinct strain types) (overall dispersion index: 1.08). VRE showed high tendency towards horizontal spread /new introductions (overall dispersion index: 3.25), and only one confirmed persistence episode.

INTERPRETATION

The emerging picture of high diversity among contaminating strains and high likelihood of room persistence despite terminal cleaning (MRSA) and horizontal spread between rooms (VRE) in this setting calls for improved cleaning practices, heightened contact precautions, and most of all to establish individually tailored facility screening programs to enable informed choices based on local, measurable and actionable epidemiologic parameters.

FUNDING

University of Michigan OAIC REC Scholarship to M.C. National Institutes of Health K24 AG050685 to L.M.

A point mutation in recC associated with subclonal replacement of carbapenem-resistant Klebsiella pneumoniae ST11 in China

Adaptation to selective pressures is crucial for clinically important pathogens to establish epidemics, but the underlying evolutionary drivers remain poorly understood. The current epidemic of carbapenem-resistant Klebsiella pneumoniae (CRKP) poses a significant threat to public health. In this study we analyzed the genome sequences of 794 CRKP bloodstream isolates collected in 40 hospitals in China between 2014 and 2019. We uncovered a subclonal replacement in the predominant clone ST11, where the previously prevalent subclone OL101:KL47 was replaced by O2v1:KL64 over time in a stepwise manner. O2v1:KL64 carried a higher load of mobile genetic elements, and a point mutation exclusively detected in the recC of O2v1:KL64 significantly promotes recombination proficiency. The epidemic success of O2v1:KL64 was further associated with a hypervirulent sublineage with enhanced resistance to phagocytosis, sulfamethoxazole-trimethoprim, and tetracycline. The phenotypic alterations were linked to the overrepresentation of hypervirulence determinants and antibiotic genes conferred by the acquisition of an rmpA-positive pLVPK-like virulence plasmid and an IncFII-type multidrug-resistant plasmid, respectively. The dissemination of the sublineage was further promoted by more frequent inter-hospital transmission. The results collectively demonstrate that the expansion of O2v1:KL64 is correlated to a repertoire of genomic alterations convergent in a subpopulation with evolutionary advantages.

Transmission of Carbapenem-Resistant Klebsiella pneumoniae in US Hospitals

BACKGROUND

Carbapenem-resistant Klebsiella pneumoniae (CRKp) is the most prevalent carbapenem-resistant Enterobacterales in the United States. We evaluated CRKp clustering in patients in US hospitals.

METHODS

From April 2016 to August 2017, 350 patients with clonal group 258 CRKp were enrolled in the Consortium on Resistance Against Carbapenems in Klebsiella and other Enterobacteriaceae, a prospective, multicenter, cohort study. A maximum likelihood tree was constructed using RAxML. Static clusters shared ≤21 single-nucleotide polymorphisms (SNP) and a most recent common ancestor. Dynamic clusters incorporated SNP distance, culture timing, and rates of SNP accumulation and transmission using the R program TransCluster.

RESULTS

Most patients were admitted from home (n = 150, 43%) or long-term care facilities (n = 115, 33%). Urine (n = 149, 43%) was the most common isolation site. Overall, 55 static and 47 dynamics clusters were identified involving 210 of 350 (60%) and 194 of 350 (55%) patients, respectively. Approximately half of static clusters were identical to dynamic clusters. Static clusters consisted of 33 (60%) intrasystem and 22 (40%) intersystem clusters. Dynamic clusters consisted of 32 (68%) intrasystem and 15 (32%) intersystem clusters and had fewer SNP differences than static clusters (8 vs 9; P = .045; 95% confidence interval [CI]: -4 to 0). Dynamic intersystem clusters contained more patients than dynamic intrasystem clusters (median [interquartile range], 4 [2, 7] vs 2 [2, 2]; P = .007; 95% CI: -3 to 0).

CONCLUSIONS

Widespread intrasystem and intersystem transmission of CRKp was identified in hospitalized US patients. Use of different methods for assessing genetic similarity resulted in only minor differences in interpretation.