Rapid inference of antibiotic resistance and susceptibility by genomic neighbour typing

Molecular rapid diagnostic testing for bloodstream infections: Nanopore targeted sequencing with pathogen-specific primers

BACKGROUND

Nanopore sequencing, known for real-time analysis, shows promise for rapid clinical infection diagnosis but lacks effective assays for bloodstream infections (BSIs).

METHODS

We prospectively assessed the performance of a novel nanopore targeted sequencing (NTS) assay in identifying pathogens and predicting antibiotic resistance in BSIs, analyzing 387 blood samples from December 2021 to April 2023.

RESULTS

The positivity rate for NTS (69.5 %, 269/387) nearly matches that of metagenomic next-generation sequencing (mNGS) (74.7 %, 289/387; p = 0.128) and surpasses the positivity rate of conventional blood culture (BC) (33.9 %, 131/387; p < 0.01). Frequent pathogens detected by NTS included Klebsiella pneumoniae (n = 54), Pseudomonas aeruginosa (n = 36), Escherichia coli (n = 36), Enterococcus faecium(n = 30), Acinetobacter baumannii(n = 26), Staphylococcus aureus(n = 23), and Human cytomegalovirus (n = 37). Against a composite BSI diagnostic standard, NTS demonstrated a sensitivity and specificity of 84.0 % (95 % CI 79.5 %-87.7 %) and 90.1 % (95 % CI 81.7 %-88.5 %), respectively. The concordance between NTS and mNGS results (the percentage of total cases where both either detected BSI-related pathogens or were both negative) was 90.2 % (359/387), whereas the consistency between NTS and BC was only 60.2 % (233/387). In 80.6 % (50/62) of the samples with identical pathogens identified by both NTS tests and BCs, the genotypic resistance identified by NTS correlated with culture-confirmed phenotypic resistance. Using NTS, 95 % of samples can be tested and analyzed in approximately 7 h, allowing for early patient diagnosis.

CONCLUSIONS

NTS is rapid, sensitive, and efficient for detecting BSIs and drug-resistant genes, making it a potential preferred diagnostic tool for early infection identification in critically ill patients.

FARPA-based tube array coupled with quick DNA extraction enables ultra-fast bedside detection of antibiotic-resistant pathogens

Rapid and accurate detection of pathogens and antimicrobial-resistant (AMR) genes of the pathogens are crucial for the clinical diagnosis and effective treatment of infectious diseases. However, the time-consuming steps of conventional culture-based methods inhibit the precise and early application of anti-infection therapy. For the prompt treatment of pathogen-infected patients, we have proposed a novel tube array strategy based on our previously reported FARPA (FEN1-aided recombinase polymerase amplification) principle for the ultra-fast detection of antibiotic-resistant pathogens on site. The entire process from "sample to result" can be completed in 25 min by combining quick DNA extraction from a urine sample with FARPA to avoid the usually complicated DNA extraction step. Furthermore, a 36-tube array made from commercial 384-well titre plates was efficiently introduced to perform FARPA in a portable analyser, achieving an increase in the loading sample throughput (from several to several tens), which is quite suitable for the point-of-care testing (POCT) of multiple pathogens and multiple samples. Finally, we tested 92 urine samples to verify the performance of our proposed method. The sensitivities for the detection of E. coli, K. pneumoniae, E. faecium, and E. faecalis were 92.7%, 93.8%, 100% and 88.9%, respectively. The specificities for the detection of the four pathogens were 100%. Consequently, our rapid, low-cost and user-friendly POCT method holds great potential for guiding the rational use of antibiotics and reducing bacterial resistance.

Efficient and Robust Search of Microbial Genomes via Phylogenetic Compression

Comprehensive collections approaching millions of sequenced genomes have become central information sources in the life sciences. However, the rapid growth of these collections has made it effectively impossible to search these data using tools such as BLAST and its successors. Here, we present a technique called phylogenetic compression, which uses evolutionary history to guide compression and efficiently search large collections of microbial genomes using existing algorithms and data structures. We show that, when applied to modern diverse collections approaching millions of genomes, lossless phylogenetic compression improves the compression ratios of assemblies, de Bruijn graphs, and k -mer indexes by one to two orders of magnitude. Additionally, we develop a pipeline for a BLAST-like search over these phylogeny-compressed reference data, and demonstrate it can align genes, plasmids, or entire sequencing experiments against all sequenced bacteria until 2019 on ordinary desktop computers within a few hours. Phylogenetic compression has broad applications in computational biology and may provide a fundamental design principle for future genomics infrastructure.

Laboratory tools for the direct detection of bacterial respiratory infections and antimicrobial resistance: a scoping review

Rapid laboratory tests are urgently required to inform antimicrobial use in food animals. Our objective was to synthesize knowledge on the direct application of long-read metagenomic sequencing to respiratory samples to detect bacterial pathogens and antimicrobial resistance genes (ARGs) compared to PCR, loop-mediated isothermal amplification, and recombinase polymerase amplification. Our scoping review protocol followed the Joanna Briggs Institute and PRISMA Scoping Review reporting guidelines. Included studies reported on the direct application of these methods to respiratory samples from animals or humans to detect bacterial pathogens ±ARGs and included turnaround time (TAT) and analytical sensitivity. We excluded studies not reporting these or that were focused exclusively on bioinformatics. We identified 5,636 unique articles from 5 databases. Two-reviewer screening excluded 3,964, 788, and 784 articles at 3 levels, leaving 100 articles (19 animal and 81 human), of which only 7 studied long-read sequencing (only 1 in animals). Thirty-two studies investigated ARGs (only one in animals). Reported TATs ranged from minutes to 2 d; steps did not always include sample collection to results, and analytical sensitivity varied by study. Our review reveals a knowledge gap in research for the direct detection of bacterial respiratory pathogens and ARGs in animals using long-read metagenomic sequencing. There is an opportunity to harness the rapid development in this space to detect multiple pathogens and ARGs on a single sequencing run. Long-read metagenomic sequencing tools show potential to address the urgent need for research into rapid tests to support antimicrobial stewardship in food animal production.

Antimicrobial resistance prediction by clinical metagenomics in pediatric severe pneumonia patients

BACKGROUND

Antimicrobial resistance (AMR) is a major threat to children's health, particularly in respiratory infections. Accurate identification of pathogens and AMR is crucial for targeted antibiotic treatment. Metagenomic next-generation sequencing (mNGS) shows promise in directly detecting microorganisms and resistance genes in clinical samples. However, the accuracy of AMR prediction through mNGS testing needs further investigation for practical clinical decision-making.

METHODS

We aimed to evaluate the performance of mNGS in predicting AMR for severe pneumonia in pediatric patients. We conducted a retrospective analysis at a tertiary hospital from May 2022 to May 2023. Simultaneous mNGS and culture were performed on bronchoalveolar lavage fluid samples obtained from pediatric patients with severe pneumonia. By comparing the results of mNGS detection of microorganisms and antibiotic resistance genes with those of culture, sensitivity, specificity, positive predictive value, and negative predictive value were calculated.

RESULTS

mNGS detected bacterial in 71.7% cases (86/120), significantly higher than culture (58/120, 48.3%). Compared to culture, mNGS demonstrated a sensitivity of 96.6% and a specificity of 51.6% in detecting pathogenic microorganisms. Phenotypic susceptibility testing (PST) of 19 antibiotics revealed significant variations in antibiotics resistance rates among different bacteria. Sensitivity prediction of mNGS for carbapenem resistance was higher than penicillins and cephalosporin (67.74% vs. 28.57%, 46.15%), while specificity showed no significant difference (85.71%, 75.00%, 75.00%). mNGS also showed a high sensitivity of 94.74% in predicting carbapenem resistance in Acinetobacter baumannii.

CONCLUSIONS

mNGS exhibits variable predictive performance among different pathogens and antibiotics, indicating its potential as a supplementary tool to conventional PST. However, mNGS currently cannot replace conventional PST.

Assessing computational predictions of antimicrobial resistance phenotypes from microbial genomes

The advent of rapid whole-genome sequencing has created new opportunities for computational prediction of antimicrobial resistance (AMR) phenotypes from genomic data. Both rule-based and machine learning (ML) approaches have been explored for this task, but systematic benchmarking is still needed. Here, we evaluated four state-of-the-art ML methods (Kover, PhenotypeSeeker, Seq2Geno2Pheno and Aytan-Aktug), an ML baseline and the rule-based ResFinder by training and testing each of them across 78 species-antibiotic datasets, using a rigorous benchmarking workflow that integrates three evaluation approaches, each paired with three distinct sample splitting methods. Our analysis revealed considerable variation in the performance across techniques and datasets. Whereas ML methods generally excelled for closely related strains, ResFinder excelled for handling divergent genomes. Overall, Kover most frequently ranked top among the ML approaches, followed by PhenotypeSeeker and Seq2Geno2Pheno. AMR phenotypes for antibiotic classes such as macrolides and sulfonamides were predicted with the highest accuracies. The quality of predictions varied substantially across species-antibiotic combinations, particularly for beta-lactams; across species, resistance phenotyping of the beta-lactams compound, aztreonam, amoxicillin/clavulanic acid, cefoxitin, ceftazidime and piperacillin/tazobactam, alongside tetracyclines demonstrated more variable performance than the other benchmarked antibiotics. By organism, Campylobacter jejuni and Enterococcus faecium phenotypes were more robustly predicted than those of Escherichia coli, Staphylococcus aureus, Salmonella enterica, Neisseria gonorrhoeae, Klebsiella pneumoniae, Pseudomonas aeruginosa, Acinetobacter baumannii, Streptococcus pneumoniae and Mycobacterium tuberculosis. In addition, our study provides software recommendations for each species-antibiotic combination. It furthermore highlights the need for optimization for robust clinical applications, particularly for strains that diverge substantially from those used for training.

