Evaluation of methods for detecting human reads in microbial sequencing datasets

Sequencing data from host-associated microbes can often be contaminated by the body of the investigator or research subject. Human DNA is typically removed from microbial reads either by subtractive alignment (dropping all reads that map to the human genome) or by using a read classification tool to predict those of human origin, and then discarding them. To inform best practice guidelines, we benchmarked eight alignment-based and two classification-based methods of human read detection using simulated data from 10 clinically prevalent bacteria and three viruses, into which contaminating human reads had been added. While the majority of methods successfully detected >99 % of the human reads, they were distinguishable by variance. The most precise methods, with negligible variance, were Bowtie2 and SNAP, both of which misidentified few, if any, bacterial reads (and no viral reads) as human. While correctly detecting a similar number of human reads, methods based on taxonomic classification, such as Kraken2 and Centrifuge, could misclassify bacterial reads as human, although the extent of this was species-specific. Among the most sensitive methods of human read detection was BWA, although this also made the greatest number of false positive classifications. Across all methods, the set of human reads not identified as such, although often representing 300 bp) bacterial reads, the highest performing approaches were classification-based, using Kraken2 or Centrifuge. For shorter (c. 150 bp) bacterial reads, combining multiple methods of human read detection maximized the recovery of human reads from contaminated short read datasets without being compromised by false positives. A particularly high-performance approach with shorter bacterial reads was a two-stage classification using Bowtie2 followed by SNAP. Using this approach, we re-examined 11 577 publicly archived bacterial read sets for hitherto undetected human contamination. We were able to extract a sufficient number of reads to call known human SNPs, including those with clinical significance, in 6 % of the samples. These results show that phenotypically distinct human sequence is detectable in publicly archived microbial read datasets

Extended Sequence Typing of Campylobacter spp., United Kingdom

Supplementing Campylobacter spp. multilocus sequence typing with nucleotide sequence typing of 3 antigen genes increased the discriminatory index achieved from 0.975 to 0.992 among 620 clinical isolates from Oxfordshire, United Kingdom. This enhanced typing scheme enabled identification of clusters and retained data required for long-range epidemiologic comparisons of isolates

Same-day diagnostic and surveillance data for tuberculosis via whole genome sequencing of direct respiratory samples

Routine full characterization of Mycobacterium tuberculosis (TB) is culture-based, taking many weeks. Whole-genome sequencing (WGS) can generate antibiotic susceptibility profiles to inform treatment, augmented with strain information for global surveillance; such data could be transformative if provided at or near point of care. We demonstrate a low-cost DNA extraction method for TB WGS direct from patient samples. We initially evaluated the method using the Illumina MiSeq sequencer (40 smear-positive respiratory samples, obtained after routine clinical testing, and 27 matched liquid cultures). M. tuberculosis was identified in all 39 samples from which DNA was successfully extracted. Sufficient data for antibiotic susceptibility prediction was obtained from 24 (62%) samples; all results were concordant with reference laboratory phenotypes. Phylogenetic placement was concordant between direct and cultured samples. Using an 70 Illumina MiSeq/MiniSeq the workflow from patient sample to results can be completed in 44/16 hours at a reagent cost of £96/£198 per sample. We then employed a non-specific PCR-based library preparation method for sequencing on an Oxford Nanopore Technologies MinION sequencer. We applied this to cultured Mycobacterium bovis BCG strain (BCG), and to combined culture negative sputum DNA and BCG DNA. For flowcell version R9.4, the estimated turnaround time from patient to identification of BCG, detection of pyrazinamide resistance, and phylogenetic placement was 7.5 hours, with full susceptibility results 5 hours later. Antibiotic susceptibility predictions were fully concordant. A critical advantage of the MinION is the ability to continue sequencing until sufficient coverage is obtained, providing a potential solution to the problem of variable amounts of M. tuberculosis in direct samples

Population-level faecal metagenomic profiling as a tool to predict antimicrobial resistance in Enterobacterales isolates causing invasive infections: an exploratory study across Cambodia, Kenya, and the UK

