Comparison of long-read sequencing technologies in the hybrid assembly of complex bacterial genomes.

Illumina sequencing allows rapid, cheap and accurate whole genome bacterial analyses, but short reads (<300 bp) do not usually enable complete genome assembly. Long-read sequencing greatly assists with resolving complex bacterial genomes, particularly when combined with short-read Illumina data (hybrid assembly). However, it is not clear how different long-read sequencing methods affect hybrid assembly accuracy. Relative automation of the assembly process is also crucial to facilitating high-throughput complete bacterial genome reconstruction, avoiding multiple bespoke filtering and data manipulation steps. In this study, we compared hybrid assemblies for 20 bacterial isolates, including two reference strains, using Illumina sequencing and long reads from either Oxford Nanopore Technologies (ONT) or SMRT Pacific Biosciences (PacBio) sequencing platforms. We chose isolates from the family Enterobacteriaceae, as these frequently have highly plastic, repetitive genetic structures, and complete genome reconstruction for these species is relevant for a precise understanding of the epidemiology of antimicrobial resistance. We de novo assembled genomes using the hybrid assembler Unicycler and compared different read processing strategies, as well as comparing to long-read-only assembly with Flye followed by short-read polishing with Pilon. Hybrid assembly with either PacBio or ONT reads facilitated high-quality genome reconstruction, and was superior to the long-read assembly and polishing approach evaluated with respect to accuracy and completeness. Combining ONT and Illumina reads fully resolved most genomes without additional manual steps, and at a lower consumables cost per isolate in our setting. Automated hybrid assembly is a powerful tool for complete and accurate bacterial genome assembly.

Evaluation of Whole-Genome Sequencing for Mycobacterial Species Identification and Drug Susceptibility Testing in a Clinical Setting: a Large-Scale Prospective Assessment of Performance against Line Probe Assays and Phenotyping.

Use of whole-genome sequencing (WGS) for routine mycobacterial species identification and drug susceptibility testing (DST) is becoming a reality. We compared the performances of WGS and standard laboratory workflows prospectively, by parallel processing at a major mycobacterial reference service over the course of 1 year, for species identification, first-line Mycobacterium tuberculosis resistance prediction, and turnaround time. Among 2,039 isolates with line probe assay results for species identification, 74 (3.6%) failed sequencing or WGS species identification. Excluding these isolates, clinically important species were identified for 1,902 isolates, of which 1,825 (96.0%) were identified as the same species by WGS and the line probe assay. A total of 2,157 line probe test results for detection of resistance to the first-line drugs isoniazid and rifampin were available for 728 M. tuberculosis complex isolates. Excluding 216 (10.0%) cases where there were insufficient sequencing data for WGS to make a prediction, overall concordance was 99.3% (95% confidence interval [CI], 98.9 to 99.6%), sensitivity was 97.6% (91.7 to 99.7%), and specificity was 99.5% (99.0 to 99.7%). A total of 2,982 phenotypic DST results were available for 777 M. tuberculosis complex isolates. Of these, 356 (11.9%) had no WGS comparator due to insufficient sequencing data, and in 154 (5.2%) cases the WGS prediction was indeterminate due to discovery of novel, previously uncharacterized mutations. Excluding these data, overall concordance was 99.2% (98.7 to 99.5%), sensitivity was 94.2% (88.4 to 97.6%), and specificity was 99.4% (99.0 to 99.7%). Median processing times for the routine laboratory tests versus WGS were similar overall, i.e., 20 days (interquartile range [IQR], 15 to 31 days) and 21 days (15 to 29 days), respectively (P = 0.41). In conclusion, WGS predicts species and drug susceptibility with great accuracy, but work is needed to increase the proportion of predictions made.

Resolving plasmid structures in Enterobacteriaceae using the MinION nanopore sequencer: assessment of MinION and MinION/Illumina hybrid data assembly approaches.