Optimising machine learning prediction of minimum inhibitory concentrations in Klebsiella pneumoniae

Minimum Inhibitory Concentrations (MICs) are the gold standard for quantitatively measuring antibiotic resistance. However, lab-based MIC determination can be time-consuming and suffers from low reproducibility, and interpretation as sensitive or resistant relies on guidelines which change over time. Genome sequencing and machine learning promise to allow in silico MIC prediction as an alternative approach which overcomes some of these difficulties, albeit the interpretation of MIC is still needed. Nevertheless, precisely how we should handle MIC data when dealing with predictive models remains unclear, since they are measured semi-quantitatively, with varying resolution, and are typically also left- and right-censored within varying ranges. We therefore investigated genome-based prediction of MICs in the pathogen Klebsiella pneumoniae using 4367 genomes with both simulated semi-quantitative traits and real MICs. As we were focused on clinical interpretation, we used interpretable rather than black-box machine learning models, namely, Elastic Net, Random Forests, and linear mixed models. Simulated traits were generated accounting for oligogenic, polygenic, and homoplastic genetic effects with different levels of heritability. Then we assessed how model prediction accuracy was affected when MICs were framed as regression and classification. Our results showed that treating the MICs differently depending on the number of concentration levels of antibiotic available was the most promising learning strategy. Specifically, to optimise both prediction accuracy and inference of the correct causal variants, we recommend considering the MICs as continuous and framing the learning problem as a regression when the number of observed antibiotic concentration levels is large, whereas with a smaller number of concentration levels they should be treated as a categorical variable and the learning problem should be framed as a classification. Our findings also underline how predictive models can be improved when prior biological knowledge is taken into account, due to the varying genetic architecture of each antibiotic resistance trait. Finally, we emphasise that incrementing the population database is pivotal for the future clinical implementation of these models to support routine machine-learning based diagnostics.

Applicability of Bronchoalveolar Lavage Fluid and Plasma Metagenomic Next-Generation Sequencing Assays in the Diagnosis of Pneumonia

Background

Metagenomic next-generation sequencing (mNGS) provides innovative solutions for predicting complex infections. A comprehensive understanding of its strengths and limitations in real-world clinical settings is necessary to ensure that it is not overused or misinterpreted.

Methods

Two hundred nine cases with suspected pneumonia were recruited to compare the capabilities of 2 available mNGS assays (bronchoalveolar lavage fluid [BALF] mNGS and plasma mNGS) to identify pneumonia-associated DNA/RNA pathogens and predict antibiotic resistance.

Results

Compared to clinical diagnosis, BALF mNGS demonstrated a high positive percent agreement (95.3%) but a low negative percent agreement (63.1%). Plasma mNGS revealed a low proportion of true negatives (30%) in predicting pulmonary infection. BALF mNGS independently diagnosed 65.6% (61/93) of coinfections and had a remarkable advantage in detecting caustic, rare, or atypical pathogens. Pathogens susceptible to invasive infection or bloodstream transmission, such as Aspergillus spp, Rhizopus spp, Chlamydia psittaci, and human herpesviruses, are prone to be detected by plasma mNGS. BALF mNGS tests provided a positive impact on the diagnosis and treatment of 128 (61.2%) patients. Plasma mNGS, on the other hand, turned out to be more suitable for diagnosing patients who received mechanical ventilation, developed severe pneumonia, or developed sepsis (all P < .01). BALF mNGS was able to identify resistance genes that matched the phenotypic resistance of 69.4% (25/36) of multidrug-resistant pathogens.

Conclusions

Our data reveal new insights into the advantages and disadvantages of 2 different sequencing modalities in pathogen identification and antibiotic resistance prediction for patients with suspected pneumonia.

A dual-process of targeted and unbiased Nanopore sequencing enables accurate and rapid diagnosis of lower respiratory infections

BACKGROUND

Nanopore metagenomics has been used for infectious disease diagnosis for bacterial pathogens. However, this technology currently lacks comprehensive performance studies in clinical settings for simultaneous detection of bacteria, fungi, and viruses.

METHODS

We developed a dual-process of Nanopore sequencing for one sample, with unbiased metagenomics in Meta process and target enrichment in Panel process (Nanopore Meta-Panel process, NanoMP) and prospectively enrolled 450 respiratory specimens from multiple centers. The filter system of pathogen detection was established with machine learning and receiver operator characteristic (ROC) curve to optimize the detection accuracy based on orthogonal test of 21 species. Antimicrobial resistance (AMR) genes were identified based on the Comprehensive Antibiotic Resistance Database (CARD) and single-nucleotide polymorphism matrix.

FINDINGS

Our approach showed high sensitivity in Meta process, with 82.9%, 88.7%, and 75.0% for bacteria, fungi (except Aspergillus), and Mycobacterium tuberculosis groups, respectively. Moreover, target amplification improved the sensitivity of virus (>80.0% vs. 39.4%) and Aspergillus (81.8% vs. 42.3%) groups in Panel process compared with Meta process. Overall, NanoMP achieved 80.2% sensitivity and 98.8% specificity compared with the composite reference standard, and we were able to accurately detect AMR genes including blaKPC-2, blaOXA-23 and mecA and distinguish their parent organisms in patients with mixed infections.

INTERPRETATION

We combined metagenomic and enriched Nanopore sequencing for one sample in parallel. Our NanoMP approach simultaneously covered bacteria, viruses and fungi in respiratory specimens and demonstrated good diagnostic performance in real clinical settings.

FUNDING

National Key Research and Development Program of China and National Natural Science Foundation of China.

Time Required for Nanopore Whole-Genome Sequencing of Neisseria gonorrhoeae for Identification of Phylogenetic Relationships

BACKGROUND

Antimicrobial resistance (AMR) in Neisseria gonorrhoeae is a global health challenge. Limitations to AMR surveillance reporting, alongside reduction in culture-based susceptibility testing, has resulted in a need for rapid diagnostics and strain detection. We investigated Nanopore sequencing time, and depth, to accurately identify closely related N. gonorrhoeae isolates, compared to Illumina sequencing.

METHODS

N. gonorrhoeae strains collected from a London sexual health clinic were cultured and sequenced with MiSeq and MinION sequencing platforms. Accuracy was determined by comparing variant calls at 68 nucleotide positions (37 resistance-associated markers). Accuracy at varying MinION sequencing depths was determined through retrospective time-stamped read analysis.

RESULTS

Of 22 MinION-MiSeq pairs reaching sufficient sequencing depth, agreement of variant call positions passing quality control criteria was 185/185 (100%; 95% confidence interval [CI], 98.0%-100.0%), 502/503 (99.8%; 95% CI, 98.9%-99.9%), and 564/565 (99.8%; 95% CI, 99.0%-100.0%) at 10x, 30x, and 40x MinION depth, respectively. Isolates identified as closely related by MiSeq, within one yearly evolutionary distance of ≤5 single nucleotide polymorphisms, were accurately identified via MinION.

CONCLUSIONS

Nanopore sequencing shows utility as a rapid surveillance tool, identifying closely related N. gonorrhoeae strains, with just 10x sequencing depth, taking a median time of 29 minutes. This highlights its potential for tracking local transmission and AMR markers.

Urinary tract infections: a review of the current diagnostics landscape

Urinary tract infections are the most common bacterial infections worldwide. Infections can range from mild, recurrent (rUTI) to complicated (cUTIs), and are predominantly caused by uropathogenic Escherichia coli (UPEC). Antibiotic therapy is important to tackle infection; however, with the continued emergence of antibiotic resistance there is an urgent need to monitor the use of effective antibiotics through better stewardship measures. Currently, clinical diagnosis of UTIs relies on empiric methods supported by laboratory testing including cellular analysis (of both human and bacterial cells), dipstick analysis and phenotypic culture. Therefore, development of novel, sensitive and specific diagnostics is an important means to rationalise antibiotic therapy in patients. This review discusses the current diagnostic landscape and highlights promising novel diagnostic technologies in development that could aid in treatment and management of antibiotic-resistant UTIs.

Predicting variable gene content in Escherichia coli using conserved genes

