Rifampicin resistant Mycobacterium tuberculosis in Vietnam, 2020-2022.

Objective: We conducted a descriptive analysis of multi-drug resistant tuberculosis (MDR-TB) in Vietnam's two largest cities, Hanoi and Ho Chi Minh city. Methods: All patients with rifampicin resistant tuberculosis were recruited from Hanoi and surrounding provinces between 2020 and 2022. Additional patients were recruited from Ho Chi Minh city over the same time period. Demographic data were recorded from all patients, and samples collected, cultured, whole genome sequenced and analysed for drug resistance mutations. Genomic susceptibility predictions were made on the basis of the World Health Organization's catalogue of mutations in Mycobacterium tuberculosis associated with drug resistance, version 2. Comparisons were made against phenotypic drug susceptibility test results where these were available. Multivariable logistic regression was used to assess risk factors for previous episodes of tuberculosis. Results: 233/265 sequenced isolates were of sufficient quality for analysis, 146 (63 %) from Ho Chi Minh City and 87 (37 %) from Hanoi. 198 (85 %) were lineage 2, 20 (9 %) were lineage 4, and 15 (6 %) were lineage 1. 17/211 (8 %) for whom HIV status was known were infected, and 109/214 (51 %) patients had had a previous episode of tuberculosis. The main risk factor for a previous episode was HIV infection (odds ratio 5.1 (95 % confidence interval 1.3-20.0); p = 0.021). Sensitivity for predicting first-line drug resistance from whole genome sequencing data was over 90 %, with the exception of pyrazinamide (85 %). For moxifloxacin and amikacin it was 50 % or less. Among rifampicin-resistant isolates, prevalence of resistance to each non-first-line drug was < 20 %. Conclusions: Drug resistance among most MDR-TB strains in Vietnam's two largest cities is confined largely to first-line drugs. Living with HIV is the main risk factor among patients with MDR-TB for having had a previous episode of tuberculosis.

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.

Hostile: accurate decontamination of microbial host sequences.

MOTIVATION: Microbial sequences generated from clinical samples are often contaminated with human host sequences that must be removed for ethical and legal reasons. Care must be taken to excise host sequences without inadvertently removing target microbial sequences to the detriment of downstream analyses such as variant calling and de novo assembly. RESULTS: To facilitate accurate host decontamination of both short and long sequencing reads, we developed Hostile, a tool capable of accurate host read removal using a laptop. We demonstrate that our approach removes at least 99.6% of real human reads and retains at least 99.989% of simulated bacterial reads. Using Hostile with a masked reference genome further increases bacterial read retention (≥99.997%) with negligible (≤0.001%) reduction in human read removal performance. Compared with an existing tool, Hostile removes 21%-23% more human short reads and 21-43 times fewer bacterial reads, typically in less time. AVAILABILITY AND IMPLEMENTATION: Hostile is implemented as an MIT-licensed Python package available from https://github.com/bede/hostile together with supplementary material.

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.

Persistence of SARS-CoV-2 antibodies over 18 months following infection: UK Biobank COVID-19 Serology Study.

BACKGROUND: Little is known about the persistence of antibodies after the first year following SARS-CoV-2 infection. We aimed to determine the proportion of individuals that maintain detectable levels of SARS-CoV-2 antibodies over an 18-month period following infection. METHODS: Population-based prospective study of 20 000 UK Biobank participants and their adult relatives recruited in May 2020. The proportion of SARS-CoV-2 cases testing positive for immunoglobulin G (IgG) antibodies against the spike protein (IgG-S), and the nucleocapsid protein (IgG-N), was calculated at varying intervals following infection. RESULTS: Overall, 20 195 participants were recruited. Their median age was 56 years (IQR 39-68), 56% were female and 88% were of white ethnicity. The proportion of SARS-CoV-2 cases with IgG-S antibodies following infection remained high (92%, 95% CI 90%-93%) at 6 months after infection. Levels of IgG-N antibodies following infection gradually decreased from 92% (95% CI 88%-95%) at 3 months to 72% (95% CI 70%-75%) at 18 months. There was no strong evidence of heterogeneity in antibody persistence by age, sex, ethnicity or socioeconomic deprivation. CONCLUSION: This study adds to the limited evidence on the long-term persistence of antibodies following SARS-CoV-2 infection, with likely implications for waning immunity following infection and the use of IgG-N in population surveys.

