Prediction of pyrazinamide resistance in Mycobacterium tuberculosis using structure-based machine-learning approaches.

Background: Pyrazinamide is one of four first-line antibiotics used to treat tuberculosis; however, antibiotic susceptibility testing for pyrazinamide is challenging. Resistance to pyrazinamide is primarily driven by genetic variation in pncA, encoding an enzyme that converts pyrazinamide into its active form. Methods: We curated a dataset of 664 non-redundant, missense amino acid mutations in PncA with associated high-confidence phenotypes from published studies and then trained three different machine-learning models to predict pyrazinamide resistance. All models had access to a range of protein structural-, chemical- and sequence-based features. Results: The best model, a gradient-boosted decision tree, achieved a sensitivity of 80.2% and a specificity of 76.9% on the hold-out test dataset. The clinical performance of the models was then estimated by predicting the binary pyrazinamide resistance phenotype of 4027 samples harbouring 367 unique missense mutations in pncA derived from 24 231 clinical isolates. Conclusions: This work demonstrates how machine learning can enhance the sensitivity/specificity of pyrazinamide resistance prediction in genetics-based clinical microbiology workflows, highlights novel mutations for future biochemical investigation, and is a proof of concept for using this approach in other drugs.

Discordance between different bioinformatic methods for identifying resistance genes from short-read genomic data, with a focus on Escherichia coli.

Several bioinformatics genotyping algorithms are now commonly used to characterize antimicrobial resistance (AMR) gene profiles in whole-genome sequencing (WGS) data, with a view to understanding AMR epidemiology and developing resistance prediction workflows using WGS in clinical settings. Accurately evaluating AMR in Enterobacterales, particularly Escherichia coli, is of major importance, because this is a common pathogen. However, robust comparisons of different genotyping approaches on relevant simulated and large real-life WGS datasets are lacking. Here, we used both simulated datasets and a large set of real E. coli WGS data (n=1818 isolates) to systematically investigate genotyping methods in greater detail. Simulated constructs and real sequences were processed using four different bioinformatic programs (ABRicate, ARIBA, KmerResistance and SRST2, run with the ResFinder database) and their outputs compared. For simulation tests where 3079 AMR gene variants were inserted into random sequence constructs, KmerResistance was correct for 3076 (99.9 %) simulations, ABRicate for 3054 (99.2 %), ARIBA for 2783 (90.4 %) and SRST2 for 2108 (68.5 %). For simulation tests where two closely related gene variants were inserted into random sequence constructs, KmerResistance identified the correct alleles in 35 338/46 318 (76.3 %) simulations, ABRicate identified them in 11 842/46 318 (25.6 %) simulations, ARIBA identified them in 1679/46 318 (3.6 %) simulations and SRST2 identified them in 2000/46 318 (4.3 %) simulations. In real data, across all methods, 1392/1818 (76 %) isolates had discrepant allele calls for at least 1 gene. In addition to highlighting areas for improvement in challenging scenarios, (e.g. identification of AMR genes at <10× coverage, identifying multiple closely related AMR genes present in the same sample), our evaluations identified some more systematic errors that could be readily soluble, such as repeated misclassification (i.e. naming) of genes as shorter variants of the same gene present within the reference resistance gene database. Such naming errors accounted for at least 2530/4321 (59 %) of the discrepancies seen in real data. Moreover, many of the remaining discrepancies were likely 'artefactual', with reporting of cut-off differences accounting for at least 1430/4321 (33 %) discrepants. Whilst we found that comparing outputs generated by running multiple algorithms on the same dataset could identify and resolve these algorithmic artefacts, the results of our evaluations emphasize the need for developing new and more robust genotyping algorithms to further improve accuracy and performance.

A crowd of BashTheBug volunteers reproducibly and accurately measure the minimum inhibitory concentrations of 13 antitubercular drugs from photographs of 96-well broth microdilution plates.