Having the ability to predict the protein-encoding gene content of an incomplete genome or metagenome-assembled genome is important for a variety of bioinformatic tasks. In this study, as a proof of concept, we built machine learning classifiers for predicting variable gene content in Escherichia coli genomes using only the nucleotide k-mers from a set of 100 conserved genes as features. Protein families were used to define orthologs, and a single classifier was built for predicting the presence or absence of each protein family occurring in 10%-90% of all E. coli genomes. The resulting set of 3,259 extreme gradient boosting classifiers had a per-genome average macro F1 score of 0.944 [0.943-0.945, 95% CI]. We show that the F1 scores are stable across multi-locus sequence types and that the trend can be recapitulated by sampling a smaller number of core genes or diverse input genomes. Surprisingly, the presence or absence of poorly annotated proteins, including "hypothetical proteins" was accurately predicted (F1 = 0.902 [0.898-0.906, 95% CI]). Models for proteins with horizontal gene transfer-related functions had slightly lower F1 scores but were still accurate (F1s = 0.895, 0.872, 0.824, and 0.841 for transposon, phage, plasmid, and antimicrobial resistance-related functions, respectively). Finally, using a holdout set of 419 diverse E. coli genomes that were isolated from freshwater environmental sources, we observed an average per-genome F1 score of 0.880 [0.876-0.883, 95% CI], demonstrating the extensibility of the models. Overall, this study provides a framework for predicting variable gene content using a limited amount of input sequence data. IMPORTANCE Having the ability to predict the protein-encoding gene content of a genome is important for assessing genome quality, binning genomes from shotgun metagenomic assemblies, and assessing risk due to the presence of antimicrobial resistance and other virulence genes. In this study, we built a set of binary classifiers for predicting the presence or absence of variable genes occurring in 10%-90% of all publicly available E. coli genomes. Overall, the results show that a large portion of the E. coli variable gene content can be predicted with high accuracy, including genes with functions relating to horizontal gene transfer. This study offers a strategy for predicting gene content using limited input sequence data.

Real-time genomics for One Health

The ongoing degradation of natural systems and other environmental changes has put our society at a crossroad with respect to our future relationship with our planet. While the concept of One Health describes how human health is inextricably linked with environmental health, many of these complex interdependencies are still not well-understood. Here, we describe how the advent of real-time genomic analyses can benefit One Health and how it can enable timely, in-depth ecosystem health assessments. We introduce nanopore sequencing as the only disruptive technology that currently allows for real-time genomic analyses and that is already being used worldwide to improve the accessibility and versatility of genomic sequencing. We showcase real-time genomic studies on zoonotic disease, food security, environmental microbiome, emerging pathogens, and their antimicrobial resistances, and on environmental health itself - from genomic resource creation for wildlife conservation to the monitoring of biodiversity, invasive species, and wildlife trafficking. We stress why equitable access to real-time genomics in the context of One Health will be paramount and discuss related practical, legal, and ethical limitations.

Towards monitoring of antimicrobial resistance in the environment: For what reasons, how to implement it, and what are the data needs?

Antimicrobial resistance (AMR) is a global threat to human and animal health and well-being. To understand AMR dynamics, it is important to monitor resistant bacteria and resistance genes in all relevant settings. However, while monitoring of AMR has been implemented in clinical and veterinary settings, comprehensive monitoring of AMR in the environment is almost completely lacking. Yet, the environmental dimension of AMR is critical for understanding the dissemination routes and selection of resistant microorganisms, as well as the human health risks related to environmental AMR. Here, we outline important knowledge gaps that impede implementation of environmental AMR monitoring. These include lack of knowledge of the 'normal' background levels of environmental AMR, definition of high-risk environments for transmission, and a poor understanding of the concentrations of antibiotics and other chemical agents that promote resistance selection. Furthermore, there is a lack of methods to detect resistance genes that are not already circulating among pathogens. We conclude that these knowledge gaps need to be addressed before routine monitoring for AMR in the environment can be implemented on a large scale. Yet, AMR monitoring data bridging different sectors is needed in order to fill these knowledge gaps, which means that some level of national, regional and global AMR surveillance in the environment must happen even without all scientific questions answered. With the possibilities opened up by rapidly advancing technologies, it is time to fill these knowledge gaps. Doing so will allow for specific actions against environmental AMR development and spread to pathogens and thereby safeguard the health and wellbeing of humans and animals.

The diagnostic utility of nanopore targeted sequencing in suspected endophthalmitis

PURPOSE

This paper aimed to assess the diagnostic utility of a newly developed gene-based technology-nanopore targeted sequencing (NTS) in suspected endophthalmitis patients.

METHODS

This retrospective study included 43 patients (44 eyes) with suspected endophthalmitis. NTS was applied along with microbiological culture to detect unknown pathogens in intraocular fluid samples. The diagnostic utility of NTS was mainly evaluated from three aspects, including the positivity rate of bacterial/fungal presence, diagnostic turnaround time and the frequency of change in treatment based on etiology test results. Non-parametric, two-sided Wilcoxon rank sum test, the McNemar's test and the kappa statistic were used for statistical comparisons.

RESULTS

NTS showed significant advantages over traditional culture in positivity rates and diagnostic time (P < 0.001, kappa = 0.082; Z = -5.805, P < 0. 001). As regards antibiotic strategy, 17 patients (39.53%) and 5 patients (11.63%) underwent medication change following NTS and culture results respectively (P < 0.001, kappa = 0.335). With reasonable use of antibiotic and surgical intervention, most patients responded favorably, judged by significantly improved visual acuity (Z = -4.249, P < 0.001). The mean duration of hospitalization was 8.49 ± 2.45 days (range, 1-16 days).

CONCLUSION

The high efficiency feature of NTS in pathogen detection renders it a valuable supplementary to traditional culture. Additionally, it has facilitated patients' management for the early and precise diagnosis of endophthalmitis.

Rapid metagenomic sequencing for diagnosis and antimicrobial sensitivity prediction of canine bacterial infections

Antimicrobial resistance is a major threat to human and animal health. There is an urgent need to ensure that antimicrobials are used appropriately to limit the emergence and impact of resistance. In the human and veterinary healthcare setting, traditional culture and antimicrobial sensitivity testing typically requires 48-72 h to identify appropriate antibiotics for treatment. In the meantime, broad-spectrum antimicrobials are often used, which may be ineffective or impact non-target commensal bacteria. Here, we present a rapid, culture-free, diagnostics pipeline, involving metagenomic nanopore sequencing directly from clinical urine and skin samples of dogs. We have planned this pipeline to be versatile and easily implementable in a clinical setting, with the potential for future adaptation to different sample types and animals. Using our approach, we can identify the bacterial pathogen present within 5 h, in some cases detecting species which are difficult to culture. For urine samples, we can predict antibiotic sensitivity with up to 95 % accuracy. Skin swabs usually have lower bacterial abundance and higher host DNA, confounding antibiotic sensitivity prediction; an additional host depletion step will likely be required during the processing of these, and other types of samples with high levels of host cell contamination. In summary, our pipeline represents an important step towards the design of individually tailored veterinary treatment plans on the same day as presentation, facilitating the effective use of antibiotics and promoting better antimicrobial stewardship.

Genealogical inference and more flexible sequence clustering using iterative-PopPUNK

Bacterial genome data are accumulating at an unprecedented speed due to the routine use of sequencing in clinical diagnoses, public health surveillance, and population genetics studies. Genealogical reconstruction is fundamental to many of these uses; however, inferring genealogy from large-scale genome data sets quickly, accurately, and flexibly is still a challenge. Here, we extend an alignment- and annotation-free method, PopPUNK, to increase its flexibility and interpretability across data sets. Our method, iterative-PopPUNK, rapidly produces multiple consistent cluster assignments across a range of sequence identities. By constructing a partially resolved genealogical tree with respect to these clusters, users can select a resolution most appropriate for their needs. We showed the accuracy of clusters at all levels of similarity and genealogical inference of iterative-PopPUNK based on simulated data and obtained phylogenetically concordant results in real data sets from seven bacterial species. Using two example sets of Escherichia/Shigella and Vibrio parahaemolyticus genomes, we show that iterative-PopPUNK can achieve cluster resolutions ranging from phylogroup down to sequence typing (ST). The iterative-PopPUNK algorithm is implemented in the "PopPUNK_iterate" program, available as part of the PopPUNK package.

Real-Time Nanopore Q20+ Sequencing Enables Extremely Fast and Accurate Core Genome MLST Typing and Democratizes Access to High-Resolution Bacterial Pathogen Surveillance

Next-generation whole-genome sequencing is essential for high-resolution surveillance of bacterial pathogens, for example, during outbreak investigations or for source tracking and escape variant analysis. However, current global sequencing and bioinformatic bottlenecks and a long time to result with standard technologies demand new approaches. In this study, we investigated whether novel nanopore Q20+ long-read chemistry enables standardized and easily accessible high-resolution typing combined with core genome multilocus sequence typing (cgMLST). We set high requirements for discriminatory power by using the slowly evolving bacterium Bordetella pertussis as a model pathogen. Our results show that the increased raw read accuracy enables the description of epidemiological scenarios and phylogenetic linkages at the level of gold-standard short reads. The same was true for our variant analysis of vaccine antigens, resistance genes, and virulence factors, demonstrating that nanopore sequencing is a legitimate competitor in the area of next-generation sequencing (NGS)-based high-resolution bacterial typing. Furthermore, we evaluated the parameters for the fastest possible analysis of the data. By combining the optimized processing pipeline with real-time basecalling, we established a workflow that allows for highly accurate and extremely fast high-resolution typing of bacterial pathogens while sequencing is still in progress. Along with advantages such as low costs and portability, the approach suggested here might democratize modern bacterial typing, enabling more efficient infection control globally.