Background: Antimicrobial resistance (AMR) in Enterobacterales is a global health threat. Capacity for individual-level surveillance remains limited in many countries, whilst population-level surveillance approaches could inform empiric antibiotic treatment guidelines. Methods: In this exploratory study, a novel approach to population-level prediction of AMR in Enterobacterales clinical isolates using metagenomic (Illumina) profiling of pooled DNA extracts from human faecal samples was developed and tested. Taxonomic and AMR gene profiles were used to derive taxonomy-adjusted population-level AMR metrics. Bayesian modelling, and model comparison based on cross-validation, were used to evaluate the capacity of each metric to predict the number of resistant Enterobacterales invasive infections at a population-level, using available bloodstream/cerebrospinal fluid infection data. Findings: Population metagenomes comprised samples from 177, 157, and 156 individuals in Kenya, the UK, and Cambodia, respectively, collected between September 2014 and April 2016. Clinical data from independent populations included 910, 3356 and 197 bacterial isolates from blood/cerebrospinal fluid infections in Kenya, the UK and Cambodia, respectively (samples collected between January 2010 and May 2017). Enterobacterales were common colonisers and pathogens, and faecal taxonomic/AMR gene distributions and proportions of antimicrobial-resistant Enterobacterales infections differed by setting. A model including terms reflecting the metagenomic abundance of the commonest clinical Enterobacterales species, and of AMR genes known to either increase the minimum inhibitory concentration (MIC) or confer clinically-relevant resistance, had a higher predictive performance in determining population-level resistance in clinical Enterobacterales isolates compared to models considering only AMR gene information, only taxonomic information, or an intercept-only baseline model (difference in expected log predictive density compared to best model, estimated using leave-one-out cross-validation: intercept-only model = -223 [95% credible interval (CI): -330,-116]; model considering only AMR gene information = -186 [95% CI: -281,-91]; model considering only taxonomic information = -151 [95% CI: -232,-69]). Interpretation: Whilst our findings are exploratory and require validation, intermittent metagenomics of pooled samples could represent an effective approach for AMR surveillance and to predict population-level AMR in clinical isolates, complementary to ongoing development of laboratory infrastructures processing individual samples

A cluster of multidrug-resistant Mycobacterium tuberculosis among patients arriving in Europe from the Horn of Africa: a molecular epidemiological study

SummaryBackground The risk of tuberculosis outbreaks among people fleeing hardship for refuge in Europe is heightened. We describe the cross-border European response to an outbreak of multidrug-resistant tuberculosis among patients from the Horn of Africa and Sudan. Methods On April 29 and May 30, 2016, the Swiss and German National Mycobacterial Reference Laboratories independently triggered an outbreak investigation after four patients were diagnosed with multidrug-resistant tuberculosis. In this molecular epidemiological study, we prospectively defined outbreak cases with 24-locus mycobacterial interspersed repetitive unit-variable number tandem repeat (MIRU-VNTR) profiles; phenotypic resistance to isoniazid, rifampicin, ethambutol, pyrazinamide, and capreomycin; and corresponding drug resistance mutations. We whole-genome sequenced all Mycobacterium tuberculosis isolates and clustered them using a threshold of five single nucleotide polymorphisms (SNPs). We collated epidemiological data from host countries from the European Centre for Disease Prevention and Control. Findings Between Feb 12, 2016, and April 19, 2017, 29 patients were diagnosed with multidrug-resistant tuberculosis in seven European countries. All originated from the Horn of Africa or Sudan, with all isolates two SNPs or fewer apart. 22 (76%) patients reported their travel routes, with clear spatiotemporal overlap between routes. We identified a further 29 MIRU-VNTR-linked cases from the Horn of Africa that predated the outbreak, but all were more than five SNPs from the outbreak. However all 58 isolates shared a capreomycin resistance-associated tlyA mutation. Interpretation Our data suggest that source cases are linked to an M tuberculosis clone circulating in northern Somalia or Djibouti and that transmission probably occurred en route before arrival in Europe. We hypothesise that the shared mutation of tlyA is a drug resistance mutation and phylogenetic marker, the first of its kind in M tuberculosis sensu stricto. Funding The Swiss Federal Office of Public Health, the University of Zurich, the Wellcome Trust, National Institute for Health Research (NIHR) Oxford Biomedical Research Centre (BRC), the Medical Research Council, BELTA-TBnet, the European Union, the German Center for Infection Research, and Leibniz Science Campus Evolutionary Medicine of the Lung (EvoLUNG)