Protection against SARS-CoV-2 Omicron BA.4/5 variant following booster vaccination or breakthrough infection in the UK.

Following primary SARS-CoV-2 vaccination, whether boosters or breakthrough infections provide greater protection against SARS-CoV-2 infection is incompletely understood. Here we investigated SARS-CoV-2 antibody correlates of protection against new Omicron BA.4/5 (re-)infections and anti-spike IgG antibody trajectories after a third/booster vaccination or breakthrough infection following second vaccination in 154,149 adults ≥18 y from the United Kingdom general population. Higher antibody levels were associated with increased protection against Omicron BA.4/5 infection and breakthrough infections were associated with higher levels of protection at any given antibody level than boosters. Breakthrough infections generated similar antibody levels to boosters, and the subsequent antibody declines were slightly slower than after boosters. Together our findings show breakthrough infection provides longer-lasting protection against further infections than booster vaccinations. Our findings, considered alongside the risks of severe infection and long-term consequences of infection, have important implications for vaccine policy.

Penicillin Binding Protein Substitutions Cooccur with Fluoroquinolone Resistance in Epidemic Lineages of Multidrug-Resistant Clostridioides difficile.

Clostridioides difficile remains a key cause of healthcare-associated infection, with multidrug-resistant (MDR) lineages causing high-mortality (≥20%) outbreaks. Cephalosporin treatment is a long-established risk factor, and antimicrobial stewardship is a key control. A mechanism underlying raised cephalosporin MICs has not been identified in C. difficile, but among other species, this is often acquired via amino acid substitutions in cell wall transpeptidases (penicillin binding proteins [PBPs]). Here, we investigated five C. difficile transpeptidases (PBP1 to PBP5) for recent substitutions, associated cephalosporin MICs, and co-occurrence with fluoroquinolone resistance. Previously published genome assemblies (n = 7,096) were obtained, representing 16 geographically widespread lineages, including healthcare-associated ST1(027). Recent amino acid substitutions were found within PBP1 (n = 50) and PBP3 (n = 48), ranging from 1 to 10 substitutions per genome. ß-Lactam MICs were measured for closely related pairs of wild-type and PBP-substituted isolates separated by 20 to 273 single nucleotide polymorphisms (SNPs). Recombination-corrected phylogenies were constructed to date substitution acquisition. Key substitutions such as PBP3 V497L and PBP1 T674I/N/V emerged independently across multiple lineages. They were associated with extremely high cephalosporin MICs; 1 to 4 doubling dilutions >wild-type, up to 1,506 µg/mL. Substitution patterns varied by lineage and clade, showed geographic structure, and occurred post-1990, coincident with the gyrA and/or gyrB substitutions conferring fluoroquinolone resistance. In conclusion, recent PBP1 and PBP3 substitutions are associated with raised cephalosporin MICs in C. difficile. Their co-occurrence with fluoroquinolone resistance hinders attempts to understand the relative importance of these drugs in the dissemination of epidemic lineages. Further controlled studies of cephalosporin and fluoroquinolone stewardship are needed to determine their relative effectiveness in outbreak control. IMPORTANCE Fluoroquinolone and cephalosporin use in healthcare settings has triggered outbreaks of high-mortality, multidrug-resistant C. difficile infection. Here, we identify a mechanism associated with raised cephalosporin MICs in C. difficile comprising amino acid substitutions in two cell wall transpeptidase enzymes (penicillin binding proteins). The higher the number of substitutions, the greater the impact on phenotype. Dated phylogenies revealed that substitutions associated with raised cephalosporin and fluoroquinolone MICs were co-acquired immediately before clinically important outbreak strains emerged. PBP substitutions were geographically structured within genetic lineages, suggesting adaptation to local antimicrobial prescribing. Antimicrobial stewardship of cephalosporins and fluoroquinolones is an effective means of C. difficile outbreak control. Genetic changes associated with raised MIC may impart a "fitness cost" after antibiotic withdrawal. Our study therefore identifies a mechanism that may explain the contribution of cephalosporin stewardship to resolving outbreak conditions. However, due to the co-occurrence of raised cephalosporin MICs and fluoroquinolone resistance, further work is needed to determine the relative importance of each.