A case for investment in clinical metagenomics in low-income and middle-income countries

Clinical metagenomics is the diagnostic approach with the broadest capacity to detect both known and novel pathogens. Clinical metagenomics is costly to run and requires infrastructure, but the use of next-generation sequencing for SARS-CoV-2 molecular epidemiology in low-income and middle-income countries (LMICs) offers an opportunity to direct this infrastructure to the establishment of clinical metagenomics programmes. Local implementation of clinical metagenomics is important to create relevant systems and evaluate cost-effective methodologies for its use, as well as to ensure that reference databases and result interpretation tools are appropriate to local epidemiology. Rational implementation, based on the needs of LMICs and the available resources, could ultimately improve individual patient care in instances in which available diagnostics are inadequate and supplement emerging infectious disease surveillance systems to ensure the next pandemic pathogen is quickly identified.

Development of an Amplicon Nanopore Sequencing Strategy for Detection of Mutations Conferring Intermediate Resistance to Vancomycin in Staphylococcus aureus Strains

Staphylococcus aureus is a major cause of bacteremia and other hospital-acquired infections. The cell-wall active antibiotic vancomycin is commonly used to treat both methicillin-resistant (MRSA) and sensitive (MSSA) infections. Vancomycin intermediate S. aureus (VISA) variants can arise through de novo mutations. Here, we performed pilot experiments to develop a combined PCR/long-read sequencing-based method for detection of previously known VISA-causing mutations. Primers were designed to generate 10 amplicons covering 16 genes associated with the VISA phenotype. We sequenced amplicon pools as long reads with Oxford Nanopore adapter ligation on Flongle flow cells. We then detected mutations by mapping reads against a parental consensus or known reference sequence and comparing called variants against a database of known VISA mutations from laboratory selection. Each amplicon in the pool was sequenced to high (>1,000×) coverage, and no relationship was found between amplicon length and coverage. We also were able to detect the causative mutation (walK 646C>G) in a VISA mutant derived from the USA300 strain (N384-3 from parental strain N384). Mixing mutant (N384-3) and parental (N384) DNA at various ratios from 0 to 1 mutant suggested a mutation detection threshold of the average minor allele frequency (6.5%) at 95% confidence (two standard errors above mean mutation frequency). The study lays the groundwork for direct S. aureus antibiotic resistance genotype inference using rapid nanopore sequencing from clinical samples. IMPORTANCE Bacteremia mortality is known to increase rapidly with time after infection, making rapid diagnostics and treatment necessary. Successful treatment depends on correct administration of antibiotics based on knowledge of strain antibiotic susceptibility. Staphylococcus aureus is a major causative agent of bacteremia that is also commonly antibiotic resistant. In this work, we develop a method to accelerate detection of a complex, polygenic antibiotic resistance phenotype in S. aureus, vancomycin-intermediate resistance (VISA), through long-read genomic sequencing of amplicons representing genes most commonly mutated in VISA selection. This method both rapidly identifies VISA genotypes and incorporates the most comprehensive database of VISA genetic determinants known to date.

Metagenomic Antimicrobial Susceptibility Testing from Simulated Native Patient Samples

Genomic antimicrobial susceptibility testing (AST) has been shown to be accurate for many pathogens and antimicrobials. However, these methods have not been systematically evaluated for clinical metagenomic data. We investigate the performance of in-silico AST from clinical metagenomes (MG-AST). Using isolate sequencing data from a multi-center study on antimicrobial resistance (AMR) as well as shotgun-sequenced septic urine samples, we simulate over 2000 complicated urinary tract infection (cUTI) metagenomes with known resistance phenotype to 5 antimicrobials. Applying rule-based and machine learning-based genomic AST classifiers, we explore the impact of sequencing depth and technology, metagenome complexity, and bioinformatics processing approaches on AST accuracy. By using an optimized metagenomics assembly and binning workflow, MG-AST achieved balanced accuracy within 5.1% of isolate-derived genomic AST. For poly-microbial infections, taxonomic sample complexity and relatedness of taxa in the sample is a key factor influencing metagenomic binning and downstream MG-AST accuracy. We show that the reassignment of putative plasmid contigs by their predicted host range and investigation of whole resistome capabilities improved MG-AST performance on poly-microbial samples. We further demonstrate that machine learning-based methods enable MG-AST with superior accuracy compared to rule-based approaches on simulated native patient samples.

Machine learning models for Neisseria gonorrhoeae antimicrobial susceptibility tests

Neisseria gonorrhoeae is an urgent public health threat due to the emergence of antibiotic resistance. As most isolates in the United States are susceptible to at least one antibiotic, rapid molecular antimicrobial susceptibility tests (ASTs) would offer the opportunity to tailor antibiotic therapy, thereby expanding treatment options. With genome sequence and antibiotic resistance phenotype data for nearly 20,000 clinical N. gonorrhoeae isolates now available, there is an opportunity to use statistical methods to develop sequence-based diagnostics that predict antibiotic susceptibility from genotype. N. gonorrhoeae, therefore, provides a useful example illustrating how to apply machine learning models to aid in the design of sequence-based ASTs. We present an overview of this framework, which begins with establishing the assay technology, the performance criteria, the population in which the diagnostic will be used, and the clinical goals, and extends to the choices that must be made to arrive at a set of features with the desired properties for predicting susceptibility phenotype from genotype. While we focus on the example of N. gonorrhoeae, the framework generalizes to other organisms for which large-scale genotype and antibiotic resistance data can be combined to aid in diagnostics development.

Performance Characteristics of Next-Generation Sequencing for the Detection of Antimicrobial Resistance Determinants in Escherichia coli Genomes and Metagenomes

Short-read sequencing can provide detection of multiple genomic determinants of antimicrobial resistance from single bacterial genomes and metagenomic samples. Despite its increasing application in human, animal, and environmental microbiology, including human clinical trials, the performance of short-read Illumina sequencing for antimicrobial resistance gene (ARG) detection, including resistance-conferring single nucleotide polymorphisms (SNPs), has not been systematically characterized. Using paired-end 2 × 150 bp (base pair) Illumina sequencing and an assembly-based method for ARG prediction, we determined sensitivity, positive predictive value (PPV), and sequencing depths required for ARG detection in an Escherichia coli isolate of sequence type (ST) 38 spiked into a synthetic microbial community at varying abundances. Approximately 300,000 reads or 15× genome coverage was sufficient to detect ARGs in E. coli ST38, with comparable sensitivity and PPV to ~100× genome coverage. Using metagenome assembly of mixed microbial communities, ARG detection at E. coli relative abundances of 1% would require assembly of approximately 30 million reads to achieve 15× target coverage. The minimum sequencing depths were validated using public data sets of 948 E. coli genomes and 10 metagenomic rectal swab samples. A read-based approach using k-mer alignment (KMA) for ARG prediction did not substantially improve minimum sequencing depths for ARG detection compared to assembly of the E. coli ST38 genome or the combined metagenomic samples. Analysis of sequencing depths from recent studies assessing ARG content in metagenomic samples demonstrated that sequencing depths had a median estimated detection frequency of 84% (interquartile range: 30%-92%) for a relative abundance of 1%. IMPORTANCE Systematically determining Illumina sequencing performance characteristics for detection of ARGs in metagenomic samples is essential to inform study design and appraisal of human, animal, and environmental metagenomic antimicrobial resistance studies. In this study, we quantified the performance characteristics of ARG detection in E. coli genomes and metagenomes and established a benchmark of ~15× coverage for ARG detection for E. coli in metagenomes. We demonstrate that for low relative abundances, sequencing depths of ~30 million reads or more may be required for adequate sensitivity for many applications.

Multiplexed and Rapid AST for Escherichia coli Infection by Simultaneously Pyrosequencing Multiple Barcodes Each Specific to an Antibiotic Exposed to a Sample

Antimicrobial susceptibility testing (AST) is an effective way to guide antibiotic selection. However, conventional culture-based phenotypic AST is time-consuming. The key point to shorten the test is to quantify the small change in the bacterial number after the antibiotic exposure. To achieve rapid AST, we proposed a combination of multiplexed PCR with barcoded pyrosequencing to significantly shorten the time for antibiotic exposure. First, bacteria exposed to each antibiotic were labeled with a unique barcode. Then, the pool of the barcoded products was amplified by PCR with a universal primer pair. Finally, barcodes in the amplicons were individually and quantitatively decoded by pyrosequencing. As pyrosequencing is able to discriminate as low as 5% variation in target concentrations, as short as 7.5 min was enough for cultivation to detect the susceptibility of Escherichia coli to an antibiotic. The barcodes enable more than six kinds of drugs or six kinds of concentrations of a drug to be tested at a time. The susceptibility of 6 antibiotics to 43 E. coli-positive samples from 482 clinical urine samples showed a consistency of 99.3% for drug-resistant samples and of 95.7% for drug-sensitive samples in comparison with the conventional method. In addition, the minimum inhibitory concentration (MIC) of 29 E. coli samples was successfully measured. The proposed AST is dye free (pyrosequencing), multiplexed (six antibiotics), fast (a half-working day for reporting the results), and able to detect the MIC, thus having a great potential for clinical use in quick antibiotic selection.