This study aimed to assess the feasibility of using the Oxford Nanopore Technologies (ONT) MinION long-read sequencer in reconstructing fully closed plasmid sequences from eight Enterobacteriaceae isolates of six different species with plasmid populations of varying complexity. Species represented were Escherichia coli, Klebsiella pneumoniae, Citrobacter freundii, Enterobacter cloacae, Serratia marcescens and Klebsiella oxytoca, with plasmid populations ranging from 1-11 plasmids with sizes of 2-330 kb. Isolates were sequenced using Illumina (short-read) and ONT's MinION (long-read) platforms, and compared with fully resolved PacBio (long-read) sequence assemblies for the same isolates. We compared the performance of different assembly approaches including SPAdes, plasmidSPAdes, hybridSPAdes, Canu, Canu+Pilon (canuPilon) and npScarf in recovering the plasmid structures of these isolates by comparing with the gold-standard PacBio reference sequences. Overall, canuPilon provided consistently good quality assemblies both in terms of assembly statistics (N50, number of contigs) and assembly accuracy [presence of single nucleotide polymorphisms (SNPs)/indels with respect to the reference sequence]. For plasmid reconstruction, Canu recovered 70 % of the plasmids in complete contigs, and combining three assembly approaches (Canu or canuPilon, hybridSPAdes and plasmidSPAdes) resulted in a total 78 % recovery rate for all the plasmids. The analysis demonstrated the potential of using MinION sequencing technology to resolve important plasmid structures in Enterobacteriaceae species independent of and in conjunction with Illumina sequencing data. A consensus assembly derived from several assembly approaches could present significant benefit in accurately resolving the greatest number of plasmid structures.

Mortality risks associated with emergency admissions during weekends and public holidays: an analysis of electronic health records.

BACKGROUND: Weekend hospital admission is associated with increased mortality, but the contributions of varying illness severity and admission time to this weekend effect remain unexplored. METHODS: We analysed unselected emergency admissions to four Oxford University National Health Service hospitals in the UK from Jan 1, 2006, to Dec 31, 2014. The primary outcome was death within 30 days of admission (in or out of hospital), analysed using Cox models measuring time from admission. The primary exposure was day of the week of admission. We adjusted for multiple confounders including demographics, comorbidities, and admission characteristics, incorporating non-linearity and interactions. Models then considered the effect of adjusting for 15 common haematology and biochemistry test results or proxies for hospital workload. FINDINGS: 257 596 individuals underwent 503 938 emergency admissions. 18 313 (4·7%) patients admitted as weekday energency admissions and 6070 (5·1%) patients admitted as weekend emergency admissions died within 30 days (p<0·0001). 9347 individuals underwent 9707 emergency admissions on public holidays. 559 (5·8%) died within 30 days (p<0·0001 vs weekday). 15 routine haematology and biochemistry test results were highly prognostic for mortality. In 271 465 (53·9%) admissions with complete data, adjustment for test results explained 33% (95% CI 21 to 70) of the excess mortality associated with emergency admission on Saturdays compared with Wednesdays, 52% (lower 95% CI 34) on Sundays, and 87% (lower 95% CI 45) on public holidays after adjustment for standard patient characteristics. Excess mortality was predominantly restricted to admissions between 1100 h and 1500 h (pinteraction=0·04). No hospital workload measure was independently associated with mortality (all p values >0·06). INTERPRETATION: Adjustment for routine test results substantially reduced excess mortality associated with emergency admission at weekends and public holidays. Adjustment for patient-level factors not available in our study might further reduce the residual excess mortality, particularly as this clustered around midday at weekends. Hospital workload was not associated with mortality. Together, these findings suggest that the weekend effect arises from patient-level differences at admission rather than reduced hospital staffing or services. FUNDING: NIHR Oxford Biomedical Research Centre.

Same-Day Diagnostic and Surveillance Data for Tuberculosis via Whole-Genome Sequencing of Direct Respiratory Samples.

Routine full characterization of Mycobacterium tuberculosis 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 the point of care. We demonstrate a low-cost method of DNA extraction directly from patient samples for M. tuberculosis WGS. We initially evaluated the method by 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 were obtained from 24 (62%) samples; all results were concordant with reference laboratory phenotypes. Phylogenetic placement was concordant between direct and cultured samples. With Illumina MiSeq/MiniSeq, the workflow from patient sample to results can be completed in 44/16 h at a reagent cost of £96/£198 per sample. We then employed a nonspecific PCR-based library preparation method for sequencing on an Oxford Nanopore Technologies MinION sequencer. We applied this to cultured Mycobacterium bovis strain BCG DNA and to combined culture-negative sputum DNA and BCG DNA. For flow cell version R9.4, the estimated turnaround time from patient to identification of BCG, detection of pyrazinamide resistance, and phylogenetic placement was 7.5 h, with full susceptibility results 5 h later. Antibiotic susceptibility predictions were fully concordant. A critical advantage of 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 DNA in direct samples.