A longitudinal study reveals persistence of antimicrobial resistance on livestock farms is not due to antimicrobial usage alone.

Introduction: There are concerns that antimicrobial usage (AMU) is driving an increase in multi-drug resistant (MDR) bacteria so treatment of microbial infections is becoming harder in humans and animals. The aim of this study was to evaluate factors, including usage, that affect antimicrobial resistance (AMR) on farm over time. Methods: A population of 14 cattle, sheep and pig farms within a defined area of England were sampled three times over a year to collect data on AMR in faecal Enterobacterales flora; AMU; and husbandry or management practices. Ten pooled samples were collected at each visit, with each comprising of 10 pinches of fresh faeces. Up to 14 isolates per visit were whole genome sequenced to determine presence of AMR genes. Results: Sheep farms had very low AMU in comparison to the other species and very few sheep isolates were genotypically resistant at any time point. AMR genes were detected persistently across pig farms at all visits, even on farms with low AMU, whereas AMR bacteria was consistently lower on cattle farms than pigs, even for those with comparably high AMU. MDR bacteria was also more commonly detected on pig farms than any other livestock species. Discussion: The results may be explained by a complex combination of factors on pig farms including historic AMU; co-selection of AMR bacteria; variation in amounts of antimicrobials used between visits; potential persistence in environmental reservoirs of AMR bacteria; or importation of pigs with AMR microbiota from supplying farms. Pig farms may also be at increased risk of AMR due to the greater use of oral routes of group antimicrobial treatment, which were less targeted than cattle treatments; the latter mostly administered to individual animals. Also, farms which exhibited either increasing or decreasing trends of AMR across the study did not have corresponding trends in their AMU. Therefore, our results suggest that factors other than AMU on individual farms are important for persistence of AMR bacteria on farms, which may be operating at the farm and livestock species level.

Enterobacterales plasmid sharing amongst human bloodstream infections, livestock, wastewater, and waterway niches in Oxfordshire, UK.

Plasmids enable the dissemination of antimicrobial resistance (AMR) in common Enterobacterales pathogens, representing a major public health challenge. However, the extent of plasmid sharing and evolution between Enterobacterales causing human infections and other niches remains unclear, including the emergence of resistance plasmids. Dense, unselected sampling is essential to developing our understanding of plasmid epidemiology and designing appropriate interventions to limit the emergence and dissemination of plasmid-associated AMR. We established a geographically and temporally restricted collection of human bloodstream infection (BSI)-associated, livestock-associated (cattle, pig, poultry, and sheep faeces, farm soils) and wastewater treatment work (WwTW)-associated (influent, effluent, waterways upstream/downstream of effluent outlets) Enterobacterales. Isolates were collected between 2008 and 2020 from sites <60 km apart in Oxfordshire, UK. Pangenome analysis of plasmid clusters revealed shared 'backbones', with phylogenies suggesting an intertwined ecology where well-conserved plasmid backbones carry diverse accessory functions, including AMR genes. Many plasmid 'backbones' were seen across species and niches, raising the possibility that plasmid movement between these followed by rapid accessory gene change could be relatively common. Overall, the signature of identical plasmid sharing is likely to be a highly transient one, implying that plasmid movement might be occurring at greater rates than previously estimated, raising a challenge for future genomic One Health studies.

