Whole-Genome Sequencing To Identify Drivers of Carbapenem-Resistant Klebsiella pneumoniae Transmission within and between Regional Long-Term Acute-Care Hospitals

Emergence and persistent spread of carbapenemase-producing Klebsiella pneumoniae high-risk clones in Greek hospitals, 2013 to 2022

BackgroundPreliminary unpublished results of the survey of carbapenem- and/or colistin-resistant Enterobacterales (CCRE survey) showed the expansion of carbapenemase-producing Klebsiella pneumoniae (CPKP) sequence type (ST) 39 in 12 of 15 participating Greek hospitals in 2019.AimWe conducted a rapid survey to determine the extent of spread of CPKP high-risk clones in Greek hospitals in 2022 and compare the distribution of circulating CPKP clones in these hospitals since 2013.MethodsWe analysed whole genome sequences and epidemiological data of 310 K. pneumoniae isolates that were carbapenem-resistant or 'susceptible, increased exposure' from Greek hospitals that participated in the European survey of carbapenemase-producing Enterobacteriaceae (EuSCAPE, 2013-2014), in the CCRE survey (2019) and in a national follow-up survey (2022) including, for the latter, an estimation of transmission events.ResultsFive K. pneumoniae STs including ST258/512 (n = 101 isolates), ST11 (n = 93), ST39 (n = 56), ST147 (n = 21) and ST323 (n = 13) accounted for more than 90% of CPKP isolates in the dataset. While ST11, ST147 and ST258/512 have been detected in participating hospitals since 2013 and 2014, KPC-2-producing ST39 and ST323 emerged in 2019 and 2022, respectively. Based on the defined genetic relatedness cut-off, 44 within-hospital transmission events were identified in the 2022 survey dataset, with 12 of 15 participating hospitals having at least one within-hospital transmission event.ConclusionThe recent emergence and rapid spread of new high-risk K. pneumoniae clones in the Greek healthcare system related to within-hospital transmission is of concern and highlights the need for molecular surveillance and enhanced infection prevention and control measures.

The ICU-CARB score: a novel clinical scoring system to predict carbapenem-resistant gram-negative bacteria carriage in critically ill patients upon ICU admission

BACKGROUND

With the widespread spread of carbapenem-resistant gram-negative bacteria (CR-GNB) in medical facilities, the carriage of CR-GNB among critically ill patients has become a significant concern in intensive care units (ICU). This study aimed to develop a scoring system to identify CR-GNB carriers upon ICU admission.

METHODS

Consecutive critically ill patients admitted to the ICU of Shanghai Ruijin Hospital between January 2017 and December 2020 were included. The patients were then divided into training and testing datasets at a 7:3 ratio. Parameters associated with CR-GNB carriage were identified using least absolute shrinkage and selection operator regression analysis. Each parameter was assigned a numerical score ranging from 0 to 100 using logistic regression analysis. Subsequently, a four-tier risk-level system was developed based on the cumulative scores, and assessed using the area under the receiver operating characteristic curve (AUC).

RESULTS

Of the 1736 patients included in this study, the prevalence of CR-GNB carriage was 10.60%. The clinical scoring system including seven variables (neurological disease, high-risk department history, length of stay ≥ 14 days, ICU history, invasive mechanical ventilation, gastrointestinal tube placement, and carbapenem usage) exhibited promising predictive capabilities. Patients were then stratified using the scoring system, resulting in CR-GNB carriage rates of 2.4%, 12.0%, 36.1%, and 57.9% at the respective risk levels (P < 0.001). Furthermore, the AUC of the developed model in the training set was calculated to be 0.82 (95% CI, 0.78-0.86), while internal validation yielded an AUC of 0.83 (95% CI, 0.77-0.89).

CONCLUSIONS

The ICU-CARB Score serves as a straightforward and precise tool that enables prompt evaluation of the risk of CR-GNB carriage at the time of ICU admission, thereby facilitating the timely implementation of targeted pre-emptive isolation.

Distinct Origins and Transmission Pathways of blaKPC Enterobacterales across Three U.S. States