Rapid, comprehensive, and affordable mycobacterial diagnosis with whole-genome sequencing: a prospective study.

BACKGROUND: Slow and cumbersome laboratory diagnostics for Mycobacterium tuberculosis complex (MTBC) risk delayed treatment and poor patient outcomes. Whole-genome sequencing (WGS) could potentially provide a rapid and comprehensive diagnostic solution. In this prospective study, we compare real-time WGS with routine MTBC diagnostic workflows. METHODS: We compared sequencing mycobacteria from all newly positive liquid cultures with routine laboratory diagnostic workflows across eight laboratories in Europe and North America for diagnostic accuracy, processing times, and cost between Sept 6, 2013, and April 14, 2014. We sequenced specimens once using local Illumina MiSeq platforms and processed data centrally using a semi-automated bioinformatics pipeline. We identified species or complex using gene presence or absence, predicted drug susceptibilities from resistance-conferring mutations identified from reference-mapped MTBC genomes, and calculated genetic distance to previously sequenced UK MTBC isolates to detect outbreaks. WGS data processing and analysis was done by staff masked to routine reference laboratory and clinical results. We also did a microcosting analysis to assess the financial viability of WGS-based diagnostics. FINDINGS: Compared with routine results, WGS predicted species with 93% (95% CI 90-96; 322 of 345 specimens; 356 mycobacteria specimens submitted) accuracy and drug susceptibility also with 93% (91-95; 628 of 672 specimens; 168 MTBC specimens identified) accuracy, with one sequencing attempt. WGS linked 15 (16% [95% CI 10-26]) of 91 UK patients to an outbreak. WGS diagnosed a case of multidrug-resistant tuberculosis before routine diagnosis was completed and discovered a new multidrug-resistant tuberculosis cluster. Full WGS diagnostics could be generated in a median of 9 days (IQR 6-10), a median of 21 days (IQR 14-32) faster than final reference laboratory reports were produced (median of 31 days [IQR 21-44]), at a cost of £481 per culture-positive specimen, whereas routine diagnosis costs £518, equating to a WGS-based diagnosis cost that is 7% cheaper annually than are present diagnostic workflows. INTERPRETATION: We have shown that WGS has a scalable, rapid turnaround, and is a financially feasible method for full MTBC diagnostics. Continued improvements to mycobacterial processing, bioinformatics, and analysis will improve the accuracy, speed, and scope of WGS-based diagnosis. FUNDING: National Institute for Health Research, Department of Health, Wellcome Trust, British Colombia Centre for Disease Control Foundation for Population and Public Health, Department of Clinical Microbiology, Trinity College Dublin.

Can rapid integrated polymerase chain reaction-based diagnostics for gastrointestinal pathogens improve routine hospital infection control practice? A diagnostic study.