The relative importance of large problems far away versus small problems closer to home:Insights into limiting the spread of antimicrobial resistance in England

Abstract Background To combat the spread of antimicrobial resistance (AMR), hospitals are advised to screen high-risk patients for carriage of antibiotic-resistant bacteria on admission. This often includes patients previously admitted to hospitals with a high AMR prevalence. However, the ability of such a strategy to identify introductions (and hence prevent onward transmission) is unclear, as it depends on AMR prevalence in each hospital, the number of patients moving between hospitals, and the number of hospitals considered ‘high risk’. Methods We tracked patient movements using data from the National Health Service of England Hospital Episode Statistics and estimated differences in regional AMR prevalences using, as an exemplar, data collected through the national reference laboratory service of Public Health England on carbapenemase-producing Enterobacteriaceae (CPE) from 2008 to 2014. Combining these datasets, we calculated expected CPE introductions into hospitals from across the hospital network to assess the effectiveness of admission screening based on defining high-prevalence hospitals as high risk. Results Based on numbers of exchanged patients, the English hospital network can be divided into 14 referral regions. England saw a sharp increase in numbers of CPE isolates referred to the national reference laboratory over 7 years, from 26 isolates in 2008 to 1649 in 2014. Large regional differences in numbers of confirmed CPE isolates overlapped with regional structuring of patient movements between hospitals. However, despite these large differences in prevalence between regions, we estimated that hospitals received only a small proportion (1.8%) of CPE-colonised patients from hospitals outside their own region, which decreased over time. Conclusions In contrast to the focus on import screening based on assigning a few hospitals as ‘high risk’, patient transfers between hospitals with small AMR problems in the same region often pose a larger absolute threat than patient transfers from hospitals in other regions with large problems, even if the prevalence in other regions is orders of magnitude higher. Because the difference in numbers of exchanged patients, between and within regions, was mostly larger than the difference in CPE prevalence, it would be more effective for hospitals to focus on their own populations or region to inform control efforts rather than focussing on problems elsewhere

Diverse sources of C. difficile infection identified on whole-genome sequencing

