Development of a novel streamlined workflow (AACRE) and database (inCREDBle) for genomic analysis of carbapenem-resistant Enterobacterales

In response to the threat of increasing antimicrobial resistance, we must increase the amount of available high-quality genomic data gathered on antibiotic-resistant bacteria. To this end, we developed an integrated pipeline for high-throughput long-read sequencing, assembly, annotation and analysis of bacterial isolates and used it to generate a large genomic data set of carbapenemase-producing Enterobacterales (CPE) isolates collected in Spain. The set of 461 isolates were sequenced with a combination of both Illumina and Oxford Nanopore Technologies (ONT) DNA sequencing technologies in order to provide genomic context for chromosomal loci and, most importantly, structural resolution of plasmids, important determinants for transmission of antimicrobial resistance. We developed an informatics pipeline called Assembly and Annotation of Carbapenem-Resistant Enterobacteriaceae (AACRE) for the full assembly and annotation of the bacterial genomes and their complement of plasmids. To explore the resulting genomic data set, we developed a new database called inCREDBle that not only stores the genomic data, but provides unique ways to filter and compare data, enabling comparative genomic analyses at the level of chromosomes, plasmids and individual genes. We identified a new sequence type, ST5000, and discovered a genomic locus unique to ST15 that may be linked to its increased spread in the population. In addition to our major objective of generating a large regional data set, we took the opportunity to compare the effects of sample quality and sequencing methods, including R9 versus R10 nanopore chemistry, on genome assembly and annotation quality. We conclude that converting short-read and hybrid microbial sequencing and assembly workflows to the latest nanopore chemistry will further reduce processing time and cost, truly enabling the routine monitoring of resistance transmission patterns at the resolution of complete chromosomes and plasmids.

Emergence of 16S rRNA methyltransferases among carbapenemase-producing Enterobacterales in Spain studied by whole-genome sequencing

The emergence of 16S rRNA methyltransferases (RMTs) in Gram-negative pathogens bearing other clinically relevant resistance mechanisms, such as carbapenemase-producing Enterobacterales (CPE), is becoming an alarming concern. We investigated the prevalence, antimicrobial susceptibility, resistance mechanisms, molecular epidemiology and genetic support of RMTs in CPE isolates from Spain. This study included a collection of 468 CPE isolates recovered during 2018 from 32 participating Spanish hospitals. MICs were determined using the broth microdilution method, the agar dilution method (fosfomycin) or MIC gradient strips (plazomicin). All isolates were subjected to hybrid whole-genome sequencing (WGS). Sequence types (STs), core genome phylogenetic relatedness, horizontally acquired resistance mechanisms, plasmid analysis and the genetic environment of RMTs were determined in silico from WGS data in all RMT-positive isolates. Among the 468 CPE isolates evaluated, 24 isolates (5.1%) recovered from nine different hospitals spanning five Spanish regions showed resistance to all aminoglycosides and were positive for an RMT (21 RmtF, 2 ArmA and 1 RmtC). All RMT-producers showed high-level resistance to all aminoglycosides, including plazomicin, and in most cases exhibited an extensively drug-resistant susceptibility profile. The RMT-positive isolates showed low genetic diversity and were global clones of Klebsiella pneumoniae (ST147, ST101, ST395) and Enterobacter cloacae (ST93) bearing bla_OXA-48, bla_NDM-1 or bla_VIM-1 carbapenemase genes. RMTs were harboured in five different multidrug resistance plasmids and linked to efficient mobile genetic elements. Our findings highlight that RMTs are emerging among clinical CPE isolates from Spain and their spread should be monitored to preserve the future clinical utility of aminoglycosides and plazomicin.

Activity of imipenem/relebactam against a Spanish nationwide collection of carbapenemase-producing Enterobacterales

BACKGROUND

Imipenem/relebactam is a novel carbapenem/β-lactamase inhibitor combination, developed to act against carbapenemase-producing Enterobacterales (CPE).

OBJECTIVES

To assess the in vitro activity of imipenem/relebactam against a Spanish nationwide collection of CPE by testing the susceptibility of these isolates to 16 widely used antimicrobials and to determine the underlying β-lactam resistance mechanisms involved and the molecular epidemiology of carbapenemases in Spain.

MATERIALS AND METHODS

Clinical CPE isolates (n = 401) collected for 2 months from 24 hospitals in Spain were tested. MIC50, MIC90 and susceptibility/resistance rates were interpreted in accordance with the EUCAST guidelines. β-Lactam resistance mechanisms and molecular epidemiology were characterized by WGS.

RESULTS

For all isolates, high rates of susceptibility to colistin (86.5%; MIC50/90 = 0.12/8 mg/L), imipenem/relebactam (85.8%; MIC50/90 = 0.5/4 mg/L) and ceftazidime/avibactam (83.8%, MIC50/90 = 1/≥256 mg/L) were observed. The subgroups of isolates producing OXA-48-like (n = 305, 75.1%) and KPC-like enzymes (n = 44, 10.8%) were highly susceptible to ceftazidime/avibactam (97.7%, MIC50/90 = 1/2 mg/L) and imipenem/relebactam (100.0%, MIC50/90 = ≤0.25/1 mg/L), respectively.The most widely disseminated high-risk clones of carbapenemase-producing Klebsiella pneumoniae across Spain were found to be ST11, ST147, ST392 and ST15 (mostly associated with OXA-48) and ST258/512 (in all cases producing KPC).