BACKGROUND: Every year approximately 5000-9000 patients are admitted to a hospital with diarrhoea, which in up to 90% of cases has a non-infectious cause. As a result, single rooms are 'blocked' by patients with non-infectious diarrhoea, while patients with infectious diarrhoea are still in open bays because of a lack of free side rooms. A rapid test for differentiating infectious from non-infectious diarrhoea could be very beneficial for patients. OBJECTIVE: To evaluate MassCode multiplex polymerase chain reaction (PCR) for the simultaneous diagnosis of multiple enteropathogens directly from stool, in terms of sensitivity/specificity to detect four common important enteropathogens: Clostridium difficile, Campylobacter spp., Salmonella spp. and norovirus. DESIGN: A retrospective study of fixed numbers of samples positive for C. difficile (n = 200), Campylobacter spp. (n = 200), Salmonella spp. (n = 100) and norovirus (n = 200) plus samples negative for all these pathogens (n = 300). Samples were sourced from NHS microbiology laboratories in Oxford and Leeds where initial diagnostic testing was performed according to Public Health England methodology. Researchers carrying out MassCode assays were blind to this information. A questionnaire survey, examining current practice for infection control teams and microbiology laboratories managing infectious diarrhoea, was also carried out. SETTING: MassCode assays were carried out at Oxford University Hospitals NHS Trust. Further multiplex assays, carried out using Luminex, were run on the same set of samples at Leeds Teaching Hospitals NHS Trust. The questionnaire was completed by various NHS trusts. MAIN OUTCOME MEASURES: Sensitivity and specificity to detect C. difficile, Campylobacter spp., Salmonella spp., and norovirus. RESULTS: Nucleic acids were extracted from 948 clinical samples using an optimised protocol (200 Campylobacter spp., 199 C. difficile, 60 S. enterica, 199 norovirus and 295 negative samples; some samples contained more than one pathogen). Using the MassCode assay, sensitivities for each organism compared with standard microbiological testing ranged from 43% to 94% and specificities from 95% to 98%, with particularly poor performance for S. enterica. Relatively large numbers of unexpected positives not confirmed with quantitative PCR were also observed, particularly for S. enterica, Giardia lamblia and Cryptosporidium spp. As the results indicated that S. enterica detection might provide generic challenges to other multiplex assays for gastrointestinal pathogens, the Luminex xTag(®) gastrointestinal assay was also run blinded on the same extracts (937/948 remaining) and on re-extracted samples (839/948 with sufficient material). For Campylobacter spp., C. difficile and norovirus, high sensitivities (> 92%) and specificities (> 96%) were observed. For S. enterica, on the original MassCode/Oxford extracts, Luminex sensitivity compared with standard microbiological testing was 84% [95% confidence interval (CI) 73% to 93%], but this dropped to 46% on a fresh extract, very similar to MassCode, with a corresponding increase in specificity from 92% to 99%. Overall agreement on the per-sample diagnosis compared with combined microbiology plus PCR for the main four/all pathogens was 85.6%/64.7%, 87.0%/82.9% and 89.8%/86.8% for the MassCode assay, Luminex assay/MassCode extract and Luminex assay/fresh extract, respectively. Luminex assay results from fresh extracts implied that 5% of samples did not represent infectious diarrhoea, even though enteropathogens were genuinely present. Managing infectious diarrhoea was a significant burden for infection control teams (taking 21% of their time) and better diagnostics were identified as having major potential benefits for patients. CONCLUSIONS: Overall, the Luminex xTag gastrointestinal panel showed similar or superior sensitivity and specificity to the MassCode assay. However, on fresh extracts, this test had low sensitivity to detect a key enteric pathogen, S. enterica; making it an unrealistic option for most microbiology laboratories. Extraction efficiency appears to be a major obstacle for nucleic acid-based tests for this organism, and possibly the whole Enterobacteriaceae family. To improve workflows in service microbiology laboratories, to reduce workload for infection control practitioners, and to improve outcomes for NHS patients, further research on deoxyribonucleic acid-based multiplex gastrointestinal diagnostics is urgently needed. FUNDING: The Health Technology Assessment programme of the National Institute for Health Research.

Assessment of Mycobacterium tuberculosis transmission in Oxfordshire, UK, 2007-12, with whole pathogen genome sequences: an observational study.

BACKGROUND: Patients born outside the UK have contributed to a 20% rise in the UK's tuberculosis incidence since 2000, but their effect on domestic transmission is not known. Here we use whole-genome sequencing to investigate the epidemiology of tuberculosis transmission in an unselected population over 6 years. METHODS: We identified all residents with Oxfordshire postcodes with a Mycobacterium tuberculosis culture or a clinical diagnosis of tuberculosis between Jan 1, 2007, and Dec 31, 2012, using local databases and checking against the national Enhanced Tuberculosis Surveillance database. We used Illumina technology to sequence all available M tuberculosis cultures from identified cases. Sequences were clustered by genetic relatedness and compared retrospectively with contact investigations. The first patient diagnosed in each cluster was defined as the index case, with links to subsequent cases assigned first by use of any epidemiological linkage, then by genetic distance, and then by timing of diagnosis. FINDINGS: Although we identified 384 patients with a diagnosis of tuberculosis, country of birth was known for 380 and we sequenced isolates from 247 of 269 cases with culture-confirmed disease. 39 cases were genomically linked within 13 clusters, implying 26 local transmission events. Only 11 of 26 possible transmissions had been previously identified through contact tracing. Of seven genomically confirmed household clusters, five contained additional genomic links to epidemiologically unidentified non-household members. 255 (67%) patients were born in a country with high tuberculosis incidence, conferring a local incidence of 109 cases per 100,000 population per year in Oxfordshire, compared with 3·5 cases per 100,000 per year for those born in low-incidence countries. However, patients born in the low-incidence countries, predominantly UK, were more likely to have pulmonary disease (adjusted odds ratio 1·8 [95% CI 1·2-2·9]; p=0·009), social risk factors (4·4 [2·0-9·4]; p<0·0001), and be part of a local transmission cluster (4·8 [1·6-14·8]; p=0·006). INTERPRETATION: Although inward migration has contributed to the overall tuberculosis incidence, our findings suggest that most patients born in high-incidence countries reactivate latent infection acquired abroad and are not involved in local onward transmission. Systematic screening of new entrants could further improve tuberculosis control, but it is important that health care remains accessible to all individuals, especially high-risk groups, if tuberculosis control is not to be jeopardised. FUNDING: UK Clinical Research Collaboration (Wellcome Trust, Medical Research Council, National Institute for Health Research [NIHR]), and NIHR Oxford Biomedical Research Centre.

