Short-read Whole-genome Sequencing Research Articles

On the (im)possibility of reconstructing plasmids from whole-genome short-read sequencing data.

To benchmark algorithms for automated plasmid sequence reconstruction from short-read sequencing data, we selected 42 publicly available complete bacterial genome sequences spanning 12 genera, containing 148 plasmids. We predicted plasmids from short-read data with four programs (PlasmidSPAdes, Recycler, cBar and PlasmidFinder) and compared the outcome to the reference sequences. PlasmidSPAdes reconstructs plasmids based on coverage differences in the assembly graph. It reconstructed most of the reference plasmids (recall=0.82), but approximately a quarter of the predicted plasmid contigs were false positives (precision=0.75). PlasmidSPAdes merged 84 % of the predictions from genomes with multiple plasmids into a single bin. Recycler searches the assembly graph for sub-graphs corresponding to circular sequences and correctly predicted small plasmids, but failed with long plasmids (recall=0.12, precision=0.30). cBar, which applies pentamer frequency analysis to detect plasmid-derived contigs, showed a recall and precision of 0.76 and 0.62, respectively. However, cBar categorizes contigs as plasmid-derived and does not bin the different plasmids. PlasmidFinder, which searches for replicons, had the highest precision (1.0), but was restricted by the contents of its database and the contig length obtained from de novo assembly (recall=0.36). PlasmidSPAdes and Recycler detected putative small plasmids (<10 kbp), which were also predicted as plasmids by cBar, but were absent in the original assembly. This study shows that it is possible to automatically predict small plasmids. Prediction of large plasmids (>50 kbp) containing repeated sequences remains challenging and limits the high-throughput analysis of plasmids from short-read whole-genome sequencing data.

Covert dissemination of carbapenemase-producing Klebsiella pneumoniae (KPC) in a successfully controlled outbreak: long- and short-read whole-genome sequencing demonstrate multiple genetic modes of transmission.

BackgroundCarbapenemase-producing Enterobacteriaceae (CPE), including KPC-producing Klebsiella pneumoniae (KPC-Kpn), are an increasing threat to patient safety.ObjectivesTo use WGS to investigate the extent and complexity of carbapenemase gene dissemination in a controlled KPC outbreak.Materials and methodsEnterobacteriaceae with reduced ertapenem susceptibility recovered from rectal screening swabs/clinical samples, during a 3 month KPC outbreak (2013–14), were investigated for carbapenemase production, antimicrobial susceptibility, variable-number-tandem-repeat profile and WGS [short-read (Illumina), long-read (MinION)]. Short-read sequences were used for MLST and plasmid/Tn4401 fingerprinting, and long-read sequence assemblies for plasmid identification. Phylogenetic analysis used IQTree followed by ClonalFrameML, and outbreak transmission dynamics were inferred using SCOTTI.ResultsTwenty patients harboured KPC-positive isolates (6 infected, 14 colonized), and 23 distinct KPC-producing Enterobacteriaceae were identified. Four distinct KPC plasmids were characterized but of 20 KPC-Kpn (from six STs), 17 isolates shared a single pKpQIL-D2 KPC plasmid. All isolates had an identical transposon (Tn4401a), except one KPC-Kpn (ST661) with a single nucleotide variant. A sporadic case of KPC-Kpn (ST491) with Tn4401a-carrying pKpQIL-D2 plasmid was identified 10 months before the outbreak. This plasmid was later seen in two other species and other KPC-Kpn (ST14,ST661) including clonal spread of KPC-Kpn (ST661) from a symptomatic case to nine ward contacts.ConclusionsWGS of outbreak KPC isolates demonstrated blaKPC dissemination via horizontal transposition (Tn4401a), plasmid spread (pKpQIL-D2) and clonal spread (K. pneumoniae ST661). Despite rapid outbreak control, considerable dissemination of blaKPC still occurred among K. pneumoniae and other Enterobacteriaceae, emphasizing its high transmission potential and the need for enhanced control efforts.

Corrigendum: In silico serotyping of E. coli from short read data identifies limited novel O-loci but extensive diversity of O:H serotype combinations within and between pathogenic lineages.