Systems-Based Approach for Optimization of Assembly-Free Bacterial MLST Mapping

Epidemiological surveillance of bacterial pathogens requires real-time data analysis with a fast turnaround, while aiming at generating two main outcomes: (1) species-level identification and (2) variant mapping at different levels of genotypic resolution for population-based tracking and surveillance, in addition to predicting traits such as antimicrobial resistance (AMR). Multi-locus sequence typing (MLST) aids this process by identifying sequence types (ST) based on seven ubiquitous genome-scattered loci. In this paper, we selected one assembly-dependent and one assembly-free method for ST mapping and applied them with the default settings and ST schemes they are distributed with, and systematically assessed their accuracy and scalability across a wide array of phylogenetically divergent Public Health-relevant bacterial pathogens with available MLST databases. Our data show that the optimal k-mer length for stringMLST is species-specific and that genome-intrinsic and -extrinsic features can affect the performance and accuracy of the program. Although suitable parameters could be identified for most organisms, there were instances where this program may not be directly deployable in its current format. Next, we integrated stringMLST into our freely available and scalable hierarchical-based population genomics platform, ProkEvo, and further demonstrated how the implementation facilitates automated, reproducible bacterial population analysis.

Clinical Metagenomic Sequencing for Species Identification and Antimicrobial Resistance Prediction in Orthopedic Device Infection

Diagnosis of orthopedic device-related infection is challenging, and causative pathogens may be difficult to culture. Metagenomic sequencing can diagnose infections without culture, but attempts to detect antimicrobial resistance (AMR) determinants using metagenomic data have been less successful. Human DNA depletion may maximize the amount of microbial DNA sequence data available for analysis. Human DNA depletion by saponin was tested in 115 sonication fluid samples generated following revision arthroplasty surgery, comprising 67 where pathogens were detected by culture and 48 culture-negative samples. Metagenomic sequencing was performed on the Oxford Nanopore Technologies GridION platform. Filtering thresholds for detection of true species versus contamination or taxonomic misclassification were determined. Mobile and chromosomal genetic AMR determinants were identified in Staphylococcus aureus-positive samples. Of 114 samples generating sequence data, species-level positive percent agreement between metagenomic sequencing and culture was 50/65 (77%; 95% confidence interval [CI], 65 to 86%) and negative percent agreement was 103/114 (90%; 95% CI, 83 to 95%). Saponin treatment reduced the proportion of human bases sequenced in comparison to 5-μm filtration from a median (interquartile range [IQR]) of 98.1% (87.0% to 99.9%) to 11.9% (0.4% to 67.0%), improving reference genome coverage at a 10-fold depth from 18.7% (0.30% to 85.7%) to 84.3% (12.9% to 93.8%). Metagenomic sequencing predicted 13/15 (87%) resistant and 74/74 (100%) susceptible phenotypes where sufficient data were available for analysis. Metagenomic nanopore sequencing coupled with human DNA depletion has the potential to detect AMR in addition to species detection in orthopedic device-related infection. Further work is required to develop pathogen-agnostic human DNA depletion methods, improving AMR determinant detection and allowing its application to other infection types.

K-mer based prediction of Clostridioides difficile relatedness and ribotypes

Comparative analysis of Clostridioides difficile whole-genome sequencing (WGS) data enables fine scaled investigation of transmission and is increasingly becoming part of routine surveillance. However, these analyses are constrained by the computational requirements of the large volumes of data involved. By decomposing WGS reads or assemblies into k-mers and using the dimensionality reduction technique MinHash, it is possible to rapidly approximate genomic distances without alignment. Here we assessed the performance of MinHash, as implemented by sourmash, in predicting single nucleotide differences between genomes (SNPs) and C. difficile ribotypes (RTs). For a set of 1905 diverse C. difficile genomes (differing by 0–168 519 SNPs), using sourmash to screen for closely related genomes, at a sensitivity of 100 % for pairs ≤10 SNPs, sourmash reduced the number of pairs from 1 813 560 overall to 161 934, i.e. by 91 %, with a positive predictive value of 32 % to correctly identify pairs ≤10 SNPs (maximum SNP distance 4144). At a sensitivity of 95 %, pairs were reduced by 94 % to 108 266 and PPV increased to 45 % (maximum SNP distance 1009). Increasing the MinHash sketch size above 2000 produced minimal performance improvement. We also explored a MinHash similarity-based ribotype prediction method. Genomes with known ribotypes (n=3937) were split into a training set (2937) and test set (1000) randomly. The training set was used to construct a sourmash index against which genomes from the test set were compared. If the closest five genomes in the index had the same ribotype this was taken to predict the searched genome’s ribotype. Using our MinHash ribotype index, predicted ribotypes were correct in 780/1000 (78 %) genomes, incorrect in 20 (2 %), and indeterminant in 200 (20 %). Relaxing the classifier to 4/5 closest matches with the same RT improved the correct predictions to 87 %. Using MinHash it is possible to subsample C. difficile genome k-mer hashes and use them to approximate small genomic differences within minutes, significantly reducing the search space for further analysis.

Genomic heterogeneity underlies multidrug resistance in Pseudomonas aeruginosa: A population-level analysis beyond susceptibility testing

BACKGROUND

Pseudomonas aeruginosa is a persistent and difficult-to-treat pathogen in many patients, especially those with Cystic Fibrosis (CF). Herein, we describe a longitudinal analysis of a series of multidrug resistant (MDR) P. aeruginosa isolates recovered in a 17-month period, from a young female CF patient who underwent double lung transplantation. Our goal was to understand the genetic basis of the observed resistance phenotypes, establish the genomic population diversity, and define the nature of sequence evolution over time.

METHODS

Twenty-two sequential P. aeruginosa isolates were obtained within a 17-month period, before and after a double-lung transplant. At the end of the study period, antimicrobial susceptibility testing, whole genome sequencing (WGS), phylogenetic analyses and RNAseq were performed in order to understand the genetic basis of the observed resistance phenotypes, establish the genomic population diversity, and define the nature of sequence changes over time.

RESULTS

The majority of isolates were resistant to almost all tested antibiotics. A phylogenetic reconstruction revealed 3 major clades representing a genotypically and phenotypically heterogeneous population. The pattern of mutation accumulation and variation of gene expression suggested that a group of closely related strains was present in the patient prior to transplantation and continued to change throughout the course of treatment. A trend toward accumulation of mutations over time was observed. Different mutations in the DNA mismatch repair gene mutL consistent with a hypermutator phenotype were observed in two clades. RNAseq performed on 12 representative isolates revealed substantial differences in the expression of genes associated with antibiotic resistance and virulence traits.

CONCLUSIONS

The overwhelming current practice in the clinical laboratories setting relies on obtaining a pure culture and reporting the antibiogram from a few isolated colonies to inform therapy decisions. Our analyses revealed significant underlying genomic heterogeneity and unpredictable evolutionary patterns that were independent of prior antibiotic treatment, highlighting the need for comprehensive sampling and population-level analysis when gathering microbiological data in the context of CF P. aeruginosa chronic infection. Our findings challenge the applicability of antimicrobial stewardship programs based on single-isolate resistance profiles for the selection of antibiotic regimens in chronic infections such as CF.

A Practical Approach for Predicting Antimicrobial Phenotype Resistance in Staphylococcus aureus Through Machine Learning Analysis of Genome Data

With the reduction in sequencing price and acceleration of sequencing speed, it is particularly important to directly link the genotype and phenotype of bacteria. Here, we firstly predicted the minimum inhibitory concentrations of ten antimicrobial agents for Staphylococcus aureus using 466 isolates by directly extracting k-mer from whole genome sequencing data combined with three machine learning algorithms: random forest, support vector machine, and XGBoost. Considering one two-fold dilution, the essential agreement and the category agreement could reach >85% and >90% for most antimicrobial agents. For clindamycin, cefoxitin and trimethoprim-sulfamethoxazole, the essential agreement and the category agreement could reach >91% and >93%, providing important information for clinical treatment. The successful prediction of cefoxitin resistance showed that the model could identify methicillin-resistant S. aureus. The results suggest that small datasets available in large hospitals could bypass the existing basic research and known antimicrobial resistance genes and accurately predict the bacterial phenotype.

Rapid expansion and extinction of antibiotic resistance mutations during treatment of acute bacterial respiratory infections

Acute bacterial infections are often treated empirically, with the choice of antibiotic therapy updated during treatment. The effects of such rapid antibiotic switching on the evolution of antibiotic resistance in individual patients are poorly understood. Here we find that low-frequency antibiotic resistance mutations emerge, contract, and even go to extinction within days of changes in therapy. We analyzed Pseudomonas aeruginosa populations in sputum samples collected serially from 7 mechanically ventilated patients at the onset of respiratory infection. Combining short- and long-read sequencing and resistance phenotyping of 420 isolates revealed that while new infections are near-clonal, reflecting a recent colonization bottleneck, resistance mutations could emerge at low frequencies within days of therapy. We then measured the in vivo frequencies of select resistance mutations in intact sputum samples with resistance-targeted deep amplicon sequencing (RETRA-Seq), which revealed that rare resistance mutations not detected by clinically used culture-based methods can increase by nearly 40-fold over 5–12 days in response to antibiotic changes. Conversely, mutations conferring resistance to antibiotics not administered diminish and even go to extinction. Our results underscore how therapy choice shapes the dynamics of low-frequency resistance mutations at short time scales, and the findings provide a possibility for driving resistance mutations to extinction during early stages of infection by designing patient-specific antibiotic cycling strategies informed by deep genomic surveillance.