Tuberculosis is a respiratory disease that is treatable with antibiotics. An increasing prevalence of resistance means that to ensure a good treatment outcome it is desirable to test the susceptibility of each infection to different antibiotics. Conventionally, this is done by culturing a clinical sample and then exposing aliquots to a panel of antibiotics, each being present at a pre-determined concentration, thereby determining if the sample isresistant or susceptible to each sample. The minimum inhibitory concentration (MIC) of a drug is the lowestconcentration that inhibits growth and is a more useful quantity but requires each sample to be tested at a range ofconcentrations for each drug. Using 96-well broth micro dilution plates with each well containing a lyophilised pre-determined amount of an antibiotic is a convenient and cost-effective way to measure the MICs of several drugs at once for a clinical sample. Although accurate, this is still an expensive and slow process that requires highly-skilled and experienced laboratory scientists. Here we show that, through the BashTheBug project hosted on the Zooniverse citizen science platform, a crowd of volunteers can reproducibly and accurately determine the MICs for 13 drugs and that simply taking the median or mode of 11-17 independent classifications is sufficient. There is therefore a potential role for crowds to support (but not supplant) the role of experts in antibiotic susceptibility testing.

Tuberculosis is a bacterial respiratory infection that kills about 1.4 million people worldwide each year. While antibiotics can cure the condition, the bacterium responsible for this disease, Mycobacterium tuberculosis, is developing resistance to these treatments. Choosing which antibiotics to use to treat the infection more carefully may help to combat the growing threat of drug-resistant bacteria. One way to find the best choice is to test how an antibiotic affects the growth of M. tuberculosis in the laboratory. To speed up this process, laboratories test multiple drugs simultaneously. They do this by growing bacteria on plates with 96 wells and injecting individual antibiotics in to each well at different concentrations. The Comprehensive Resistance Prediction for Tuberculosis (CRyPTIC) consortium has used this approach to collect and analyse bacteria from over 20,000 tuberculosis patients. An image of the 96-well plate is then captured and the level of bacterial growth in each well is assessed by laboratory scientists. But this work is difficult, time-consuming, and subjective, even for tuberculosis experts. Here, Fowler et al. show that enlisting citizen scientists may help speed up this process and reduce errors that arise from analysing such a large dataset. In April 2017, Fowler et al. launched the project 'BashTheBug' on the Zooniverse citizen science platform where anyone can access and analyse the images from the CRyPTIC consortium. They found that a crowd of inexperienced volunteers were able to consistently and accurately measure the concentration of antibiotics necessary to inhibit the growth of M. tuberculosis. If the concentration is above a pre-defined threshold, the bacteria are considered to be resistant to the treatment. A consensus result could be reached by calculating the median value of the classifications provided by as few as 17 different BashTheBug participants. The work of BashTheBug volunteers has reduced errors in the CRyPTIC project data, which has been used for several other studies. For instance, the World Health Organization (WHO) has also used the data to create a catalogue of genetic mutations associated with antibiotics resistance in M. tuberculosis. Enlisting citizen scientists has accelerated research on tuberculosis and may help with other pressing public health concerns.

Mycobacterium tuberculosis transmission in Birmingham, UK, 2009-19: An observational study.

Background: Over 10-years of whole-genome sequencing (WGS) of Mycobacterium tuberculosis in Birmingham presents an opportunity to explore epidemiological trends and risk factors for transmission in new detail. Methods: Between 1st January 2009 and 15th June 2019, we obtained the first WGS isolate from every patient resident in a postcode district covered by Birmingham's centralised tuberculosis service. Data on patients' sex, country of birth, social risk-factors, anatomical locus of disease, and strain lineage were collected. Poisson harmonic regression was used to assess seasonal variation in case load and a mixed-effects multivariable Cox proportionate hazards model was used to assess risk factors for a future case arising in clusters defined by a 5 single nucleotide polymorphism (SNP) threshold, and by 12 SNPs in a sensitivity analysis. Findings: 511/1653 (31%) patients were genomically clustered with another. A seasonal variation in diagnoses was observed, peaking in spring, but only among clustered cases. Risk-factors for a future clustered case included UK-birth (aHR=2·03 (95%CI 1·35-3·04), p < 0·001), infectious (pulmonary/laryngeal/miliary) tuberculosis (aHR=3·08 (95%CI 1·98-4·78), p < 0·001), and M. tuberculosis lineage 3 (aHR=1·91 (95%CI 1·03-3·56), p = 0·041) and 4 (aHR=2·27 (95%CI 1·21-4·26), p = 0·011), vs. lineage 1. Similar results pertained to 12 SNP clusters, for which social risk-factors were also significant (aHR 1·72 (95%CI 1·02-2·93), p = 0·044). There was marked heterogeneity in transmission patterns between postcode districts. Interpretation: There is seasonal variation in the diagnosis of genomically clustered, but not non-clustered, cases. Risk factors for clustering include UK-birth, infectious forms of tuberculosis, and infection with lineage 3 or 4. Funding: Wellcome Trust, MRC, UKHSA.