[This corrects the article DOI: 10.1099/mgen.0.000064.].

Genomic epidemiology of global Klebsiella pneumoniae carbapenemase (KPC)-producing Escherichia coli

The dissemination of carbapenem resistance in Escherichia coli has major implications for the management of common infections. blaKPC, encoding a transmissible carbapenemase (KPC), has historically largely been associated with Klebsiella pneumoniae, a predominant plasmid (pKpQIL), and a specific transposable element (Tn4401, ~10 kb). Here we characterize the genetic features of blaKPC emergence in global E. coli, 2008–2013, using both long- and short-read whole-genome sequencing. Amongst 43/45 successfully sequenced blaKPC-E. coli strains, we identified substantial strain diversity (n = 21 sequence types, 18% of annotated genes in the core genome); substantial plasmid diversity (≥9 replicon types); and substantial blaKPC-associated, mobile genetic element (MGE) diversity (50% not within complete Tn4401 elements). We also found evidence of inter-species, regional and international plasmid spread. In several cases blaKPC was found on high copy number, small Col-like plasmids, previously associated with horizontal transmission of resistance genes in the absence of antimicrobial selection pressures. E. coli is a common human pathogen, but also a commensal in multiple environmental and animal reservoirs, and easily transmissible. The association of blaKPC with a range of MGEs previously linked to the successful spread of widely endemic resistance mechanisms (e.g. blaTEM, blaCTX-M) suggests that it may become similarly prevalent.

Whole-Genome Sequence of the Metastatic PC3 and LNCaP Human Prostate Cancer Cell Lines.

The bone metastasis-derived PC3 and the lymph node metastasis-derived LNCaP prostate cancer cell lines are widely studied, having been described in thousands of publications over the last four decades. Here, we report short-read whole-genome sequencing (WGS) and de novo assembly of PC3 (ATCC CRL-1435) and LNCaP (clone FGC; ATCC CRL-1740) at ∼70 × coverage. A known homozygous mutation in TP53 and homozygous loss of PTEN were robustly identified in the PC3 cell line, whereas the LNCaP cell line exhibited a larger number of putative inactivating somatic point and indel mutations (and in particular a loss of stop codon events). This study also provides preliminary evidence that loss of one or both copies of the tumor suppressor Capicua (CIC) contributes to primary tumor relapse and metastatic progression, potentially offering a treatment target for castration-resistant prostate cancer (CRPC). Our work provides a resource for genetic, genomic, and biological studies employing two commonly-used prostate cancer cell lines.

Global Molecular Epidemiology of IMP-Producing Enterobacteriaceae.

International data on the molecular epidemiology of Enterobacteriaceae with IMP carbapenemases are lacking. We performed short-read (Illumina) whole-genome sequencing on a global collection of 38 IMP-producing clinical Enterobacteriaceae (2008 to 2014). IMP-producing Enterobacteriaceae (7 varieties within 11 class 1 integrons) were mainly present in the South Pacific and Asia. Specific blaIMP-containing integrons (In809 with blaIMP-4, In722 with blaIMP-6, and In687 with blaIMP-14) were circulating among different bacteria in countries such as Australia, Japan, and Thailand. In1312 with blaIMP-1 was present in Klebsiella pneumoniae from Japan and Citrobacter freundii from Brazil. Klebsiella pneumoniae (n = 22) was the most common species; clonal complex 14 (CC14) from Philippines and Japan was the most common clone and contained In1310 with blaIMP-26 and In1321 with blaIMP-6 The Enterobacter cloacae complex (n = 9) consisted of Enterobacter hormaechei and E. cloacae cluster III. CC78 (from Taiwan) containing In73 with blaIMP-8 was the most common clone among the E. cloacae complex. This study highlights the importance of surveillance programs using the latest molecular techniques for providing insight into the characteristics and global distribution of Enterobacteriaceae with blaIMP genes.

Short-Read Whole-Genome Sequencing for Laboratory-Based Surveillance of Bordetella pertussis