It remains unclear how rapid antibiotic switching affects the evolution of antibiotic resistance in individual patients. Here, Chung et al. combine short- and long-read sequencing and resistance phenotyping of 420 serial isolates of Pseudomonas aeruginosa collected from the onset of respiratory infection, and show that rare resistance mutations can increase by nearly 40-fold over 5–12 days in response to antibiotic changes, while mutations conferring resistance to antibiotics not administered diminish and even go to extinction.

Surveillance and epidemiology of syphilis, gonorrhoea and chlamydia in the non-European Union countries of the World Health Organization European Region, 2015 to 2020

Background

Epidemics of sexually transmitted infections (STI) are a major public health challenge in the World Health Organization (WHO) European Region.

Aim

We aimed to provide an overview of case reporting and other surveillance data for syphilis, gonorrhoea and chlamydia for the non-European Union (EU)/European Economic Area (EEA) countries of the Centre and East part of the WHO European Region as per classification used by the WHO Regional Office for Europe (WHO/Europe) and the European Centre for Disease Prevention and Control.

Methods

Data were provided by the surveillance agencies of the Member States for the period 2015 to 2019 through the WHO/Europe Communicable Diseases Annual Reporting Form. We analysed reported cases, explored data reported to the WHO Gonococcal Antimicrobial Surveillance Programme (GASP) and performed a review of publications on antimicrobial resistance (AMR) in gonorrhoea in the period 2015 to 2020 using systematic methodology.

Results

From 2015 to 2019, in most of the countries with three or more data points, there was a pattern of decrease in reported syphilis, gonorrhoea and chlamydia cases, which is in contrast to the EU/EEA. The number of reported cases per 100,000 population was 0.4–26.5 for syphilis, 0–18.5 for gonorrhoea and 0–43.3 for chlamydia. Four countries reported recent data on AMR in gonorrhoea to GASP, and we identified further publications from Georgia, Russia and Ukraine.

Conclusion

We found wide heterogeneity in reported rates of STI. There is a strong need to improve availability and quality of STI surveillance data in the non-EU/EEA countries.

Bacterial keratitis: identifying the areas of clinical uncertainty

Bacterial keratitis is a common corneal infection that is treated with topical antimicrobials. By the time of presentation there may already be severe visual loss from corneal ulceration and opacity, which may persist despite treatment. There are significant differences in the associated risk factors and the bacterial isolates between high income and low- or middle-income countries, so that general management guidelines may not be appropriate. Although the diagnosis of bacterial keratitis may seem intuitive there are multiple uncertainties about the criteria that are used, which impacts the interpretation of investigations and recruitment to clinical studies. Importantly, the concept that bacterial keratitis can only be confirmed by culture ignores the approximately 50% of cases clinically consistent with bacterial keratitis in which investigations are negative. The aetiology of these culture-negative cases is unknown. Currently, the estimation of bacterial susceptibility to antimicrobials is based on data from systemic administration and achievable serum or tissue concentrations, rather than relevant corneal concentrations and biological activity in the cornea. The provision to the clinician of minimum inhibitory concentrations of the antimicrobials for the isolated bacteria would be an important step forward. An increase in the prevalence of antimicrobial resistance is a concern, but the effect this has on disease outcomes is yet unclear. Virulence factors are not routinely assessed although they may affect the pathogenicity of bacteria within species and affect outcomes. New technologies have been developed to detect and kill bacteria, and their application to bacterial keratitis is discussed. In this review we present the multiple areas of clinical uncertainty that hamper research and the clinical management of bacterial keratitis, and we address some of the assumptions and dogma that have become established in the literature.

Antimicrobial resistance prediction in Neisseria gonorrhoeae: Current status and future prospects

INTRODUCTION

Several nucleic acid amplification tests (NAATs), mostly real-time PCRs, to detect antimicrobial resistance (AMR) determinants and predict AMR in Neisseria gonorrhoeae are promising, and some may be ready to apply at the point-of-care (POC), but important limitations remain with most NAATs. Next-generation sequencing (NGS) can overcome many of these limitations.Areas covered: Recent advances, with main focus on publications since 2017, in the development and use of NAATs and NGS to predict gonococcal AMR for surveillance and clinical use, and pros and cons of these tests as well as future perspectives for appropriate use of molecular AMR prediction for N. gonorrhoeae.Expert Commentary: NAATs and/or NGS for AMR prediction should supplement culture-based AMR surveillance, which will remain because it detects also AMR due to unknown AMR determinants, and translation into POC tests is imperative for the end-goal of individualized treatment, sparing ceftriaxone±azithromycin. Several challenges for direct testing of clinical, especially pharyngeal, specimens and for accurate prediction of cephalosporins and azithromycin resistance, especially using NAATs, remain. The choice of AMR prediction assay needs to carefully consider the intended use of the assay; limitations intrinsic to the AMR prediction technology, algorithms and specific to chosen methodology; specimen types analyzed; and cost-effectiveness.

Identification of isolated or mixed strains from long reads: a challenge met on Streptococcus thermophilus using a MinION sequencer

This study aimed to provide efficient recognition of bacterial strains on personal computers from MinION (Nanopore) long read data. Thanks to the fall in sequencing costs, the identification of bacteria can now proceed by whole genome sequencing. MinION is a fast, but highly error-prone sequencing device and it is a challenge to successfully identify the strain content of unknown simple or complex microbial samples. It is heavily constrained by memory management and fast access to the read and genome fragments. Our strategy involves three steps: indexing of known genomic sequences for a given or several bacterial species; a request process to assign a read to a strain by matching it to the closest reference genomes; and a final step looking for a minimum set of strains that best explains the observed reads. We have applied our method, called ORI, on 77 strains of Streptococcus thermophilus. We worked on several genomic distances and obtained a detailed classification of the strains, together with a criterion that allows merging of what we termed 'sibling' strains, only separated by a few mutations. Overall, isolated strains can be safely recognized from MinION data. For mixtures of several non-sibling strains, results depend on strain abundance.

Evaluating the potential for respiratory metagenomics to improve treatment of secondary infection and detection of nosocomial transmission on expanded COVID-19 intensive care units

BACKGROUND

Clinical metagenomics (CMg) has the potential to be translated from a research tool into routine service to improve antimicrobial treatment and infection control decisions. The SARS-CoV-2 pandemic provides added impetus to realise these benefits, given the increased risk of secondary infection and nosocomial transmission of multi-drug-resistant (MDR) pathogens linked with the expansion of critical care capacity.

METHODS

CMg using nanopore sequencing was evaluated in a proof-of-concept study on 43 respiratory samples from 34 intubated patients across seven intensive care units (ICUs) over a 9-week period during the first COVID-19 pandemic wave.

RESULTS

An 8-h CMg workflow was 92% sensitive (95% CI, 75-99%) and 82% specific (95% CI, 57-96%) for bacterial identification based on culture-positive and culture-negative samples, respectively. CMg sequencing reported the presence or absence of β-lactam-resistant genes carried by Enterobacterales that would modify the initial guideline-recommended antibiotics in every case. CMg was also 100% concordant with quantitative PCR for detecting Aspergillus fumigatus from 4 positive and 39 negative samples. Molecular typing using 24-h sequencing data identified an MDR-K. pneumoniae ST307 outbreak involving 4 patients and an MDR-C. striatum outbreak involving 14 patients across three ICUs.

CONCLUSION

CMg testing provides accurate pathogen detection and antibiotic resistance prediction in a same-day laboratory workflow, with assembled genomes available the next day for genomic surveillance. The provision of this technology in a service setting could fundamentally change the multi-disciplinary team approach to managing ICU infections. The potential to improve the initial targeted treatment and rapidly detect unsuspected outbreaks of MDR-pathogens justifies further expedited clinical assessment of CMg.

A genomic data resource for predicting antimicrobial resistance from laboratory-derived antimicrobial susceptibility phenotypes

Antimicrobial resistance (AMR) is a major global health threat that affects millions of people each year. Funding agencies worldwide and the global research community have expended considerable capital and effort tracking the evolution and spread of AMR by isolating and sequencing bacterial strains and performing antimicrobial susceptibility testing (AST). For the last several years, we have been capturing these efforts by curating data from the literature and data resources and building a set of assembled bacterial genome sequences that are paired with laboratory-derived AST data. This collection currently contains AST data for over 67 000 genomes encompassing approximately 40 genera and over 100 species. In this paper, we describe the characteristics of this collection, highlighting areas where sampling is comparatively deep or shallow, and showing areas where attention is needed from the research community to improve sampling and tracking efforts. In addition to using the data to track the evolution and spread of AMR, it also serves as a useful starting point for building machine learning models for predicting AMR phenotypes. We demonstrate this by describing two machine learning models that are built from the entire dataset to show where the predictive power is comparatively high or low. This AMR metadata collection is freely available and maintained on the Bacterial and Viral Bioinformatics Center (BV-BRC) FTP site ftp://ftp.bvbrc.org/RELEASE_NOTES/PATRIC_genomes_AMR.txt.