Carbapenem-resistant Enterobacterales (CRE) are among the most concerning antibiotic resistance threats due to high rates of multidrug resistance, transmissibility in health care settings, and high mortality rates. We evaluated the potential for regional genomic surveillance to track the spread of blaKPC-carrying CRE (KPC-CRE) by using isolate collections from health care facilities in three U.S. states. Clinical isolates were collected from Connecticut (2017 to 2018), Minnesota (2012 to 2018), and Tennessee (2016 to 2017) through the U.S. Centers for Disease Control and Prevention's Multi-site Gram-negative Surveillance Initiative (MuGSI) and additional surveillance. KPC-CRE isolates were whole-genome sequenced, yielding 255 isolates from 214 patients across 96 facilities. Case report data on patient comorbidities, facility exposures, and interfacility patient transfer were extracted. We observed that in Connecticut, most KPC-CRE isolates showed evidence of importation from outside the state, with limited local transmission. In Minnesota, cases were mainly from sporadic importation and transmission of blaKPC-carrying Klebsiella pneumoniae ST258, and clonal expansion of blaKPC-carrying Enterobacter hormaechei ST171, primarily at a single focal facility and its satellite facilities. In Tennessee, we observed transmission of diverse strains of blaKPC-carrying Enterobacter and Klesbiella, with evidence that most derived from the local acquisition of blaKPC plasmids circulating in an interconnected regional health care network. Thus, the underlying processes driving KPC-CRE burden can differ substantially across regions and can be discerned through regional genomic surveillance. This study provides proof of concept that integrating genomic data with information on interfacility patient transfers can provide insights into locations and drivers of regional KPC-CRE burden that can enable targeted interventions.

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.

Coexistence of tet(A) and blaKPC-2 in the ST11 hypervirulent tigecycline- and carbapenem-resistant Klebsiella pneumoniae isolated from a blood sample

Carbapenem-resistant Klebsiella pneumoniae are distributed worldwide. This study aimed to characterize a hypervirulent tigecycline-resistant and carbapenem-resistant Klebsiella pneumoniae strain, XJ-K2, collected from a patient's blood. We tested antimicrobial susceptibility, virulence, and whole-genome sequencing (WGS) on strain XJ-K2. WGS data were used to identify virulence and resistance genes and to perform multilocus sequence typing (MLST) and phylogenetic analysis. Three novel plasmids, including a pLVPK-like virulence plasmid (pXJ-K2-p1) and two multiple resistance plasmids (pXJ-K2-KPC-2 and pXJ-K2-p3), were discovered in strain XJ-K2. The IncFII(pCRY) plasmid pXJ-K2-p3 carried the dfrA14, sul2, qnrS1, bla_LAP-2, and tet(A) resistance genes. The IncFII(pHN7A8)/IncR plasmid pXJ-K2-KPC-2 also carried a range of resistance elements, containing rmtB, bla_KPC-2, bla_TEM-1, bla_CTX-M-65, and fosA3. MLST analysis revealed that strain XJ-K2 belonged to sequence type 11 (ST11). Seven complete phage sequences and many virulence genes were found in strain XJ-K2. Meanwhile, antimicrobial susceptibility tests and G. mellonella larval infection models confirmed the extensively drug resistance (XDR) and hypervirulence of KJ-K2. To our knowledge, this is the first observation and description of the ST11 hypervirulent tigecycline- and carbapenem-resistant K. pneumoniae strain co-carrying bla_KPC-2 and the tet(A) in a patient's blood in China. Further investigation is needed to understand the resistance and virulence mechanisms of this significant hypervirulent tigecycline- and carbapenem-resistant strain.

Risk Factors for Colistin-Resistant Carbapenem-Resistant Klebsiella pneumoniae in the Postacute Care Setting

We assessed risk factors for colistin resistance among carbapenem-resistant Klebsiella pneumoniae (CRKP) from 375 patients in long-term acute care hospitals. Recent colistin or polymyxin B exposure was associated with increased odds of colistin resistance (adjusted odds ratio = 1.11 per day of exposure, 95% confidence interval = 1.03-1.19, P = .007).

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

BACKGROUND