Detection of mixed infection from bacterial whole genome sequence data allows assessment of its role in Clostridium difficile transmission.

Bacterial whole genome sequencing offers the prospect of rapid and high precision investigation of infectious disease outbreaks. Close genetic relationships between microorganisms isolated from different infected cases suggest transmission is a strong possibility, whereas transmission between cases with genetically distinct bacterial isolates can be excluded. However, undetected mixed infections-infection with ≥2 unrelated strains of the same species where only one is sequenced-potentially impairs exclusion of transmission with certainty, and may therefore limit the utility of this technique. We investigated the problem by developing a computationally efficient method for detecting mixed infection without the need for resource-intensive independent sequencing of multiple bacterial colonies. Given the relatively low density of single nucleotide polymorphisms within bacterial sequence data, direct reconstruction of mixed infection haplotypes from current short-read sequence data is not consistently possible. We therefore use a two-step maximum likelihood-based approach, assuming each sample contains up to two infecting strains. We jointly estimate the proportion of the infection arising from the dominant and minor strains, and the sequence divergence between these strains. In cases where mixed infection is confirmed, the dominant and minor haplotypes are then matched to a database of previously sequenced local isolates. We demonstrate the performance of our algorithm with in silico and in vitro mixed infection experiments, and apply it to transmission of an important healthcare-associated pathogen, Clostridium difficile. Using hospital ward movement data in a previously described stochastic transmission model, 15 pairs of cases enriched for likely transmission events associated with mixed infection were selected. Our method identified four previously undetected mixed infections, and a previously undetected transmission event, but no direct transmission between the pairs of cases under investigation. These results demonstrate that mixed infections can be detected without additional sequencing effort, and this will be important in assessing the extent of cryptic transmission in our hospitals.

Transforming clinical microbiology with bacterial genome sequencing.

Whole-genome sequencing of bacteria has recently emerged as a cost-effective and convenient approach for addressing many microbiological questions. Here, we review the current status of clinical microbiology and how it has already begun to be transformed by using next-generation sequencing. We focus on three essential tasks: identifying the species of an isolate, testing its properties, such as resistance to antibiotics and virulence, and monitoring the emergence and spread of bacterial pathogens. We predict that the application of next-generation sequencing will soon be sufficiently fast, accurate and cheap to be used in routine clinical microbiology practice, where it could replace many complex current techniques with a single, more efficient workflow.

Acquisition of Streptococcus pneumoniae and nonspecific morbidity in infants and their families: a cohort study.

BACKGROUND: Most children are believed to acquire Streptococcus pneumoniae asymptomatically, with only a few developing overt S. pneumoniae disease. This study investigates the relationship between acquisition of S. pneumoniae and mild nonspecific infection leading to general practitioner (GP) consultation. METHODS: A prospective birth cohort study of 213 infants assessed at home 9 times during 24 weeks by nasopharyngeal swab and parental interview was conducted. RESULTS: All positive S. pneumoniae swabs (including acquisition and carriage) were significantly associated with GP consultations for infection by the study infant compared with infants with negative swabs [odds ratio (OR), 1.6; 95% confidence interval (CI) 1.1-2.2; P = 0.005]. There was a stronger association with S. pneumoniae acquisition alone (OR 2.1; 95% CI 1.3-3.4; P = 0.001) than with carriage only (OR 1.4; 95% CI 0.9-2.0; P = 0.1). Multivariate analysis confirmed that S. pneumoniae acquisition by the study subject was independently associated with GP consultations: adjusted hazard ratio, 1.8 (95% CI 1.1-2.9); P = 0.02. A similar and independent association was found between S. pneumoniae acquisition by the study subject, and GP consultations for infection by the family (adjusted hazard ratio, 1.8; 95% CI 1.1-2.8; P = 0.01). CONCLUSION: Acquisition of S. pneumoniae by the study infant was significantly associated with GP consultations for infection by the infant or family.