BACKGROUND It has been thought that Clostridium difficile infection is transmitted predominantly within health care settings. However, endemic spread has hampered identification of precise sources of infection and the assessment of the efficacy of interventions. METHODS From September 2007 through March 2011, we performed whole-genome sequencing on isolates obtained from all symptomatic patients with C. difficile infection identified in health care settings or in the community in Oxfordshire, United Kingdom. We compared single-nucleotide variants (SNVs) between the isolates, using C. difficile evolution rates estimated on the basis of the first and last samples obtained from each of 145 patients, with 0 to 2 SNVs expected between transmitted isolates obtained less than 124 days apart, on the basis of a 95% prediction interval. We then identified plausible epidemiologic links among genetically related cases from data on hospital admissions and community location. RESULTS Of 1250 C. difficile cases that were evaluated, 1223 (98%) were successfully sequenced. In a comparison of 957 samples obtained from April 2008 through March 2011 with those obtained from September 2007 onward, a total of 333 isolates (35%) had no more than 2 SNVs from at least 1 earlier case, and 428 isolates (45%) had more than 10 SNVs from all previous cases. Reductions in incidence over time were similar in the two groups, a finding that suggests an effect of interventions targeting the transition from exposure to disease. Of the 333 patients with no more than 2 SNVs (consistent with transmission), 126 patients (38%) had close hospital contact with another patient, and 120 patients (36%) had no hospital or community contact with another patient. Distinct subtypes of infection continued to be identified throughout the study, which suggests a considerable reservoir of C. difficile. CONCLUSIONS Over a 3-year period, 45% of C. difficile cases in Oxfordshire were genetically distinct from all previous cases. Genetically diverse sources, in addition to symptomatic patients, play a major part in C. difficile transmission. (Funded by the U.K. Clinical Research Collaboration Translational Infection Research Initiative and others.

Insidious risk of severe mycobacterium chimaera infection in cardiac surgery patients

Background An urgent UK investigation was launched to assess risk of invasive Mycobacterium chimaera infection in cardiothoracic surgery and a possible association with cardiopulmonary bypass heater-cooler units following alerts in Switzerland and The Netherlands. Methods Parallel investigations were pursued: (1) identification of cardiopulmonary bypass–associated M. chimaera infection through national laboratory and hospital admissions data linkage; (2) cohort study to assess patient risk; (3) microbiological and aerobiological investigations of heater-coolers in situ and under controlled laboratory conditions; and (4) whole-genome sequencing of clinical and environmental isolates. Results Eighteen probable cases of cardiopulmonary bypass–associated M. chimaera infection were identified; all except one occurred in adults. Patients had undergone valve replacement in 11 hospitals between 2007 and 2015, a median of 19 months prior to onset (range, 3 months to 5 years). Risk to patients increased after 2010 from Conclusions We identified low but escalating risk of severe M. chimaera infection associated with heater-coolers with cases in a quarter of cardiothoracic centers. Our investigations strengthen etiological evidence for the role of heater-coolers in transmission and raise the possibility of an ongoing, international point-source outbreak. Active management of heater-coolers and heightened clinical awareness are imperative given the consequences of infection.</p

Effect of Amoxicillin dose and treatment duration on the need for antibiotic re-treatment in children with Community-Acquired Pneumonia: The CAP-IT randomized clinical trial

Importance: The optimal dose and duration of oral amoxicillin for children with community-acquired pneumonia (CAP) are unclear. Objective: To determine whether lower-dose amoxicillin is noninferior to higher dose and whether 3-day treatment is noninferior to 7 days. Design, Setting, and Participants: Multicenter, randomized, 2 × 2 factorial noninferiority trial enrolling 824 children, aged 6 months and older, with clinically diagnosed CAP, treated with amoxicillin on discharge from emergency departments and inpatient wards of 28 hospitals in the UK and 1 in Ireland between February 2017 and April 2019, with last trial visit on May 21, 2019. Interventions: Children were randomized 1:1 to receive oral amoxicillin at a lower dose (35-50 mg/kg/d; n = 410) or higher dose (70-90 mg/kg/d; n = 404), for a shorter duration (3 days; n = 413) or a longer duration (7 days; n = 401). Main Outcomes and Measures: The primary outcome was clinically indicated antibiotic re-treatment for respiratory infection within 28 days after randomization. The noninferiority margin was 8%. Secondary outcomes included severity/duration of 9 parent-reported CAP symptoms, 3 antibiotic-related adverse events, and phenotypic resistance in colonizing Streptococcus pneumoniae isolates. Results: Of 824 participants randomized into 1 of the 4 groups, 814 received at least 1 dose of trial medication (median [IQR] age, 2.5 years [1.6-2.7]; 421 [52%] males and 393 [48%] females), and the primary outcome was available for 789 (97%). For lower vs higher dose, the primary outcome occurred in 12.6% with lower dose vs 12.4% with higher dose (difference, 0.2% [1-sided 95% CI -∞ to 4.0%]), and in 12.5% with 3-day treatment vs 12.5% with 7-day treatment (difference, 0.1% [1-sided 95% CI -∞ to 3.9]). Both groups demonstrated noninferiority with no significant interaction between dose and duration (P =.63). Of the 14 prespecified secondary end points, the only significant differences were 3-day vs 7-day treatment for cough duration (median 12 days vs 10 days; hazard ratio [HR], 1.2 [95% CI, 1.0 to 1.4]; P =.04) and sleep disturbed by cough (median, 4 days vs 4 days; HR, 1.2 [95% CI, 1.0 to 1.4]; P =.03). Among the subgroup of children with severe CAP, the primary end point occurred in 17.3% of lower-dose recipients vs 13.5% of higher-dose recipients (difference, 3.8% [1-sided 95% CI, -∞ to10%]; P value for interaction =.18) and in 16.0% with 3-day treatment vs 14.8% with 7-day treatment (difference, 1.2% [1-sided 95% CI, -∞ to 7.4%]; P value for interaction =.73). Conclusions and Relevance: Among children with CAP discharged from an emergency department or hospital ward (within 48 hours), lower-dose outpatient oral amoxicillin was noninferior to higher dose, and 3-day duration was noninferior to 7 days, with regard to need for antibiotic re-treatment. However, disease severity, treatment setting, prior antibiotics received, and acceptability of the noninferiority margin require consideration when interpreting the findings. Trial Registration: ISRCTN Identifier: ISRCTN76888927