A crucial barrier to the routine application of whole-genome sequencing (WGS) for infection prevention is the insufficient criteria for determining whether a genomic linkage is consistent with transmission within the facility. We evaluated the use of single-nucleotide variant (SNV) thresholds, as well as a novel threshold-free approach, for inferring transmission linkages in a high-transmission setting.

METHODS

We did a retrospective genomic epidemiology analysis of samples previously collected in the context of an intervention study at a long-term acute care hospital in the USA. We performed WGS on 435 isolates of Klebsiella pneumoniae harbouring the bla_KPC carbapenemase (KPC-Kp) collected from 256 patients through admission and surveillance culturing (once every 2 weeks) of almost every patient who was admitted to hospital over a 1-year period.

FINDINGS

Our analysis showed that the standard approach of using an SNV threshold to define transmission would lead to false-positive and false-negative inferences. False-positive inferences were driven by the frequent importation of closely related strains, which were presumably linked via transmission at connected health-care facilities. False-negative inferences stemmed from the diversity of colonising populations that were spread among patients, with multiple examples of hypermutator strain emergence within patients and, as a result, putative transmission links separated by large genetic distances. Motivated by limitations of an SNV threshold, we implemented a novel threshold-free transmission cluster inference approach, in which each of the acquired KPC-Kp isolates were linked back to the imported KPC-Kp isolate with which it shared the most variants. This approach yielded clusters that varied in levels of genetic diversity but where 105 (81%) of 129 unique strain acquisition events were associated with epidemiological links in the hospital. Of 100 patients who acquired KPC-Kp isolates that were included in a cluster, 47 could be linked to a single patient who was positive for KPC-Kp at admission, compared with 31 and 25 using 10 SNV and 20 SNV thresholds, respectively. Holistic examination of clusters highlighted extensive variation in the magnitude of onward transmission stemming from more than 100 importation events and revealed patterns in cluster propagation that could inform improvements to infection prevention strategies.

INTERPRETATION

Our results show how the integration of culture surveillance data into genomic analyses can overcome limitations of cluster detection based on SNV-thresholds and improve the ability to track pathways of pathogen transmission in health-care settings.

FUNDING

US Center for Disease Control and Prevention and University of Michigan.

Inter-species geographic signatures for tracing horizontal gene transfer and long-term persistence of carbapenem resistance

BACKGROUND

Carbapenem-resistant Enterobacterales (CRE) are an urgent global health threat. Inferring the dynamics of local CRE dissemination is currently limited by our inability to confidently trace the spread of resistance determinants to unrelated bacterial hosts. Whole-genome sequence comparison is useful for identifying CRE clonal transmission and outbreaks, but high-frequency horizontal gene transfer (HGT) of carbapenem resistance genes and subsequent genome rearrangement complicate tracing the local persistence and mobilization of these genes across organisms.

METHODS

To overcome this limitation, we developed a new approach to identify recent HGT of large, near-identical plasmid segments across species boundaries, which also allowed us to overcome technical challenges with genome assembly. We applied this to complete and near-complete genome assemblies to examine the local spread of CRE in a systematic, prospective collection of all CRE, as well as time- and species-matched carbapenem-susceptible Enterobacterales, isolated from patients from four US hospitals over nearly 5 years.

RESULTS

Our CRE collection comprised a diverse range of species, lineages, and carbapenem resistance mechanisms, many of which were encoded on a variety of promiscuous plasmid types. We found and quantified rearrangement, persistence, and repeated transfer of plasmid segments, including those harboring carbapenemases, between organisms over multiple years. Some plasmid segments were found to be strongly associated with specific locales, thus representing geographic signatures that make it possible to trace recent and localized HGT events. Functional analysis of these signatures revealed genes commonly found in plasmids of nosocomial pathogens, such as functions required for plasmid retention and spread, as well survival against a variety of antibiotic and antiseptics common to the hospital environment.

CONCLUSIONS

Collectively, the framework we developed provides a clearer, high-resolution picture of the epidemiology of antibiotic resistance importation, spread, and persistence in patients and healthcare networks.

A review of epidemiological models of Clostridioides difficile transmission and control (2009-2021)