CONCLUSIONS

Imipenem/relebactam, colistin and ceftazidime/avibactam were the most active antimicrobials against all CPEs. Imipenem/relebactam is a valuable addition to the antimicrobial arsenal used in the fight against CPE, particularly against KPC-producing isolates, which in all cases were susceptible to this combination.

Improved prediction of salvage antiretroviral therapy outcomes using ultrasensitive HIV-1 drug resistance testing

BACKGROUND

The clinical relevance of ultrasensitive human immunodeficiency virus type 1 (HIV-1) genotypic resistance testing in antiretroviral treatment (ART)-experienced individuals remains unknown.

METHODS

This was a retrospective, multicentre, cohort study in ART-experienced, HIV-1-infected adults who initiated salvage ART including, at least 1 ritonavir-boosted protease inhibitor, raltegravir or etravirine. Presalvage ART Sanger and 454 sequencing of plasma HIV-1 were used to generate separate genotypic sensitivity scores (GSS) using the HIVdb, ANRS, and REGA algorithms. Virological failure (VF) was defined as 2 consecutive HIV-1 RNA levels ≥200 copies/mL at least 12 weeks after salvage ART initiation, whereas subjects remained on the same ART. The ability of Sanger and 454-GSS to predict VF was assessed by receiver operating characteristic (ROC) curves and survival analyses.

RESULTS

The study included 132 evaluable subjects; 28 (21%) developed VF. Using HIVdb, 454 predicted VF better than Sanger sequencing in the ROC curve analysis (area under the curve: 0.69 vs 0.60, Delong test P = .029). Time to VF was shorter for subjects with 454-GSS < 3 vs 454-GSS ≥ 3 (Log-rank P = .003) but not significantly different between Sanger-GSS < 3 and ≥3. Factors independently associated with increased risk of VF in multivariate Cox regression were a 454-GSS < 3 (HR = 4.6, 95 CI, [1.5, 14.0], P = .007), and the number of previous antiretrovirals received (HR = 1.2 per additional drug, 95 CI, [1.1, 1.3], P = .001). Equivalent findings were obtained with the ANRS and REGA algorithms.

CONCLUSIONS

Ultrasensitive HIV-1 genotyping improves GSS-based predictions of virological outcomes of salvage ART relative to Sanger sequencing. This may improve the clinical management of ART-experienced subjects living with HIV-1.

CLINICAL TRIALS REGISTRATION

NCT01346878.

Assessing differential expression measurements by highly parallel pyrosequencing and DNA microarrays: a comparative study

To explore the feasibility of pyrosequencing for quantitative differential gene expression analysis we have performed a comparative study of the results of the sequencing experiments to those obtained by a conventional DNA microarray platform. A conclusion from our analysis is that, over a threshold of 35 normalized reads per gene, the measurements of gene expression display a good correlation with the references. The observed concordance between pyrosequencing and DNA microarray platforms beyond the threshold was of 0.8, measured as a Pearson's correlation coefficient. In differential gene expression the initial aim is the quantification the differences among transcripts when comparing experimental conditions. Thus, even in a scenario of low coverage the concordance in the measurements is quite acceptable. On the other hand, the comparatively longer read size obtained by pyrosequencing allows detecting unconventional splicing forms.

Added value of deep sequencing relative to population sequencing in heavily pre-treated HIV-1-infected subjects

OBJECTIVE

To explore the potential of deep HIV-1 sequencing for adding clinically relevant information relative to viral population sequencing in heavily pre-treated HIV-1-infected subjects.

METHODS

In a proof-of-concept study, deep sequencing was compared to population sequencing in HIV-1-infected individuals with previous triple-class virological failure who also developed virologic failure to deep salvage therapy including, at least, darunavir, tipranavir, etravirine or raltegravir. Viral susceptibility was inferred before salvage therapy initiation and at virological failure using deep and population sequencing genotypes interpreted with the HIVdb, Rega and ANRS algorithms. The threshold level for mutant detection with deep sequencing was 1%.

RESULTS

7 subjects with previous exposure to a median of 15 antiretrovirals during a median of 13 years were included. Deep salvage therapy included darunavir, tipranavir, etravirine or raltegravir in 4, 2, 2 and 5 subjects, respectively. Self-reported treatment adherence was adequate in 4 and partial in 2; one individual underwent treatment interruption during follow-up. Deep sequencing detected all mutations found by population sequencing and identified additional resistance mutations in all but one individual, predominantly after virological failure to deep salvage therapy. Additional genotypic information led to consistent decreases in predicted susceptibility to etravirine, efavirenz, nucleoside reverse transcriptase inhibitors and indinavir in 2, 1, 2 and 1 subject, respectively. Deep sequencing data did not consistently modify the susceptibility predictions achieved with population sequencing for darunavir, tipranavir or raltegravir.

CONCLUSIONS

In this subset of heavily pre-treated individuals, deep sequencing improved the assessment of genotypic resistance to etravirine, but did not consistently provide additional information on darunavir, tipranavir or raltegravir susceptibility. These data may inform the design of future studies addressing the clinical value of minority drug-resistant variants in treatment-experienced subjects.