Flanker: a tool for comparative genomics of gene flanking regions.

Analysing the flanking sequences surrounding genes of interest is often highly relevant to understanding the role of mobile genetic elements (MGEs) in horizontal gene transfer, particular for antimicrobial-resistance genes. Here, we present Flanker, a Python package that performs alignment-free clustering of gene flanking sequences in a consistent format, allowing investigation of MGEs without prior knowledge of their structure. These clusters, known as 'flank patterns' (FPs), are based on Mash distances, allowing for easy comparison of similarity across sequences. Additionally, Flanker can be flexibly parameterized to fine-tune outputs by characterizing upstream and downstream regions separately, and investigating variable lengths of flanking sequence. We apply Flanker to two recent datasets describing plasmid-associated carriage of important carbapenemase genes (blaOXA-48 and blaKPC-2/3) and show that it successfully identifies distinct clusters of FPs, including both known and previously uncharacterized structural variants. For example, Flanker identified four Tn4401 profiles that could not be sufficiently characterized using TETyper or MobileElementFinder, demonstrating the utility of Flanker for flanking-gene characterization. Similarly, using a large (n=226) European isolate dataset, we confirm findings from a previous smaller study demonstrating association between Tn1999.2 and blaOXA-48 upregulation and demonstrate 17 FPs (compared to the 5 previously identified). More generally, the demonstration in this study that FPs are associated with geographical regions and antibiotic-susceptibility phenotypes suggests that they may be useful as epidemiological markers. Flanker is freely available under an MIT license at https://github.com/wtmatlock/flanker.

Ten-year longitudinal molecular epidemiology study of Escherichia coli and Klebsiella species bloodstream infections in Oxfordshire, UK.

BACKGROUND: The incidence of Gram-negative bloodstream infections (BSIs), predominantly caused by Escherichia coli and Klebsiella species, continues to increase; however, the causes of this are unclear and effective interventions are therefore hard to design. METHODS: In this study, we sequenced 3468 unselected isolates over a decade in Oxfordshire (UK) and linked this data to routinely collected electronic healthcare records and mandatory surveillance reports. We annotated genomes for clinically relevant genes, contrasting the distribution of these within and between species, and compared incidence trends over time using stacked negative binomial regression. RESULTS: We demonstrate that the observed increases in E. coli incidence were not driven by the success of one or more sequence types (STs); instead, four STs continue to dominate a stable population structure, with no evidence of adaptation to hospital/community settings. Conversely in Klebsiella spp., most infections are caused by sporadic STs with the exception of a local drug-resistant outbreak strain (ST490). Virulence elements are highly structured by ST in E. coli but not Klebsiella spp. where they occur in a diverse spectrum of STs and equally across healthcare and community settings. Most clinically hypervirulent (i.e. community-onset) Klebsiella BSIs have no known acquired virulence loci. Finally, we demonstrate a diverse but largely genus-restricted mobilome with close associations between antimicrobial resistance (AMR) genes and insertion sequences but not typically specific plasmid replicon types, consistent with the dissemination of AMR genes being highly contingent on smaller mobile genetic elements (MGEs). CONCLUSIONS: Our large genomic study highlights distinct differences in the molecular epidemiology of E. coli and Klebsiella BSIs and suggests that no single specific pathogen genetic factors (e.g. AMR/virulence genes/sequence type) are likely contributing to the increasing incidence of BSI overall, that association with AMR genes in E. coli is a contributor to the increasing number of E. coli BSIs, and that more attention should be given to AMR gene associations with non-plasmid MGEs to try and understand horizontal gene transfer networks.

An end-to-end heterogeneous graph attention network for Mycobacterium tuberculosis drug-resistance prediction.

Antimicrobial resistance (AMR) poses a threat to global public health. To mitigate the impacts of AMR, it is important to identify the molecular mechanisms of AMR and thereby determine optimal therapy as early as possible. Conventional machine learning-based drug-resistance analyses assume genetic variations to be homogeneous, thus not distinguishing between coding and intergenic sequences. In this study, we represent genetic data from Mycobacterium tuberculosis as a graph, and then adopt a deep graph learning method-heterogeneous graph attention network ('HGAT-AMR')-to predict anti-tuberculosis (TB) drug resistance. The HGAT-AMR model is able to accommodate incomplete phenotypic profiles, as well as provide 'attention scores' of genes and single nucleotide polymorphisms (SNPs) both at a population level and for individual samples. These scores encode the inputs, which the model is 'paying attention to' in making its drug resistance predictions. The results show that the proposed model generated the best area under the receiver operating characteristic (AUROC) for isoniazid and rifampicin (98.53 and 99.10%), the best sensitivity for three first-line drugs (94.91% for isoniazid, 96.60% for ethambutol and 90.63% for pyrazinamide), and maintained performance when the data were associated with incomplete phenotypes (i.e. for those isolates for which phenotypic data for some drugs were missing). We also demonstrate that the model successfully identifies genes and SNPs associated with drug resistance, mitigating the impact of resistance profile while considering particular drug resistance, which is consistent with domain knowledge.