Bordetella pertussis is a Gram-negative bacterium that causes respiratory infections in humans. Ongoing molecular surveillance of B. pertussis acellular vaccine (aP) antigens is critical for understanding the interaction between evolutionary pressures, disease pathogenesis, and vaccine effectiveness. Methods currently used to characterize aP components are relatively labor-intensive and low throughput. To address this challenge, we sought to derive aP antigen genotypes from minimally processed short-read whole-genome sequencing data generated from 40 clinical B. pertussis isolates and analyzed using the SRST2 bioinformatic package. SRST2 was able to identify aP antigen genotypes for all antigens with the exception of pertactin, possibly due to low read coverage in GC-rich low-complexity regions of variation. Two main genotypes were observed in addition to a singular third genotype that contained an 84-bp deletion that was identified by SRST2 despite the issues in allele calling. This method has the potential to generate large pools of B. pertussis molecular data that can be linked to clinical and epidemiological information to facilitate research of vaccine effectiveness and disease severity in the context of emerging vaccine antigen-deficient strains.

Clinical detection of deletion structural variants in whole-genome sequences.

Optimal management of acutely ill infants with monogenetic diseases requires rapid identification of causative haplotypes. Whole-genome sequencing (WGS) has been shown to identify pathogenic nucleotide variants in such infants. Deletion structural variants (DSVs, >50 nt) are implicated in many genetic diseases, and tools have been designed to identify DSVs using short-read WGS. Optimisation and integration of these tools into a WGS pipeline could improve diagnostic sensitivity and specificity of WGS. In addition, it may improve turnaround time when compared with current CNV assays, enhancing utility in acute settings. Here we describe DSV detection methods for use in WGS for rapid diagnosis in acutely ill infants: SKALD (Screening Konsensus and Annotation of Large Deletions) combines calls from two tools (Breakdancer and GenomeStrip) with calibrated filters and clinical interpretation rules. In four WGS runs, the average analytic precision (positive predictive value) of SKALD was 78%, and recall (sensitivity) was 27%, when compared with validated reference DSV calls. When retrospectively applied to a cohort of 36 families with acutely ill infants SKALD identified causative DSVs in two. The first was heterozygous deletion of exons 1–3 of MMP21 in trans with a heterozygous frame-shift deletion in two siblings with transposition of the great arteries and heterotaxy. In a newborn female with dysmorphic features, ventricular septal defect and persistent pulmonary hypertension, SKALD identified the breakpoints of a heterozygous, de novo 1p36.32p36.13 deletion. In summary, consensus DSV calling, implemented in an 8-h computational pipeline with parameterised filtering, has the potential to increase the diagnostic yield of WGS in acutely ill neonates and discover novel disease genes.

In silico serotyping of E. coli from short read data identifies limited novel O-loci but extensive diversity of O:H serotype combinations within and between pathogenic lineages.

The lipopolysaccharide (O) and flagellar (H) surface antigens of Escherichia coli are targets for serotyping that have traditionally been used to identify pathogenic lineages. These surface antigens are important for the survival of E. coli within mammalian hosts. However, traditional serotyping has several limitations, and public health reference laboratories are increasingly moving towards whole genome sequencing (WGS) to characterize bacterial isolates. Here we present a method to rapidly and accurately serotype E. coli isolates from raw, short read WGS data. Our approach bypasses the need for de novo genome assembly by directly screening WGS reads against a curated database of alleles linked to known and novel E. coli O-groups and H-types (the EcOH database) using the software package srst2. We validated the approach by comparing in silico results for 197 enteropathogenic E. coli isolates with those obtained by serological phenotyping in an independent laboratory. We then demonstrated the utility of our method to characterize isolates in public health and clinical settings, and to explore the genetic diversity of >1500 E. coli genomes from multiple sources. Importantly, we showed that transfer of O- and H-antigen loci between E. coli chromosomal backbones is common, with little evidence of constraints by host or pathotype, suggesting that E. coli ‘strain space’ may be virtually unlimited, even within specific pathotypes. Our findings show that serotyping is most useful when used in combination with strain genotyping to characterize microevolution events within an inferred population structure.