Analysing the fitness cost of antibiotic resistance to identify targets for combination antimicrobials

Mutations in the rifampicin (Rif)-binding site of RNA polymerase (RNAP) confer antibiotic resistance and often have global effects on transcription that compromise fitness and stress tolerance of resistant mutants. We suggested that the non-essential genome, through its impact on the bacterial transcription cycle, may represent an untapped source of targets for combination antimicrobial therapies. Using transposon sequencing, we carried out a genome-wide analysis of fitness cost in a clinically common rpoB H526Y mutant. We find that genes whose products enable increased transcription elongation rates compound the fitness costs of resistance whereas genes whose products function in cell wall synthesis and division mitigate it. We validate our findings by showing that the cell wall synthesis and division defects of rpoB H526Y result from an increased transcription elongation rate that is further exacerbated by the activity of the uracil salvage pathway and unresponsiveness of the mutant RNAP to the alarmone ppGpp. We applied our findings to identify drugs that inhibit more readily rpoB H526Y and other Rif^R alleles from the same phenotypic class. Thus, genome-wide analysis of fitness cost of antibiotic-resistant mutants should expedite the discovery of new combination therapies and delineate cellular pathways that underlie the molecular mechanisms of cost.

Heuristic and Hierarchical-Based Population Mining of Salmonella enterica Lineage I Pan-Genomes as a Platform to Enhance Food Safety

The recent incorporation of bacterial whole-genome sequencing (WGS) into Public Health laboratories has enhanced foodborne outbreak detection and source attribution. As a result, large volumes of publicly available datasets can be used to study the biology of foodborne pathogen populations at an unprecedented scale. To demonstrate the application of a heuristic and agnostic hierarchical population structure guided pan-genome enrichment analysis (PANGEA), we used populations of S. enterica lineage I to achieve two main objectives: (i) show how hierarchical population inquiry at different scales of resolution can enhance ecological and epidemiological inquiries; and (ii) identify population-specific inferable traits that could provide selective advantages in food production environments. Publicly available WGS data were obtained from NCBI database for three serovars of Salmonella enterica subsp. enterica lineage I (S. Typhimurium, S. Newport, and S. Infantis). Using the hierarchical genotypic classifications (Serovar, BAPS1, ST, cgMLST), datasets from each of the three serovars showed varying degrees of clonal structuring. When the accessory genome (PANGEA) was mapped onto these hierarchical structures, accessory loci could be linked with specific genotypes. A large heavy-metal resistance mobile element was found in the Monophasic ST34 lineage of S. Typhimurium, and laboratory testing showed that Monophasic isolates have on average a higher degree of copper resistance than the Biphasic ones. In S. Newport, an extra sugE gene copy was found among most isolates of the ST45 lineage, and laboratory testing of multiple isolates confirmed that isolates of S. Newport ST45 were on average less sensitive to the disinfectant cetylpyridimium chloride than non-ST45 isolates. Lastly, data-mining of the accessory genomic content of S. Infantis revealed two cryptic Ecotypes with distinct accessory genomic content and distinct ecological patterns. Poultry appears to be the major reservoir for Ecotype 1, and temporal analysis further suggested a recent ecological succession, with Ecotype 2 apparently being displaced by Ecotype 1. Altogether, the use of a heuristic hierarchical-based population structure analysis that includes bacterial pan-genomes (core and accessory genomes) can (1) improve genomic resolution for mapping populations and accessing epidemiological patterns; and (2) define lineage-specific informative loci that may be associated with survival in the food chain.

Recent advances in understanding and combatting Neisseria gonorrhoeae: a genomic perspective

The sexually transmitted infection (STI) gonorrhoea remains a major global public health concern. The World Health Organization (WHO) estimates that 87 million new cases in individuals who were 15 to 49 years of age occurred in 2016. The growing number of gonorrhoea cases is concerning given the rise in gonococci developing antimicrobial resistance (AMR). Therefore, a global action plan is needed to facilitate surveillance. Indeed, the WHO has made surveillance leading to the elimination of STIs (including gonorrhoea) a global health priority. The availability of whole genome sequence data offers new opportunities to combat gonorrhoea. This can be through (i) enhanced surveillance of the global prevalence of AMR, (ii) improved understanding of the population biology of the gonococcus, and (iii) opportunities to mine sequence data in the search for vaccine candidates. Here, we review the current status in Neisseria gonorrhoeae genomics. In particular, we explore how genomics continues to advance our understanding of this complex pathogen.

Perspectives and Benefits of High-Throughput Long-Read Sequencing in Microbial Ecology

Short-read, high-throughput sequencing (HTS) methods have yielded numerous important insights into microbial ecology and function. Yet, in many instances short-read HTS techniques are suboptimal, for example, by providing insufficient phylogenetic resolution or low integrity of assembled genomes. Single-molecule and synthetic long-read (SLR) HTS methods have successfully ameliorated these limitations. In addition, nanopore sequencing has generated a number of unique analysis opportunities, such as rapid molecular diagnostics and direct RNA sequencing, and both Pacific Biosciences (PacBio) and nanopore sequencing support detection of epigenetic modifications. Although initially suffering from relatively low sequence quality, recent advances have greatly improved the accuracy of long-read sequencing technologies. In spite of great technological progress in recent years, the long-read HTS methods (PacBio and nanopore sequencing) are still relatively costly, require large amounts of high-quality starting material, and commonly need specific solutions in various analysis steps. Despite these challenges, long-read sequencing technologies offer high-quality, cutting-edge alternatives for testing hypotheses about microbiome structure and functioning as well as assembly of eukaryote genomes from complex environmental DNA samples.

An omics-based framework for assessing the health risk of antimicrobial resistance genes

Antibiotic resistance genes (ARGs) are widespread among bacteria. However, not all ARGs pose serious threats to public health, highlighting the importance of identifying those that are high-risk. Here, we developed an ‘omics-based’ framework to evaluate ARG risk considering human-associated-enrichment, gene mobility, and host pathogenicity. Our framework classifies human-associated, mobile ARGs (3.6% of all ARGs) as the highest risk, which we further differentiate as ‘current threats’ (Rank I; 3%) - already present among pathogens - and ‘future threats’ (Rank II; 0.6%) - novel resistance emerging from non-pathogens. Our framework identified 73 ‘current threat’ ARG families. Of these, 35 were among the 37 high-risk ARGs proposed by the World Health Organization and other literature; the remaining 38 were significantly enriched in hospital plasmids. By evaluating all pathogen genomes released since framework construction, we confirmed that ARGs that recently transferred into pathogens were significantly enriched in Rank II (‘future threats’). Lastly, we applied the framework to gut microbiome genomes from fecal microbiota transplantation donors. We found that although ARGs were widespread (73% of genomes), only 8.9% of genomes contained high-risk ARGs. Our framework provides an easy-to-implement approach to identify current and future antimicrobial resistance threats, with potential clinical applications including reducing risk of microbiome-based interventions.

Antibiotic resistance genes are common but not all are of high risk to human health. Here, the authors develop an omics-based framework for ranking genes by risk that incorporates level of enrichment in human associated environments, gene mobility, and host pathogenicity.

Predictive Antibiotic Susceptibility Testing by Next-Generation Sequencing for Periprosthetic Joint Infections: Potential and Limitations

Joint replacement surgeries are one of the most frequent medical interventions globally. Infections of prosthetic joints are a major health challenge and typically require prolonged or even indefinite antibiotic treatment. As multidrug-resistant pathogens continue to rise globally, novel diagnostics are critical to ensure appropriate treatment and help with prosthetic joint infections (PJI) management. To this end, recent studies have shown the potential of molecular methods such as next-generation sequencing to complement established phenotypic, culture-based methods. Together with advanced bioinformatics approaches, next-generation sequencing can provide comprehensive information on pathogen identity as well as antimicrobial susceptibility, potentially enabling rapid diagnosis and targeted therapy of PJIs. In this review, we summarize current developments in next generation sequencing based predictive antibiotic susceptibility testing and discuss potential and limitations for common PJI pathogens.

Isothermal amplifications - a comprehensive review on current methods

The introduction of nucleic acid amplification techniques has revolutionized the field of medical diagnostics in the last decade. The advent of PCR catalyzed the increasing application of DNA, not just for molecular cloning but also for molecular based diagnostics. Since the introduction of PCR, a deeper understanding of molecular mechanisms and enzymes involved in DNA/RNA replication has spurred the development of novel methods devoid of temperature cycling. Isothermal amplification methods have since been introduced utilizing different mechanisms, enzymes, and conditions. The ease with which isothermal amplification methods have allowed nucleic acid amplification to be carried out has had a profound impact on the way molecular diagnostics are being designed after the turn of the millennium. With all the advantages isothermal amplification brings, the issues or complications surrounding each method are heterogeneous making it difficult to identify the best approach for an end-user. This review pays special attention to the various isothermal amplification methods by classifying them based on the mechanistic characteristics which include reaction formats, amplification information, promoter, strand break, and refolding mechanisms. We would also compare the efficiencies and usefulness of each method while highlighting the potential applications and detection methods involved. This review will serve as an overall outlook on the journey and development of isothermal amplification methods as a whole.