Genomic Sequencing from Sputum for Tuberculosis Disease Diagnosis, Lineage Determination, and Drug Susceptibility Prediction.

Universal access to drug susceptibility testing for newly diagnosed tuberculosis patients is recommended. Access to culture-based diagnostics remains limited, and targeted molecular assays are vulnerable to emerging resistance mutations. Improved protocols for direct-from-sputum Mycobacterium tuberculosis sequencing would accelerate access to comprehensive drug susceptibility testing and molecular typing. We assessed a thermo-protection buffer-based direct-from-sample M. tuberculosis whole-genome sequencing protocol. We prospectively analyzed 60 acid-fast bacilli smear-positive clinical sputum samples in India and Madagascar. A diversity of semiquantitative smear positivity-level samples were included. Sequencing was performed using Illumina and MinION (monoplex and multiplex) technologies. We measured the impact of bacterial inoculum and sequencing platforms on genomic read depth, drug susceptibility prediction performance, and typing accuracy. M. tuberculosis was identified by direct sputum sequencing in 45/51 samples using Illumina, 34/38 were identified using MinION-monoplex sequencing, and 20/24 were identified using MinION-multiplex sequencing. The fraction of M. tuberculosis reads from MinION sequencing was lower than from Illumina, but monoplexing grade 3+ samples on MinION produced higher read depth than Illumina (P < 0.05) and MinION multiplexing (P < 0.01). No significant differences in sensitivity and specificity of drug susceptibility predictions were seen across sequencing modalities or within each technology when stratified by smear grade. Illumina sequencing from sputum accurately identified 1/8 (rifampin) and 6/12 (isoniazid) resistant samples, compared to 2/3 (rifampin) and 3/6 (isoniazid) accurately identified with Nanopore monoplex. Lineage agreement levels between direct and culture-based sequencing were 85% (MinION-monoplex), 88% (Illumina), and 100% (MinION-multiplex). M. tuberculosis direct-from-sample whole-genome sequencing remains challenging. Improved and affordable sample treatment protocols are needed prior to clinical deployment.

Virus Detection and Identification in Minutes Using Single-Particle Imaging and Deep Learning.

The increasing frequency and magnitude of viral outbreaks in recent decades, epitomized by the COVID-19 pandemic, has resulted in an urgent need for rapid and sensitive diagnostic methods. Here, we present a methodology for virus detection and identification that uses a convolutional neural network to distinguish between microscopy images of fluorescently labeled intact particles of different viruses. Our assay achieves labeling, imaging, and virus identification in less than 5 min and does not require any lysis, purification, or amplification steps. The trained neural network was able to differentiate SARS-CoV-2 from negative clinical samples, as well as from other common respiratory pathogens such as influenza and seasonal human coronaviruses. We were also able to differentiate closely related strains of influenza, as well as SARS-CoV-2 variants. Additional and novel pathogens can easily be incorporated into the test through software updates, offering the potential to rapidly utilize the technology in future infectious disease outbreaks or pandemics. Single-particle imaging combined with deep learning therefore offers a promising alternative to traditional viral diagnostic and genomic sequencing methods and has the potential for significant impact.

High fluoroquinolone resistance proportions among multidrug-resistant tuberculosis driven by dominant L2 Mycobacterium tuberculosis clones in the Mumbai Metropolitan Region.