SARS-CoV-2 within-host diversity and transmission.

Extensive global sampling and sequencing of the pandemic virus severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) have enabled researchers to monitor its spread and to identify concerning new variants. Two important determinants of variant spread are how frequently they arise within individuals and how likely they are to be transmitted. To characterize within-host diversity and transmission, we deep-sequenced 1313 clinical samples from the United Kingdom. SARS-CoV-2 infections are characterized by low levels of within-host diversity when viral loads are high and by a narrow bottleneck at transmission. Most variants are either lost or occasionally fixed at the point of transmission, with minimal persistence of shared diversity, patterns that are readily observable on the phylogenetic tree. Our results suggest that transmission-enhancing and/or immune-escape SARS-CoV-2 variants are likely to arise infrequently but could spread rapidly if successfully transmitted.

Genomic network analysis of environmental and livestock F-type plasmid populations.

F-type plasmids are diverse and of great clinical significance, often carrying genes conferring antimicrobial resistance (AMR) such as extended-spectrum ß-lactamases, particularly in Enterobacterales. Organising this plasmid diversity is challenging, and current knowledge is largely based on plasmids from clinical settings. Here, we present a network community analysis of a large survey of F-type plasmids from environmental (influent, effluent and upstream/downstream waterways surrounding wastewater treatment works) and livestock settings. We use a tractable and scalable methodology to examine the relationship between plasmid metadata and network communities. This reveals how niche (sampling compartment and host genera) partition and shape plasmid diversity. We also perform pangenome-style analyses on network communities. We show that such communities define unique combinations of core genes, with limited overlap. Building plasmid phylogenies based on alignments of these core genes, we demonstrate that plasmid accessory function is closely linked to core gene content. Taken together, our results suggest that stable F-type plasmid backbone structures can persist in environmental settings while allowing dramatic variation in accessory gene content that may be linked to niche adaptation. The association of F-type plasmids with AMR may reflect their suitability for rapid niche adaptation.

Diagnosis of SARS-CoV-2 Infection with LamPORE, a High-Throughput Platform Combining Loop-Mediated Isothermal Amplification and Nanopore Sequencing.

LamPORE is a novel diagnostic platform for the detection of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) RNA combining loop-mediated isothermal amplification with nanopore sequencing, which could potentially be used to analyze thousands of samples per day on a single instrument. We evaluated the performance of LamPORE against reverse transcriptase PCR (RT-PCR) using RNA extracted from spiked respiratory samples and stored nose and throat swabs collected at two UK hospitals. The limit of detection of LamPORE was 10 genome copies/µl of extracted RNA, which is above the limit achievable by RT-PCR, but was not associated with a significant reduction of sensitivity in clinical samples. Positive clinical specimens came mostly from patients with acute symptomatic infection, and among them, LamPORE had a diagnostic sensitivity of 99.1% (226/228; 95% confidence interval [CI], 96.9% to 99.9%). Among negative clinical specimens, including 153 with other respiratory pathogens detected, LamPORE had a diagnostic specificity of 99.6% (278/279; 98.0% to 100.0%). Overall, 1.4% (7/514; 0.5% to 2.9%) of samples produced an indeterminate result on first testing, and repeat LamPORE testing on the same RNA extract had a reproducibility of 96.8% (478/494; 94.8% to 98.1%). LamPORE has a similar performance as RT-PCR for the diagnosis of SARS-CoV-2 infection in symptomatic patients and offers a promising approach to high-throughput testing.

Stringent thresholds in SARS-CoV-2 IgG assays lead to under-detection of mild infections.