Developing a model to predict individualised treatment for gonorrhoea: a modelling study

OBJECTIVE

To develop a tool predicting individualised treatment for gonorrhoea, enabling treatment with previously recommended antibiotics, to reduce use of last-line treatment ceftriaxone.

DESIGN

A modelling study.

SETTING

England and Wales.

PARTICIPANTS

Individuals accessing sentinel health services.

INTERVENTION

Developing an Excel model which uses participants' demographic, behavioural and clinical characteristics to predict susceptibility to legacy antibiotics. Model parameters were calculated using data for 2015-2017 from the Gonococcal Resistance to Antimicrobials Surveillance Programme.

MAIN OUTCOME MEASURES

Estimated number of doses of ceftriaxone saved, and number of people delayed effective treatment, by model use in clinical practice. Model outputs are the predicted risk of resistance to ciprofloxacin, azithromycin, penicillin and cefixime, in groups of individuals with different combinations of characteristics (gender, sexual orientation, number of recent sexual partners, age, ethnicity), and a treatment recommendation.

RESULTS

Between 2015 and 2017, 8013 isolates were collected: 64% from men who have sex with men, 18% from heterosexual men and 18% from women. Across participant subgroups, stratified by all predictors, resistance prevalence was high for ciprofloxacin (range: 11%-51%) and penicillin (range: 6%-33%). Resistance prevalence for azithromycin and cefixime ranged from 0% to 13% and for ceftriaxone it was 0%. Simulating model use, 88% of individuals could be given cefixime and 10% azithromycin, saving 97% of ceftriaxone doses, with 1% of individuals delayed effective treatment.

CONCLUSIONS

Using demographic and behavioural characteristics, we could not reliably identify a participant subset in which ciprofloxacin or penicillin would be effective. Cefixime resistance was almost universally low; however, substituting ceftriaxone for near-uniform treatment with cefixime risks re-emergence of resistance to cefixime and ceftriaxone. Several subgroups had low azithromycin resistance, but widespread azithromycin monotherapy risks resistance at population level. However, this dataset had limitations; further exploration of individual characteristics to predict resistance to a wider range of legacy antibiotics may still be appropriate.

Predicting Antimicrobial Resistance Using Partial Genome Alignments

Antimicrobial resistance (AMR) is an important global health threat that impacts millions of people worldwide each year. Developing methods that can detect and predict AMR phenotypes can help to mitigate the spread of AMR by informing clinical decision making and appropriate mitigation strategies. Many bioinformatic methods have been developed for predicting AMR phenotypes from whole-genome sequences and AMR genes, but recent studies have indicated that predictions can be made from incomplete genome sequence data. In order to more systematically understand this, we built random forest-based machine learning classifiers for predicting susceptible and resistant phenotypes for Klebsiella pneumoniae (1,640 strains), Mycobacterium tuberculosis (2,497 strains), and Salmonella enterica (1,981 strains). We started by building models from alignments that were based on a reference chromosome for each species. We then subsampled each chromosomal alignment and built models for the resulting subalignments, finding that very small regions, representing approximately 0.1 to 0.2% of the chromosome, are predictive. In K. pneumoniae, M. tuberculosis, and S. enterica, the subalignments are able to predict multiple AMR phenotypes with at least 70% accuracy, even though most do not encode an AMR-related function. We used these models to identify regions of the chromosome with high and low predictive signals. Finally, subalignments that retain high accuracy across larger phylogenetic distances were examined in greater detail, revealing genes and intergenic regions with potential links to AMR, virulence, transport, and survival under stress conditions.

IMPORTANCE Antimicrobial resistance causes thousands of deaths annually worldwide. Understanding the regions of the genome that are involved in antimicrobial resistance is important for developing mitigation strategies and preventing transmission. Machine learning models are capable of predicting antimicrobial resistance phenotypes from bacterial genome sequence data by identifying resistance genes, mutations, and other correlated features. They are also capable of implicating regions of the genome that have not been previously characterized as being involved in resistance. In this study, we generated global chromosomal alignments for Klebsiella pneumoniae, Mycobacterium tuberculosis, and Salmonella enterica and systematically searched them for small conserved regions of the genome that enable the prediction of antimicrobial resistance phenotypes. In addition to known antimicrobial resistance genes, this analysis identified genes involved in virulence and transport functions, as well as many genes with no previous implication in antimicrobial resistance.

Genomic Prediction of Antimicrobial Resistance: Ready or Not, Here It Comes!

BACKGROUND

Next-generation sequencing (NGS) technologies are being used to predict antimicrobial resistance. The field is evolving rapidly and transitioning out of the research setting into clinical use. Clinical laboratories are evaluating the accuracy and utility of genomic resistance prediction, including methods for NGS, downstream bioinformatic pipeline components, and the clinical settings in which this type of testing should be offered.

CONTENT

We describe genomic sequencing as it pertains to predicting antimicrobial resistance in clinical isolates and samples. We elaborate on current methodologies and workflows to perform this testing and summarize the current state of genomic resistance prediction in clinical settings. To highlight this aspect, we include 3 medically relevant microorganism exemplars: Mycobacterium tuberculosis, Staphylococcus aureus, and Neisseria gonorrhoeae. Last, we discuss the future of genomic-based resistance detection in clinical microbiology laboratories.

SUMMARY

Antimicrobial resistance prediction by genomic approaches is in its infancy for routine patient care. Genomic approaches have already added value to the current diagnostic testing landscape in specific circumstances and will play an increasingly important role in diagnostic microbiology. Future advancements will shorten turnaround time, reduce costs, and improve our analysis and interpretation of clinically actionable results.

Antibiotic resistance: Time of synthesis in a post-genomic age

Antibiotic resistance has been highlighted by international organizations, including World Health Organization, World Bank and United Nations, as one of the most relevant global health problems. Classical approaches to study this problem have focused in infected humans, mainly at hospitals. Nevertheless, antibiotic resistance can expand through different ecosystems and geographical allocations, hence constituting a One-Health, Global-Health problem, requiring specific integrative analytic tools. Antibiotic resistance evolution and transmission are multilayer, hierarchically organized processes with several elements (from genes to the whole microbiome) involved. However, their study has been traditionally gene-centric, each element independently studied. The development of robust-economically affordable whole genome sequencing approaches, as well as other -omic techniques as transcriptomics and proteomics, is changing this panorama. These technologies allow the description of a system, either a cell or a microbiome as a whole, overcoming the problems associated with gene-centric approaches. We are currently at the time of combining the information derived from -omic studies to have a more holistic view of the evolution and spread of antibiotic resistance. This synthesis process requires the accurate integration of -omic information into computational models that serve to analyse the causes and the consequences of acquiring AR, fed by curated databases capable of identifying the elements involved in the acquisition of resistance. In this review, we analyse the capacities and drawbacks of the tools that are currently in use for the global analysis of AR, aiming to identify the more useful targets for effective corrective interventions.

Rapid metagenomics for diagnosis of bloodstream and respiratory tract nosocomial infections: current status and future prospects

Introduction: Nosocomial infections represent a major problem for the health-care systems worldwide. Currently, diagnosis relies on microbiological culture, which is slow and has poor sensitivity. While waiting for a diagnosis, patients are treated with empiric broad spectrum antimicrobials, which are often inappropriate for the infecting pathogen. This results in poor patient outcomes, poor antimicrobial stewardship and increased costs for health-care systems.Areas covered: Clinical metagenomics (CMg), the application of metagenomic sequencing for the diagnosis of infection, has the potential to become a viable alternative to culture that can offer rapid results with high accuracy. In this article, we review current CMg methods for the diagnosis of nosocomial bloodstream (BSI) and lower respiratory-tract infections (LRTI).Expert opinion: CMg approaches are more accurate in LRTI compared to BSI. This is because BSIs are caused by low pathogen numbers in a high background of human cells. To overcome this, most approaches focus on cell-free DNA, but, to date, these tests are not accurate enough yet to replace blood culture. The higher pathogen numbers in LRTI samples make this a more suitable for CMg and accurate approaches have been developed, which are likely to be implemented in hospitals within the next 2-5 years.

Simplitigs as an efficient and scalable representation of de Bruijn graphs

de Bruijn graphs play an essential role in bioinformatics, yet they lack a universal scalable representation. Here, we introduce simplitigs as a compact, efficient, and scalable representation, and ProphAsm, a fast algorithm for their computation. For the example of assemblies of model organisms and two bacterial pan-genomes, we compare simplitigs to unitigs, the best existing representation, and demonstrate that simplitigs provide a substantial improvement in the cumulative sequence length and their number. When combined with the commonly used Burrows-Wheeler Transform index, simplitigs reduce memory, and index loading and query times, as demonstrated with large-scale examples of GenBank bacterial pan-genomes.