BACKGROUND: Multidrug-resistant (MDR) Mycobacterium tuberculosis complex (MTBC) strains are a serious health problem in India, also contributing to one-fourth of the global MDR tuberculosis (TB) burden. About 36% of the MDR MTBC strains are reported fluoroquinolone (FQ) resistant leading to high pre-extensively drug-resistant (pre-XDR) and XDR-TB (further resistance against bedaquiline and/or linezolid) rates. Still, factors driving the MDR/pre-XDR epidemic in India are not well defined. METHODS: In a retrospective study, we analyzed 1852 consecutive MTBC strains obtained from patients from a tertiary care hospital laboratory in Mumbai by whole genome sequencing (WGS). Univariate and multivariate statistics was used to investigate factors associated with pre-XDR. Core genome multi locus sequence typing, time scaled haplotypic density (THD) method and homoplasy analysis were used to analyze epidemiological success, and positive selection in different strain groups, respectively. RESULTS: In total, 1016 MTBC strains were MDR, out of which 703 (69.2%) were pre-XDR and 45 (4.4%) were XDR. Cluster rates were high among MDR (57.8%) and pre-XDR/XDR (79%) strains with three dominant L2 (Beijing) strain clusters (Cl 1-3) representing half of the pre-XDR and 40% of the XDR-TB cases. L2 strains were associated with pre-XDR/XDR-TB (P < 0.001) and, particularly Cl 1-3 strains, had high first-line and FQ resistance rates (81.6-90.6%). Epidemic success analysis using THD showed that L2 strains outperformed L1, L3, and L4 strains in short- and long-term time scales. More importantly, L2 MDR and MDR + strains had higher THD success indices than their not-MDR counterparts. Overall, compensatory mutation rates were highest in L2 strains and positive selection was detected in genes of L2 strains associated with drug tolerance (prpB and ppsA) and virulence (Rv2828c). Compensatory mutations in L2 strains were associated with a threefold increase of THD indices, suggesting improved transmissibility. CONCLUSIONS: Our data indicate a drastic increase of FQ resistance, as well as emerging bedaquiline resistance which endangers the success of newly endorsed MDR-TB treatment regimens. Rapid changes in treatment and control strategies are required to contain transmission of highly successful pre-XDR L2 strains in the Mumbai Metropolitan region but presumably also India-wide.

Mortality risks associated with empirical antibiotic activity in Escherichia coli bacteraemia: an analysis of electronic health records.

BACKGROUND: Reported bacteraemia outcomes following inactive empirical antibiotics (based on in vitro testing) are conflicting, potentially reflecting heterogeneity in causative species, MIC breakpoints defining resistance/susceptibility, and times to rescue therapy. METHODS: We investigated adult inpatients with Escherichia coli bacteraemia at Oxford University Hospitals, UK, from 4 February 2014 to 30 June 2021 who were receiving empirical amoxicillin/clavulanate with/without other antibiotics. We used Cox regression to analyse 30 day all-cause mortality by in vitro amoxicillin/clavulanate susceptibility (activity) using the EUCAST resistance breakpoint (>8/2 mg/L), categorical MIC, and a higher resistance breakpoint (>32/2 mg/L), adjusting for other antibiotic activity and confounders including comorbidities, vital signs and blood tests. RESULTS: A total of 1720 E. coli bacteraemias (1626 patients) were treated with empirical amoxicillin/clavulanate. Thirty-day mortality was 193/1400 (14%) for any active baseline therapy and 52/320 (16%) for inactive baseline therapy (P = 0.17). With EUCAST breakpoints, there was no evidence that mortality differed for inactive versus active amoxicillin/clavulanate [adjusted HR (aHR) = 1.27 (95% CI 0.83-1.93); P = 0.28], nor of an association with active aminoglycoside (P = 0.93) or other active antibiotics (P = 0.18). Considering categorical amoxicillin/clavulanate MIC, MICs > 32/2 mg/L were associated with mortality [aHR = 1.85 versus MIC = 2/2 mg/L (95% CI 0.99-3.73); P = 0.054]. A higher resistance breakpoint (>32/2 mg/L) was independently associated with higher mortality [aHR = 1.82 (95% CI 1.07-3.10); P = 0.027], as were MICs > 32/2 mg/L with active empirical aminoglycosides [aHR = 2.34 (95% CI 1.40-3.89); P = 0.001], but not MICs > 32/2 mg/L with active non-aminoglycoside antibiotic(s) [aHR = 0.87 (95% CI 0.40-1.89); P = 0.72]. CONCLUSIONS: We found no evidence that EUCAST-defined amoxicillin/clavulanate resistance was associated with increased mortality, but a higher resistance breakpoint (MIC > 32/2 mg/L) was. Additional active baseline non-aminoglycoside antibiotics attenuated amoxicillin/clavulanate resistance-associated mortality, but aminoglycosides did not. Granular phenotyping and comparison with clinical outcomes may improve AMR breakpoints.