BACKGROUND: Thresholds for SARS-CoV-2 antibody assays have typically been determined using samples from symptomatic, often hospitalised, patients. In this setting the sensitivity and specificity of the best performing assays can both exceed 98%. However, antibody assay performance following mild infection is less clear. METHODS: We assessed quantitative IgG responses in a cohort of healthcare workers in Oxford, UK, with a high pre-test probability of Covid-19, in particular the 991/11,475(8.6%) who reported loss of smell/taste. We use anosmia/ageusia and other risk factors as probes for Covid-19 infection potentially undiagnosed by immunoassays by investigating their relationship with antibody readings either side of assay thresholds. RESULTS: The proportion of healthcare workers reporting anosmia/ageusia increased at antibody readings below diagnostic thresholds using an in-house ELISA (n = 9324) and the Abbott Architect chemiluminescent microparticle immunoassay (CMIA; n = 11,324): 426/906 (47%) reported anosmia/ageusia with a positive ELISA, 59/449 (13.1%) with high-negative and 326/7969 (4.1%) with low-negative readings. Similarly, by CMIA, 518/1093 (47.4%) with a positive result reported anosmia/ageusia, 106/686 (15.5%) with a high-negative and 358/9563 (3.7%) with a low-negative result. Adjusting for the proportion of staff reporting anosmia/ageusia suggests the sensitivity of both assays in mild infection is lower than previously reported: Oxford ELISA 89.8% (95%CI 86.6-92.8%) and Abbott CMIA 79.3% (75.9-82.7%). CONCLUSION: Following mild SARS-CoV-2 infection 10-30% of individuals may have negative immunoassay results. While lowered diagnostic thresholds may result in unacceptable specificity, our findings have implications for epidemiological analyses and result interpretation in individuals with a high pre-test probability. Samples from mild PCR-confirmed infections should be included in SARS-CoV-2 immunoassay evaluations.

Ten Years of Population-Level Genomic Escherichia coli and Klebsiella pneumoniae Serotype Surveillance Informs Vaccine Development for Invasive Infections.

BACKGROUND: The incidence of bloodstream infections (BSIs) caused by Escherichia coli and Klebsiella pneumoniae is increasing, with substantial associated morbidity, mortality, and antimicrobial resistance. Unbiased serotyping studies to guide vaccine target selection are limited. METHODS: We conducted unselected, population-level genomic surveillance of bloodstream E. coli and Klebsiella pneumoniae isolates from 2008 to 2018 in Oxfordshire, United Kingdom. We supplemented this with an analysis of publicly available global sequencing data (n = 3678). RESULTS: We sequenced 3478 E. coli isolates (3278 passed quality control) and 556 K. pneumoniae isolates (535 [K-antigen] and 549 [O-antigen] passed quality control). The 4 most common E. coli O-antigens (O1/O2/O6/O25) were identified in 1499/3278 isolates; the incidence of these O-types increased over time (incidence rate ratio per year [IRRy] = 1.14, 95% confidence interval [CI]: 1.11-1.16). These O-types accounted for 616/1434 multidrug-resistant (MDR) and 173/256 extended-spectrum beta-lactamase (ESBL)-resistant isolates in Oxfordshire but only 19/90 carbapenem-resistant isolates across all studies. For Klebsiella pneumoniae, the most common O-antigens (O2v2/O1v1/O3b/O1v2) accounted for 410/549 isolates; the incidence of BSIs caused by these also increased annually (IRRy = 1.09; 95% CI: 1.05-1.12). These O-types accounted for 122/148 MDR and 106/123 ESBL isolates in Oxfordshire and 557/734 carbapenem-resistant isolates across all studies. Conversely we observed substantial capsular antigen diversity. Analysis of 3678 isolates from global studies demonstrated the generalizability of these findings. For E. coli, based on serotyping, the ExPEC4V and ExPEC10V vaccines under investigation would cover 46% and 72% of Oxfordshire isolates respectively, and 47% and 71% of MDR isolates. CONCLUSIONS: O-antigen targeted vaccines may be useful in reducing the morbidity, mortality, and antimicrobial resistance associated with E. coli and K. pneumoniae BSIs.

Optimized use of Oxford Nanopore flowcells for hybrid assemblies.