Clostridioides difficile is the leading cause of infectious diarrhea and one of the most common healthcare-acquired infections worldwide. We performed a systematic search and a bibliometric analysis of mathematical and computational models for Clostridioides difficile transmission. We identified 33 publications from 2009 to 2021. Models have underscored the importance of asymptomatic colonized patients in maintaining transmission in health-care settings. Infection control, antimicrobial stewardship, active testing, and vaccination have often been evaluated in models. Despite active testing and vaccination being not currently implemented, they are the most commonly evaluated interventions. Some aspects of C. difficile transmission, such community transmission and interventions in health-care settings other than in acute-care hospitals, remained less evaluated through modeling.

regentrans: a framework and R package for using genomics to study regional pathogen transmission

Increasing evidence of regional pathogen transmission networks highlights the importance of investigating the dissemination of multidrug-resistant organisms (MDROs) across a region to identify where transmission is occurring and how pathogens move across regions. We developed a framework for investigating MDRO regional transmission dynamics using whole-genome sequencing data and created regentrans, an easy-to-use, open source R package that implements these methods (https://github.com/Snitkin-Lab-Umich/regentrans). Using a dataset of over 400 carbapenem-resistant isolates of Klebsiella pneumoniae collected from patients in 21 long-term acute care hospitals over a one-year period, we demonstrate how to use our framework to gain insights into differences in inter- and intra-facility transmission across different facilities and over time. This framework and corresponding R package will allow investigators to better understand the origins and transmission patterns of MDROs, which is the first step in understanding how to stop transmission at the regional level.

High-Resolution Genomic Profiling of Carbapenem-Resistant Klebsiella pneumoniae Isolates: A Multicentric Retrospective Indian Study

BACKGROUND

Carbapenem-resistant Klebsiella pneumoniae (CRKP) is a threat to public health in India because of its high dissemination, mortality, and limited treatment options. Its genomic variability is reflected in the diversity of sequence types, virulence factors, and antimicrobial resistance (AMR) mechanisms. This study aims to characterize the clonal relationships and genetic mechanisms of resistance and virulence in CRKP isolates in India.

MATERIALS AND METHODS

We characterized 344 retrospective K. pneumoniae clinical isolates collected from 8 centers across India collected in 2013-2019. Susceptibility to antibiotics was tested with VITEK 2. Capsular types, multilocus sequence type, virulence genes, AMR determinants, plasmid replicon types, and a single-nucleotide polymorphism phylogeny were inferred from their whole genome sequences.

RESULTS

Phylogenetic analysis of the 325 Klebsiella isolates that passed quality control revealed 3 groups: K. pneumoniae sensu stricto (n = 307), K. quasipneumoniae (n = 17), and K. variicola (n = 1). Sequencing and capsular diversity analysis of the 307 K. pneumoniae sensu stricto isolates revealed 28 sequence types, 26 K-locus types, and 11 O-locus types, with ST231, KL51, and O1V2 being predominant. blaOXA-48-like and blaNDM-1/5 were present in 73.2% and 24.4% of isolates, respectively. The major plasmid replicon types associated with carbapenase genes were IncF (51.0%) and Col group (35.0%).

CONCLUSION

Our study documents for the first time the genetic diversity of K and O antigens circulating in India. The results demonstrate the practical applicability of genomic surveillance and its utility in tracking the population dynamics of CRKP. It alerts us to the urgency for longitudinal surveillance of these transmissible lineages.

Epidemiology and Transmission of Carbapenemase-Producing Enterobacteriaceae in a Health Care Network of an Acute-Care Hospital and Its Affiliated Intermediate- and Long-Term-Care Facilities in Singapore

Movement of patients in a health care network poses challenges for the control of carbapenemase-producing Enterobacteriaceae (CPE). We aimed to identify intra- and interfacility transmission events and facility type-specific risk factors of CPE in an acute-care hospital (ACH) and its intermediate-term and long-term-care facilities (ILTCFs). Serial cross-sectional studies were conducted in June and July of 2014 to 2016 to screen for CPE. Whole-genome sequencing was done to identify strain relatedness and CPE genes (bla_IMI, bla_IMP-1, bla_KPC-2, bla_NDM-1, and bla_OXA-48). Multivariable logistic regression models, stratified by facility type, were used to determine independent risk factors. Of 5,357 patients, half (55%) were from the ACH. CPE prevalence was 1.3% in the ACH and 0.7% in ILTCFs (P = 0.029). After adjusting for sociodemographics, screening year, and facility type, the odds of CPE colonization increased significantly with a hospital stay of ≥3 weeks (adjusted odds ratio [aOR], 2.67; 95% confidence interval [CI], 1.17 to 6.05), penicillin use (aOR, 3.00; 95% CI, 1.05 to 8.56), proton pump inhibitor use (aOR, 3.20; 95% CI, 1.05 to 9.80), dementia (aOR, 3.42; 95% CI, 1.38 to 8.49), connective tissue disease (aOR, 5.10; 95% CI, 1.19 to 21.81), and prior carbapenem-resistant Enterobacteriaceae (CRE) carriage (aOR, 109.02; 95% CI, 28.47 to 417.44) in the ACH. For ILTCFs, presence of wounds (aOR, 5.30; 95% CI, 1.01 to 27.72), respiratory procedures (aOR, 4.97; 95% CI, 1.09 to 22.71), vancomycin-resistant enterococcus carriage (aOR, 16.42; 95% CI, 1.52 to 177.48), and CRE carriage (aOR, 758.30; 95% CI, 33.86 to 16,982.52) showed significant association. Genomic analysis revealed only possible intra-ACH transmission and no evidence for ACH-to-ILTCF transmission. Although CPE colonization was predominantly in the ACH, risk factors varied between facilities. Targeted screening and precautionary measures are warranted.

Validation and Extrapolation of a Multimodal Infection Prevention and Control Intervention on Carbapenem-Resistant Klebsiella pneumoniae in an Epidemic Region: A Historical Control Quasi-Experimental Study

Objective: To verify the effects of comprehensive infection prevention and control (IPC) interventions for the prevention of the cross-transmission of carbapenem-resistant Klebsiella pneumoniae (CRKP) within intensive care units (ICUs) in an epidemic region.

Methods: A historical control, quasi-experimental design was performed. The study was conducted between January 2017 and December 2019, following the implementation of a multimodal IPC bundle. The baseline period was established from January 2013 to June 2013, when only basic IPC measures were applied.

Results: A total of 748 patients were enrolled during the entire study. The incidence of ICU-acquired CRKP colonization/infection was 1.16 per 1,000 patient-days during the intervention period, compared with 10.19 per 1,000 patient-days during the baseline period (p = 0.002). The slope of the monthly incidence of CRKP at admission showed an increasing trend (p = 0.03). The incidence of ICU-acquired catheter-related bloodstream infections caused by CRKP decreased from 2.54 to 0.96 per 1,000 central-line-days (p = 0.08). Compliance with contact precautions and terminal room disinfection improved during the intervention period. All environmental surface culture samples acquired after terminal room disinfection were negative for CRKP.

Conclusion: Our findings suggest that in epidemic settings, multimodal IPC intervention strategies and consistent monitoring of compliance, may limit the spread of CRKP in ICUs.

A multisite genomic epidemiology study of Clostridioides difficile infections in the USA supports differential roles of healthcare versus community spread for two common strains

Clostridioides difficile is the leading cause of healthcare-associated infectious diarrhoea. However, it is increasingly appreciated that healthcare-associated infections derive from both community and healthcare environments, and that the primary sites of C. difficile transmission may be strain-dependent. We conducted a multisite genomic epidemiology study to assess differential genomic evidence of healthcare vs community spread for two of the most common C. difficile strains in the USA: sequence type (ST) 1 (associated with ribotype 027) and ST2 (associated with ribotype 014/020). We performed whole-genome sequencing and phylogenetic analyses on 382 ST1 and ST2 C. difficile isolates recovered from stool specimens collected during standard clinical care at 3 geographically distinct US medical centres between 2010 and 2017. ST1 and ST2 isolates both displayed some evidence of phylogenetic clustering by study site, but clustering was stronger and more apparent in ST1, consistent with our healthcare-based study more comprehensively sampling local transmission of ST1 compared to ST2 strains. Analyses of pairwise single-nucleotide variant (SNV) distance distributions were also consistent with more evidence of healthcare transmission of ST1 compared to ST2, with 44 % of ST1 isolates being within two SNVs of another isolate from the same geographical collection site compared to 5.5 % of ST2 isolates (P-value=<0.001). Conversely, ST2 isolates were more likely to have close genetic neighbours across disparate geographical sites compared to ST1 isolates, further supporting non-healthcare routes of spread for ST2 and highlighting the potential for misattributing genomic similarity among ST2 isolates to recent healthcare transmission. Finally, we estimated a lower evolutionary rate for the ST2 lineage compared to the ST1 lineage using Bayesian timed phylogenomic analyses, and hypothesize that this may contribute to observed differences in geographical concordance among closely related isolates. Together, these findings suggest that ST1 and ST2, while both common causes of C. difficile infection in hospitals, show differential reliance on community and hospital spread. This conclusion supports the need for strain-specific criteria for interpreting genomic linkages and emphasizes the importance of considering differences in the epidemiology of circulating strains when devising interventions to reduce the burden of C. difficile infections.

Regional Spread of blaNDM-1-Containing Klebsiella pneumoniae ST147 in Post-Acute Care Facilities

BACKGROUND

Carbapenem-resistant Enterobacterales (CRE) harboring blaKPC have been endemic in Chicago-area healthcare networks for more than a decade. During 2016-2019, a series of regional point-prevalence surveys identified increasing prevalence of blaNDM-containing CRE in multiple long-term acute care hospitals (LTACHs) and ventilator-capable skilled nursing facilities (vSNFs). We performed a genomic epidemiology investigation of blaNDM-producing CRE to understand their regional emergence and spread.

METHODS

We performed whole-genome sequencing on New Delhi metallo-beta-lactamase (NDM)+ CRE isolates from 4 point-prevalence surveys across 35 facilities (LTACHs, vSNFs, and acute care hospital medical intensive care units) in the Chicago area and investigated the genomic relatedness and transmission dynamics of these isolates over time.

RESULTS

Genomic analyses revealed that the rise of NDM+ CRE was due to the clonal dissemination of an sequence type (ST) 147 Klebsiella pneumoniae strain harboring blaNDM-1 on an IncF plasmid. Dated phylogenetic reconstructions indicated that ST147 was introduced into the region around 2013 and likely acquired NDM around 2015. Analyzing the relatedness of strains within and between facilities supported initial increases in prevalence due to intrafacility transmission in certain vSNFs, with evidence of subsequent interfacility spread among LTACHs and vSNFs connected by patient transfer.

CONCLUSIONS

We identified a regional outbreak of blaNDM-1 ST147 that began in and disseminated across Chicago area post-acute care facilities. Our findings highlight the importance of performing genomic surveillance at post-acute care facilities to identify emerging threats.

Patient and Microbial Genomic Factors Associated with Carbapenem-Resistant Klebsiella pneumoniae Extraintestinal Colonization and Infection

Carbapenem-resistant Klebsiella pneumoniae (CRKP) is a critical-priority antibiotic resistance threat that has emerged over the past several decades, spread across the globe, and accumulated resistance to last-line antibiotic agents. While CRKP infections are associated with high mortality, only a subset of patients acquiring CRKP extraintestinal colonization will develop clinical infection. Here, we sought to ascertain the relative importance of patient characteristics and CRKP genetic background in determining patient risk of infection. Machine learning models classifying colonization versus infection were built using whole-genome sequences and clinical metadata from a comprehensive set of 331 CRKP extraintestinal isolates collected across 21 long-term acute-care hospitals over the course of a year. Model performance was evaluated based on area under the receiver operating characteristic curve (AUROC) on held-out test data. We found that patient and genomic features were predictive of clinical CRKP infection to similar extents (AUROC interquartile ranges [IQRs]: patient = 0.59 to 0.68, genomic = 0.55 to 0.61, combined = 0.62 to 0.68). Patient predictors of infection included the presence of indwelling devices, kidney disease, and length of stay. Genomic predictors of infection included presence of the ICEKp10 mobile genetic element carrying the yersiniabactin iron acquisition system and disruption of an O-antigen biosynthetic gene in a sublineage of the epidemic ST258 clone. Altered O-antigen biosynthesis increased association with the respiratory tract, and subsequent ICEKp10 acquisition was associated with increased virulence. These results highlight the potential of integrated models including both patient and microbial features to provide a more holistic understanding of patient clinical trajectories and ongoing within-lineage pathogen adaptation.IMPORTANCE Multidrug-resistant organisms, such as carbapenem-resistant Klebsiella pneumoniae (CRKP), colonize alarmingly large fractions of patients in regions of endemicity, but only a subset of patients develop life-threatening infections. While patient characteristics influence risk for infection, the relative contribution of microbial genetic background to patient risk remains unclear. We used machine learning to determine whether patient and/or microbial characteristics can discriminate between CRKP extraintestinal colonization and infection across multiple health care facilities and found that both patient and microbial factors were predictive. Examination of informative microbial genetic features revealed variation within the ST258 epidemic lineage that was associated with respiratory tract colonization and increased rates of infection. These findings indicate that circulating genetic variation within a highly prevalent epidemic lineage of CRKP influences patient clinical trajectories. In addition, this work supports the need for future studies examining the microbial genetic determinants of clinical outcomes in human populations, as well as epidemiologic and experimental follow-ups of identified features to discern generalizability and biological mechanisms.

Cohorting KPC+ Klebsiella pneumoniae (KPC-Kp)-positive patients: A genomic exposé of cross-colonization hazards in a long-term acute-care hospital (LTACH)

OBJECTIVE

Cohorting patients who are colonized or infected with multidrug-resistant organisms (MDROs) protects uncolonized patients from acquiring MDROs in healthcare settings. The potential for cross transmission within the cohort and the possibility of colonized patients acquiring secondary isolates with additional antibiotic resistance traits is often neglected. We searched for evidence of cross transmission of KPC+ Klebsiella pneumoniae (KPC-Kp) colonization among cohorted patients in a long-term acute-care hospital (LTACH), and we evaluated the impact of secondary acquisitions on resistance potential.

DESIGN

Genomic epidemiological investigation.

SETTING

A high-prevalence LTACH during a bundled intervention that included cohorting KPC-Kp-positive patients.

METHODS

Whole-genome sequencing (WGS) and location data were analyzed to identify potential cases of cross transmission between cohorted patients.

RESULTS

Secondary KPC-Kp isolates from 19 of 28 admission-positive patients were more closely related to another patient's isolate than to their own admission isolate. Of these 19 cases, 14 showed strong genomic evidence for cross transmission (<10 single nucleotide variants or SNVs), and most of these patients occupied shared cohort floors (12 patients) or rooms (4 patients) at the same time. Of the 14 patients with strong genomic evidence of acquisition, 12 acquired antibiotic resistance genes not found in their primary isolates.

CONCLUSIONS

Acquisition of secondary KPC-Kp isolates carrying distinct antibiotic resistance genes was detected in nearly half of cohorted patients. These results highlight the importance of healthcare provider adherence to infection prevention protocols within cohort locations, and they indicate the need for future studies to assess whether multiple-strain acquisition increases risk of adverse patient outcomes.

prewas: data pre-processing for more informative bacterial GWAS

While variant identification pipelines are becoming increasingly standardized, less attention has been paid to the pre-processing of variants prior to their use in bacterial genome-wide association studies (bGWAS). Three nuances of variant pre-processing that impact downstream identification of genetic associations include the separation of variants at multiallelic sites, separation of variants in overlapping genes, and referencing of variants relative to ancestral alleles. Here we demonstrate the importance of these variant pre-processing steps on diverse bacterial genomic datasets and present prewas, an R package, that standardizes the pre-processing of multiallelic sites, overlapping genes, and reference alleles before bGWAS. This package facilitates improved reproducibility and interpretability of bGWAS results. prewas enables users to extract maximal information from bGWAS by implementing multi-line representation for multiallelic sites and variants in overlapping genes. prewas outputs a binary SNP matrix that can be used for SNP-based bGWAS and will prevent the masking of minor alleles during bGWAS analysis. The optional binary gene matrix output can be used for gene-based bGWAS, which will enable users to maximize the power and evolutionary interpretability of their bGWAS studies. prewas is available for download from GitHub.