SARS-CoV-2 antibody trajectories after a single COVID-19 vaccination with and without prior infection.

Given high SARS-CoV-2 incidence, coupled with slow and inequitable vaccine roll-out in many settings, there is a need for evidence to underpin optimum vaccine deployment, aiming to maximise global population immunity. We evaluate whether a single vaccination in individuals who have already been infected with SARS-CoV-2 generates similar initial and subsequent antibody responses to two vaccinations in those without prior infection. We compared anti-spike IgG antibody responses after a single vaccination with ChAdOx1, BNT162b2, or mRNA-1273 SARS-CoV-2 vaccines in the COVID-19 Infection Survey in the UK general population. In 100,849 adults median (50 (IQR: 37-63) years) receiving at least one vaccination, 13,404 (13.3%) had serological/PCR evidence of prior infection. Prior infection significantly boosted antibody responses, producing higher peak levels and/or longer half-lives after one dose of all three vaccines than those without prior infection receiving one or two vaccinations. In those with prior infection, the median time above the positivity threshold was >1 year after the first vaccination. Single-dose vaccination targeted to those previously infected may provide at least as good protection to two-dose vaccination among those without previous infection.

Catwalk: identifying closely related sequences in large microbial sequence databases.

There is a need to identify microbial sequences that may form part of transmission chains, or that may represent importations across national boundaries, amidst large numbers of SARS-CoV-2 and other bacterial or viral sequences. Reference-based compression is a sequence analysis technique that allows both a compact storage of sequence data and comparisons between sequences. Published implementations of the approach are being challenged by the large sample collections now being generated. Our aim was to develop a fast software detecting highly similar sequences in large collections of microbial genomes, including millions of SARS-CoV-2 genomes. To do so, we developed Catwalk, a tool that bypasses bottlenecks in the generation, comparison and in-memory storage of microbial genomes generated by reference mapping. It is a compiled solution, coded in Nim to increase performance. It can be accessed via command line, rest api or web server interfaces. We tested Catwalk using both SARS-CoV-2 and Mycobacterium tuberculosis genomes generated by prospective public-health sequencing programmes. Pairwise sequence comparisons, using clinically relevant similarity cut-offs, took about 0.39 and 0.66 µs, respectively; in 1 s, between 1 and 2 million sequences can be searched. Catwalk operates about 1700 times faster than, and uses about 8 % of the RAM of, a Python reference-based compression and comparison tool in current use for outbreak detection. Catwalk can rapidly identify close relatives of a SARS-CoV-2 or M. tuberculosis genome amidst millions of samples.

The bacteriology of pleural infection (TORPIDS): an exploratory metagenomics analysis through next generation sequencing.