Hybrid assemblies are highly valuable for studies of Enterobacteriaceae due to their ability to fully resolve the structure of mobile genetic elements, such as plasmids, which are involved in the carriage of clinically important genes (e.g. those involved in antimicrobial resistance/virulence). The widespread application of this technique is currently primarily limited by cost. Recent data have suggested that non-inferior, and even superior, hybrid assemblies can be produced using a fraction of the total output from a multiplexed nanopore [Oxford Nanopore Technologies (ONT)] flowcell run. In this study we sought to determine the optimal minimal running time for flowcells when acquiring reads for hybrid assembly. We then evaluated whether the ONT wash kit might allow users to exploit shorter running times by sequencing multiple libraries per flowcell. After 24 h of sequencing, most chromosomes and plasmids had circularized and there was no benefit associated with longer running times. Quality was similar at 12 h, suggesting that shorter running times are likely to be acceptable for certain applications (e.g. plasmid genomics). The ONT wash kit was highly effective in removing DNA between libraries. Contamination between libraries did not appear to affect subsequent hybrid assemblies, even when the same barcodes were used successively on a single flowcell. Utilizing shorter run times in combination with between-library nuclease washes allows at least 36 Enterobacteriaceae isolates to be sequenced per flowcell, significantly reducing the per-isolate sequencing cost. Ultimately this will facilitate large-scale studies utilizing hybrid assembly, advancing our understanding of the genomics of key human pathogens.

Risk Factors Associated with Carbapenemase-Producing Enterobacterales (CPE) Positivity in the Hospital Wastewater Environment.

Hospital wastewater is an increasingly recognized reservoir for resistant Gram-negative organisms. Factors involved in establishment and persistence of Klebsiella pneumoniae carbapenemase-producing organisms (KPCOs) in hospital wastewater plumbing are unclear. This study was conducted at a hospital with endemic KPCOs linked to wastewater reservoirs and robust patient perirectal screening for silent KPCO carriage. Over 5 months, both rooms occupied and rooms not occupied by KPCO-positive patients were sampled at three wastewater sites within each room (sink drain, sink P-trap, and toilet or hopper). Risk factors for KPCO positivity were assessed using logistic regression. Whole-genome sequencing (WGS) identified environmental seeding by KPCO-positive patients. A total of 219/475 (46%) room sampling events were KPCO positive in at least one wastewater site. KPCO-positive patient exposure was associated with increased risk of environmental positivity for the room and toilet/hopper. Previous positivity and intensive care unit room type were consistently associated with increased risk. Tube feeds were associated with increased risk for the drain, while exposure to patients with Clostridioides difficile was associated with decreased risk. Urinary catheter exposure was associated with increased risk of P-trap positivity. P-trap heaters reduced risk of P-trap and sink drain positivity. WGS identified genomically linked environmental seeding in 6 of 99 room occupations by 40 KPCO-positive patients. In conclusion, KPCO-positive patients seed the environment in at least 6% of opportunities; once positive for KPCOs, wastewater sites are at greater risk of being positive subsequently. Increased nutrient exposure, e.g., due to tube food disposal down sinks, may increase risk; frequent flushing may be protective.IMPORTANCEKlebsiella pneumoniae carbapenemase-producing organisms (KPCOs) are bacteria that are resistant to most antibiotics and thus are challenging to treat when they cause infections in patients. These organisms can be acquired by patients who are hospitalized for other reasons, complicating their hospital stay and even leading to death. Hospital wastewater sites, such as sink drains and toilets, have played a role in many reported outbreaks over the past decade. The significance of our research is in identifying risk factors for environmental positivity for KPCOs, which will facilitate further work to prevent transmission of these organisms to patients from the hospital environment.

Route and duration of antibiotic therapy in acute cellulitis: A systematic review and meta-analysis of the effectiveness and harms of antibiotic treatment.

OBJECTIVES: Compared with guideline recommendations, antibiotic overuse is common in treating cellulitis. We conducted a systematic review and meta-analyses on antibiotic route and duration of treatment for cellulitis in adults and children. METHODS: We searched MEDLINE, EMBASE and trial registries from inception to Dec 11, 2019 for interventional and observational studies of antibiotic treatment for cellulitis. Exclusions included case series/reports, pre-septal/orbital cellulitis and non-English language articles. Random-effects meta-analyses were used to produce summary relative risk (RR) estimates for our primary outcome of clinical response. PROSPERO: CRD42018100602. RESULTS: We included 47/8423 articles, incorporating data from eleven trials (1855 patients) in two meta-analyses. The overall risk of bias was moderate. Only two trials compared the same antibiotic agent in each group. We found no evidence of difference in clinical response rates for antibiotic route or duration (RR(oral:IV)=1.12, 95%CI 0.98-1.27, I2=32% and RR(shorter:longer)=0.99, 95%CI 0•96-1.03, I2 = 0%, respectively). Findings were consistent in observational studies. Follow-up data beyond 30 days were sparse. CONCLUSIONS: The evidence base for antibiotic treatment decisions in cellulitis is flawed by biased comparisons, short follow-up and lack of data around harms of antibiotic overuse. Future research should focus on developing patient-tailored antibiotic prescribing for cellulitis to reduce unnecessary antibiotic use.

