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.

Dissemination and persistence of antimicrobial resistance (AMR) along the wastewater-river continuum.

Antimicrobial resistance (AMR) is a global health hazard. Although clinical and agricultural environments are well-established contributors to the evolution and dissemination of AMR, research on wastewater treatment works (WwTWs) has highlighted their potential role as disseminators of AMR in freshwater environments. Using metagenomic sequencing and analysis, we investigated the changes in resistomes and associated mobile genetic elements within untreated wastewater influents and treated effluents of five WwTWs, and sediments collected from corresponding river environments in Oxfordshire, UK, across three seasonal periods within a year. Our analysis demonstrated a high diversity and abundance of antimicrobial resistance genes (ARGs) in untreated wastewater influents, reflecting the varied anthropogenic and environmental origins of wastewater. WwTWs effectively reduced AMR in the final effluent, with an average 87 % reduction in normalised ARG abundance and an average 63 % reduction in richness. However, wastewater effluents significantly impacted the antimicrobial resistome of the receiving rivers, with an average 543 % increase in ARG abundance and a 164 % increase in richness from upstream sediments to downstream sediments. The normalised abundance of the human gut-associated bacteriophage crAssphage was highly associated with both ARG abundance and richness. We observed seasonal variation in the resistome of raw influent which was not found in the effluent-receiving sediments. We illustrate the potential of WwTWs as focal points for disseminating ARGs and resistance-selecting chemicals, contributing to the elevation of environmental AMR. Our study emphasises the need for a comprehensive understanding of the anthropogenic impacts on AMR evolution and dissemination in wastewater and river environments, informing efforts to mitigate this growing public health crisis.

Improving the representativeness of UK's national COVID-19 Infection Survey through spatio-temporal regression and post-stratification.

Population-representative estimates of SARS-CoV-2 infection prevalence and antibody levels in specific geographic areas at different time points are needed to optimise policy responses. However, even population-wide surveys are potentially impacted by biases arising from differences in participation rates across key groups. Here, we used spatio-temporal regression and post-stratification models to UK's national COVID-19 Infection Survey (CIS) to obtain representative estimates of PCR positivity (6,496,052 tests) and antibody prevalence (1,941,333 tests) for different regions, ages and ethnicities (7-December-2020 to 4-May-2022). Not accounting for vaccination status through post-stratification led to small underestimation of PCR positivity, but more substantial overestimations of antibody levels in the population (up to 21 percentage points), particularly in groups with low vaccine uptake in the general population. There was marked variation in the relative contribution of different areas and age-groups to each wave. Future analyses of infectious disease surveys should take into account major drivers of outcomes of interest that may also influence participation, with vaccination being an important factor to consider.

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.

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.