BACKGROUND: Pleural infection is a common and severe disease with high morbidity and mortality worldwide. The knowledge of pleural infection bacteriology remains incomplete, as pathogen detection methods based on culture have insufficient sensitivity and are biased to selected microbes. We designed a study with the aim to discover and investigate the total microbiome of pleural infection and assess the correlation between bacterial patterns and 1-year survival of patients. METHODS: We assessed 243 pleural fluid samples from the PILOT study, a prospective observational study on pleural infection, with 16S rRNA next generation sequencing. 20 pleural fluid samples from patients with pleural effusion due to a non-infectious cause and ten PCR-grade water samples were used as controls. Downstream analysis was done with the DADA2 pipeline. We applied multivariate Cox regression analyses to investigate the association between bacterial patterns and 1-year survival of patients with pleural infection. FINDINGS: Pleural infection was predominately polymicrobial (192 [79%] of 243 samples), with diverse bacterial frequencies observed in monomicrobial and polymicrobial disease and in both community-acquired and hospital-acquired infection. Mixed anaerobes and other Gram-negative bacteria predominated in community-acquired polymicrobial infection whereas Streptococcus pneumoniae prevailed in monomicrobial cases. The presence of anaerobes (hazard ratio 0·46, 95% CI 0·24-0·86, p=0·015) or bacteria of the Streptococcus anginosus group (0·43, 0·19-0·97, p=0·043) was associated with better patient survival, whereas the presence (5·80, 2·37-14·21, p<0·0001) or dominance (3·97, 1·20-13·08, p=0·024) of Staphylococcus aureus was linked with lower survival. Moreover, dominance of Enterobacteriaceae was associated with higher risk of death (2·26, 1·03-4·93, p=0·041). INTERPRETATION: Pleural infection is a predominantly polymicrobial infection, explaining the requirement for broad spectrum antibiotic cover in most individuals. High mortality infection associated with S aureus and Enterobacteriaceae favours more aggressive, with a narrower spectrum, antibiotic strategies. FUNDING: UK Medical Research Council, National Institute for Health Research Oxford Biomedical Research Centre, Wellcome Trust, Oxfordshire Health Services Research Committee, Chinese Academy of Medical Sciences, and John Fell Fund.

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.

The 2021 WHO catalogue of Mycobacterium tuberculosis complex mutations associated with drug resistance: A genotypic analysis.

Background: Molecular diagnostics are considered the most promising route to achieving rapid, universal drug susceptibility testing for Mycobacterium tuberculosiscomplex (MTBC). We aimed to generate a WHO endorsed catalogue of mutations to serve as a global standard for interpreting molecular information for drug resistance prediction. Methods: A candidate gene approach was used to identify mutations as associated with resistance, or consistent with susceptibility, for 13 WHO endorsed anti-tuberculosis drugs. 38,215 MTBC isolates with paired whole-genome sequencing and phenotypic drug susceptibility testing data were amassed from 45 countries. For each mutation, a contingency table of binary phenotypes and presence or absence of the mutation computed positive predictive value, and Fisher's exact tests generated odds ratios and Benjamini-Hochberg corrected p-values. Mutations were graded as Associated with Resistance if present in at least 5 isolates, if the odds ratio was >1 with a statistically significant corrected p-value, and if the lower bound of the 95% confidence interval on the positive predictive value for phenotypic resistance was >25%. A series of expert rules were applied for final confidence grading of each mutation. Findings: 15,667 associations were computed for 13,211 unique mutations linked to one or more drugs. 1,149/15,667 (7·3%) mutations were classified as associated with phenotypic resistance and 107/15,667 (0·7%) were deemed consistent with susceptibility. For rifampicin, isoniazid, ethambutol, fluoroquinolones, and streptomycin, the mutations' pooled sensitivity was >80%. Specificity was over 95% for all drugs except ethionamide (91·4%), moxifloxacin (91·6%) and ethambutol (93·3%). Only two resistance mutations were classified for bedaquiline, delamanid, clofazimine, and linezolid as prevalence of phenotypic resistance was low for these drugs. Interpretation: This first WHO endorsed catalogue of molecular targets for MTBC drug susceptibility testing provides a global standard for resistance interpretation. Its existence should encourage the implementation of molecular diagnostics by National Tuberculosis Programmes. Funding: UNITAID, Wellcome, MRC, BMGF.

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.