DNA Thermo-Protection Facilitates Whole-Genome Sequencing of Mycobacteria Direct from Clinical Samples.

Mycobacterium tuberculosis is the leading cause of death from bacterial infection. Improved rapid diagnosis and antimicrobial resistance determination, such as by whole-genome sequencing, are required. Our aim was to develop a simple, low-cost method of preparing DNA for sequencing direct from M. tuberculosis-positive clinical samples (without culture). Simultaneous sputum liquefaction, bacteria heat inactivation (99°C/30 min), and enrichment for mycobacteria DNA were achieved using an equal volume of thermo-protection buffer (4 M KCl, 0.05 M HEPES buffer, pH 7.5, 0.1% dithiothreitol [DTT]). The buffer emulated intracellular conditions found in hyperthermophiles, thus protecting DNA from rapid thermodegradation, which renders it a poor template for sequencing. Initial validation experiments employed mycobacteria DNA, either extracted or intracellular. Next, mock clinical samples (infection-negative human sputum spiked with 0 to 105Mycobacterium bovis BCG cells/ml) underwent liquefaction in thermo-protection buffer and heat inactivation. DNA was extracted and sequenced. Human DNA degraded faster than mycobacteria DNA, resulting in target enrichment. Four replicate experiments achieved M. tuberculosis detection at 101 BCG cells/ml, with 31 to 59 M. tuberculosis complex reads. Maximal genome coverage (>97% at 5× depth) occurred at 104 BCG cells/ml; >91% coverage (1× depth) occurred at 103 BCG cells/ml. Final validation employed M. tuberculosis-positive clinical samples (n = 20), revealing that initial sample volumes of ≥1 ml typically yielded higher mean depths of M. tuberculosis genome coverage, with an overall range of 0.55 to 81.02. A mean depth of 3 gave >96% 1-fold tuberculosis (TB) genome coverage (in 15/20 clinical samples). A mean depth of 15 achieved >99% 5-fold genome coverage (in 9/20 clinical samples). In summary, direct-from-sample sequencing of M. tuberculosis genomes was facilitated by a low-cost thermo-protection buffer.

Genomic surveillance of Escherichia coli and Klebsiella spp. in hospital sink drains and patients.

Escherichia coli and Klebsiella spp. are important human pathogens that cause a wide spectrum of clinical disease. In healthcare settings, sinks and other wastewater sites have been shown to be reservoirs of antimicrobial-resistant E. coli and Klebsiella spp., particularly in the context of outbreaks of resistant strains amongst patients. Without focusing exclusively on resistance markers or a clinical outbreak, we demonstrate that many hospital sink drains are abundantly and persistently colonized with diverse populations of E. coli, Klebsiella pneumoniae and Klebsiella oxytoca, including both antimicrobial-resistant and susceptible strains. Using whole-genome sequencing of 439 isolates, we show that environmental bacterial populations are largely structured by ward and sink, with only a handful of lineages, such as E. coli ST635, being widely distributed, suggesting different prevailing ecologies, which may vary as a result of different inputs and selection pressures. Whole-genome sequencing of 46 contemporaneous patient isolates identified one (2 %; 95 % CI 0.05-11 %) E. coli urine infection-associated isolate with high similarity to a prior sink isolate, suggesting that sinks may contribute to up to 10 % of infections caused by these organisms in patients on the ward over the same timeframe. Using metagenomics from 20 sink-timepoints, we show that sinks also harbour many clinically relevant antimicrobial resistance genes including blaCTX-M, blaSHV and mcr, and may act as niches for the exchange and amplification of these genes. Our study reinforces the potential role of sinks in contributing to Enterobacterales infection and antimicrobial resistance in hospital patients, something that could be amenable to intervention. This article contains data hosted by Microreact.

Multi-Label Random Forest Model for Tuberculosis Drug Resistance Classification and Mutation Ranking.

Resistance prediction and mutation ranking are important tasks in the analysis of Tuberculosis sequence data. Due to standard regimens for the use of first-line antibiotics, resistance co-occurrence, in which samples are resistant to multiple drugs, is common. Analysing all drugs simultaneously should therefore enable patterns reflecting resistance co-occurrence to be exploited for resistance prediction. Here, multi-label random forest (MLRF) models are compared with single-label random forest (SLRF) for both predicting phenotypic resistance from whole genome sequences and identifying important mutations for better prediction of four first-line drugs in a dataset of 13402 Mycobacterium tuberculosis isolates. Results confirmed that MLRFs can improve performance compared to conventional clinical methods (by 18.10%) and SLRFs (by 0.91%). In addition, we identified a list of candidate mutations that are important for resistance prediction or that are related to resistance co-occurrence. Moreover, we found that retraining our analysis to a subset of top-ranked mutations was sufficient to achieve satisfactory performance. The source code can be found at http://www.robots.ox.ac.uk/~davidc/code.php.

Genomic Epidemiology of Complex, Multispecies, Plasmid-Borne blaKPC Carbapenemase in Enterobacterales in the United Kingdom from 2009 to 2014.

Carbapenem resistance in Enterobacterales is a public health threat. Klebsiella pneumoniae carbapenemase (encoded by alleles of the blaKPC family) is one of the most common transmissible carbapenem resistance mechanisms worldwide. The dissemination of blaKPC historically has been associated with distinct K. pneumoniae lineages (clonal group 258 [CG258]), a particular plasmid family (pKpQIL), and a composite transposon (Tn4401). In the United Kingdom, blaKPC has represented a large-scale, persistent management challenge for some hospitals, particularly in North West England. The dissemination of blaKPC has evolved to be polyclonal and polyspecies, but the genetic mechanisms underpinning this evolution have not been elucidated in detail; this study used short-read whole-genome sequencing of 604 blaKPC-positive isolates (Illumina) and long-read assembly (PacBio)/polishing (Illumina) of 21 isolates for characterization. We observed the dissemination of blaKPC (predominantly blaKPC-2; 573/604 [95%] isolates) across eight species and more than 100 known sequence types. Although there was some variation at the transposon level (mostly Tn4401a, 584/604 [97%] isolates; predominantly with ATTGA-ATTGA target site duplications, 465/604 [77%] isolates), blaKPC spread appears to have been supported by highly fluid, modular exchange of larger genetic segments among plasmid populations dominated by IncFIB (580/604 isolates), IncFII (545/604 isolates), and IncR (252/604 isolates) replicons. The subset of reconstructed plasmid sequences (21 isolates, 77 plasmids) also highlighted modular exchange among non-blaKPC and blaKPC plasmids and the common presence of multiple replicons within blaKPC plasmid structures (>60%). The substantial genomic plasticity observed has important implications for our understanding of the epidemiology of transmissible carbapenem resistance in Enterobacterales for the implementation of adequate surveillance approaches and for control.

Metagenomic Nanopore Sequencing of Influenza Virus Direct from Clinical Respiratory Samples.

Influenza is a major global public health threat as a result of its highly pathogenic variants, large zoonotic reservoir, and pandemic potential. Metagenomic viral sequencing offers the potential for a diagnostic test for influenza virus which also provides insights on transmission, evolution, and drug resistance and simultaneously detects other viruses. We therefore set out to apply the Oxford Nanopore Technologies sequencing method to metagenomic sequencing of respiratory samples. We generated influenza virus reads down to a limit of detection of 102 to 103 genome copies/ml in pooled samples, observing a strong relationship between the viral titer and the proportion of influenza virus reads (P = 4.7 × 10-5). Applying our methods to clinical throat swabs, we generated influenza virus reads for 27/27 samples with mid-to-high viral titers (cycle threshold [CT ] values, <30) and 6/13 samples with low viral titers (CT values, 30 to 40). No false-positive reads were generated from 10 influenza virus-negative samples. Thus, Nanopore sequencing operated with 83% sensitivity (95% confidence interval [CI], 67 to 93%) and 100% specificity (95% CI, 69 to 100%) compared to the current diagnostic standard. Coverage of full-length virus was dependent on sample composition, being negatively influenced by increased host and bacterial reads. However, at high influenza virus titers, we were able to reconstruct >99% complete sequences for all eight gene segments. We also detected a human coronavirus coinfection in one clinical sample. While further optimization is required to improve sensitivity, this approach shows promise for the Nanopore platform to be used in the diagnosis and genetic analysis of influenza virus and other respiratory viruses.