hybridSPAdes: an algorithm for hybrid assembly of short and long reads

A fatal outbreak of ST11 carbapenem-resistant hypervirulent Klebsiella pneumoniae in a Chinese hospital: a molecular epidemiological study

BACKGROUND

Hypervirulent Klebsiella pneumoniae strains often cause life-threatening community-acquired infections in young and healthy hosts, but are usually sensitive to antibiotics. In this study, we investigated a fatal outbreak of ventilator-associated pneumonia caused by a new emerging hypervirulent K pneumoniae strain.

METHODS

The outbreak occurred in the integrated intensive care unit of a new branch of the Second Affiliated Hospital of Zhejiang University (Hangzhou, China). We collected 21 carbapenem-resistant K pneumoniae strains from five patients and characterised these strains for their antimicrobial susceptibility, multilocus sequence types, and genetic relatedness using VITEK-2 compact system, multilocus sequence typing, and whole genome sequencing. We selected one representative isolate from each patient to establish the virulence potential using a human neutrophil assay and Galleria mellonella model and to establish the genetic basis of their hypervirulence phenotype.

FINDINGS

All five patients had undergone surgery for multiple trauma and subsequently received mechanical ventilation. The patients were aged 53-73 years and were admitted to the intensive care unit between late February and April, 2016. They all had severe pneumonia, carbapenem-resistant K pneumoniae infections, and poor responses to antibiotic treatment and died due to severe lung infection, multiorgan failure, or septic shock. All five representative carbapenem-resistant K pneumoniae strains belonged to the ST11 type, which is the most prevalent carbapenem-resistant K pneumoniae type in China, and originated from the same clone. The strains were positive on the string test, had survival of about 80% after 1 h incubation in human neutrophils, and killed 100% of wax moth larvae (G mellonella) inoculated with 1 × 10⁶ colony-forming units of the specimens within 24 h, suggesting that they were hypervirulent K pneumoniae. Genomic analyses showed that the emergence of these ST11 carbapenem-resistant hypervirulent K pneumoniae strains was due to the acquisition of a roughly 170 kbp pLVPK-like virulence plasmid by classic ST11 carbapenem-resistant K pneumoniae strains. We also detected these strains in specimens collected in other regions of China.

INTERPRETATION

The ST11 carbapenem-resistant hypervirulent K pneumoniae strains pose a substantial threat to human health because they are simultaneously hypervirulent, multidrug resistant, and highly transmissible. Control measures should be implemented to prevent further dissemination of such organisms in the hospital setting and the community.

FUNDING

Chinese National Key Basic Research and Development Program and Collaborative Research Fund of Hong Kong Research Grant Council.

Field-based species identification of closely-related plants using real-time nanopore sequencing

Advances in DNA sequencing and informatics have revolutionised biology over the past four decades, but technological limitations have left many applications unexplored. Recently, portable, real-time, nanopore sequencing (RTnS) has become available. This offers opportunities to rapidly collect and analyse genomic data anywhere. However, generation of datasets from large, complex genomes has been constrained to laboratories. The portability and long DNA sequences of RTnS offer great potential for field-based species identification, but the feasibility and accuracy of these technologies for this purpose have not been assessed. Here, we show that a field-based RTnS analysis of closely-related plant species (Arabidopsis spp.) has many advantages over laboratory-based high-throughput sequencing (HTS) methods for species level identification and phylogenomics. Samples were collected and sequenced in a single day by RTnS using a portable, “al fresco” laboratory. Our analyses demonstrate that correctly identifying unknown reads from matches to a reference database with RTnS reads enables rapid and confident species identification. Individually annotated RTnS reads can be used to infer the evolutionary relationships of A. thaliana. Furthermore, hybrid genome assembly with RTnS and HTS reads substantially improved upon a genome assembled from HTS reads alone. Field-based RTnS makes real-time, rapid specimen identification and genome wide analyses possible.

Hybrid assembly with long and short reads improves discovery of gene family expansions

Background

Long-read and short-read sequencing technologies offer competing advantages for eukaryotic genome sequencing projects. Combinations of both may be appropriate for surveys of within-species genomic variation.

Methods

We developed a hybrid assembly pipeline called “Alpaca” that can operate on 20X long-read coverage plus about 50X short-insert and 50X long-insert short-read coverage. To preclude collapse of tandem repeats, Alpaca relies on base-call-corrected long reads for contig formation.

Results

Compared to two other assembly protocols, Alpaca demonstrated the most reference agreement and repeat capture on the rice genome. On three accessions of the model legume Medicago truncatula, Alpaca generated the most agreement to a conspecific reference and predicted tandemly repeated genes absent from the other assemblies.

Conclusion

Our results suggest Alpaca is a useful tool for investigating structural and copy number variation within de novo assemblies of sampled populations.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3927-8) contains supplementary material, which is available to authorized users.

The long reads ahead: de novo genome assembly using the MinION

Nanopore technology provides a novel approach to DNA sequencing that yields long, label-free reads of constant quality. The first commercial implementation of this approach, the MinION, has shown promise in various sequencing applications. This review gives an up-to-date overview of the MinION's utility as a de novo sequencing device. It is argued that the MinION may allow for portable and affordable de novo sequencing of even complex genomes in the near future, despite the currently error-prone nature of its reads. Through continuous updates to the MinION hardware and the development of new assembly pipelines, both sequencing accuracy and assembly quality have already risen rapidly. However, this fast pace of development has also lead to a lack of overview of the expanding landscape of analysis tools, as performance evaluations are outdated quickly. As the MinION is approaching a state of maturity, its user community would benefit from a thorough comparative benchmarking effort of de novo assembly pipelines in the near future. An earlier version of this article can be found on  bioRxiv.

Complete Genome Sequence of the Olive-Infecting Strain Xylella fastidiosa subsp. pauca De Donno

We report here the complete and annotated genome sequence of the plant-pathogenic bacterium Xylella fastidiosa subsp. pauca strain De Donno. This strain was recovered from an olive tree severely affected by olive quick decline syndrome (OQDS), a devastating olive disease associated with X. fastidiosa infections in susceptible olive cultivars.

De novo yeast genome assemblies from MinION, PacBio and MiSeq platforms

Long-read sequencing technologies such as Pacific Biosciences and Oxford Nanopore MinION are capable of producing long sequencing reads with average fragment lengths of over 10,000 base-pairs and maximum lengths reaching 100,000 base- pairs. Compared with short reads, the assemblies obtained from long-read sequencing platforms have much higher contig continuity and genome completeness as long fragments are able to extend paths into problematic or repetitive regions. Many successful assembly applications of the Pacific Biosciences technology have been reported ranging from small bacterial genomes to large plant and animal genomes. Recently, genome assemblies using Oxford Nanopore MinION data have attracted much attention due to the portability and low cost of this novel sequencing instrument. In this paper, we re-sequenced a well characterized genome, the Saccharomyces cerevisiae S288C strain using three different platforms: MinION, PacBio and MiSeq. We present a comprehensive metric comparison of assemblies generated by various pipelines and discuss how the platform associated data characteristics affect the assembly quality. With a given read depth of 31X, the assemblies from both Pacific Biosciences and Oxford Nanopore MinION show excellent continuity and completeness for the 16 nuclear chromosomes, but not for the mitochondrial genome, whose reconstruction still represents a significant challenge.

Characterization of the Polymyxin D Synthetase Biosynthetic Cluster and Product Profile of Paenibacillus polymyxa ATCC 10401

The increasing prevalence of polymyxin-resistant bacteria has stimulated the search for improved polymyxin lipopeptides. Here we describe the sequence and product profile for polymyxin D nonribosomal peptide synthetase from Paenibacillus polymyxa ATCC 10401. The polymyxin D synthase gene cluster comprised five genes that encoded ABC transporters (pmxC and pmxD) and enzymes responsible for the biosynthesis of polymyxin D (pmxA, pmxB, and pmxE). Unlike polymyxins B and E, polymyxin D contains d-Ser at position 3 as opposed to l-α,γ-diaminobutyric acid and has an l-Thr at position 7 rather than l-Leu. Module 3 of pmxE harbored an auxiliary epimerization domain that catalyzes the conversion of l-Ser to the d-form. Structural modeling suggested that the adenylation domains of module 3 in PmxE and modules 6 and 7 in PmxA could bind amino acids with larger side chains than their preferred substrate. Feeding individual amino acids into the culture media not only affected production of polymyxins D₁ and D₂ but also led to the incorporation of different amino acids at positions 3, 6, and 7 of polymyxin D. Interestingly, the unnatural polymyxin analogues did not show antibiotic activity against a panel of Gram-negative clinical isolates, while the natural polymyxins D₁ and D₂ exhibited excellent in vitro antibacterial activity and were efficacious against Klebsiella pneumoniae and Acinetobacter baumannii in a mouse blood infection model. The results demonstrate the excellent antibacterial activity of these unusual d-Ser³ polymxyins and underscore the possibility of incorporating alternate amino acids at positions 3, 6, and 7 of polymyxin D via manipulation of the polymyxin nonribosomal biosynthetic machinery.

Genome sequence of Shigella flexneri strain SP1, a diarrheal isolate that encodes an extended-spectrum β-lactamase (ESBL)

BACKGROUND

Shigellosis is the most common cause of gastrointestinal infections in developing countries. In China, the species most frequently responsible for shigellosis is Shigella flexneri. S. flexneri remains largely unexplored from a genomic standpoint and is still described using a vocabulary based on biochemical and serological properties. Moreover, increasing numbers of ESBL-producing Shigella strains have been isolated from clinical samples. Despite this, only a few cases of ESBL-producing Shigella have been described in China. Therefore, a better understanding of ESBL-producing Shigella from a genomic standpoint is required. In this study, a S. flexneri type 1a isolate SP1 harboring bla_CTX-M-14, which was recovered from the patient with diarrhea, was subjected to whole genome sequencing.

RESULTS

The draft genome assembly of S. flexneri strain SP1 consisted of 4,592,345 bp with a G+C content of 50.46%. RAST analysis revealed the genome contained 4798 coding sequences (CDSs) and 100 RNA-encoding genes. We detected one incomplete prophage and six candidate CRISPR loci in the genome. In vitro antimicrobial susceptibility testing demonstrated that strain SP1 is resistant to ampicillin, amoxicillin/clavulanic acid, cefazolin, ceftriaxone and trimethoprim. In silico analysis detected genes mediating resistance to aminoglycosides, β-lactams, phenicol, tetracycline, sulphonamides, and trimethoprim. The bla _CTX-M-14 gene was located on an IncFII2 plasmid. A series of virulence factors were identified in the genome.

CONCLUSIONS

In this study, we report the whole genome sequence of a bla_CTX-M-14-encoding S. flexneri strain SP1. Dozens of resistance determinants were detected in the genome and may be responsible for the multidrug-resistance of this strain, although further confirmation studies are warranted. Numerous virulence factors identified in the strain suggest that isolate SP1 is potential pathogenic. The availability of the genome sequence and comparative analysis with other S. flexneri strains provides the basis to further address the evolution of drug resistance mechanisms and pathogenicity in S. flexneri.

metaSPAdes: a new versatile metagenomic assembler

While metagenomics has emerged as a technology of choice for analyzing bacterial populations, the assembly of metagenomic data remains challenging, thus stifling biological discoveries. Moreover, recent studies revealed that complex bacterial populations may be composed from dozens of related strains, thus further amplifying the challenge of metagenomic assembly. metaSPAdes addresses various challenges of metagenomic assembly by capitalizing on computational ideas that proved to be useful in assemblies of single cells and highly polymorphic diploid genomes. We benchmark metaSPAdes against other state-of-the-art metagenome assemblers and demonstrate that it results in high-quality assemblies across diverse data sets.

Canu: scalable and accurate long-read assembly via adaptive k-mer weighting and repeat separation

Long-read single-molecule sequencing has revolutionized de novo genome assembly and enabled the automated reconstruction of reference-quality genomes. However, given the relatively high error rates of such technologies, efficient and accurate assembly of large repeats and closely related haplotypes remains challenging. We address these issues with Canu, a successor of Celera Assembler that is specifically designed for noisy single-molecule sequences. Canu introduces support for nanopore sequencing, halves depth-of-coverage requirements, and improves assembly continuity while simultaneously reducing runtime by an order of magnitude on large genomes versus Celera Assembler 8.2. These advances result from new overlapping and assembly algorithms, including an adaptive overlapping strategy based on tf-idf weighted MinHash and a sparse assembly graph construction that avoids collapsing diverged repeats and haplotypes. We demonstrate that Canu can reliably assemble complete microbial genomes and near-complete eukaryotic chromosomes using either Pacific Biosciences (PacBio) or Oxford Nanopore technologies and achieves a contig NG50 of >21 Mbp on both human and Drosophila melanogaster PacBio data sets. For assembly structures that cannot be linearly represented, Canu provides graph-based assembly outputs in graphical fragment assembly (GFA) format for analysis or integration with complementary phasing and scaffolding techniques. The combination of such highly resolved assembly graphs with long-range scaffolding information promises the complete and automated assembly of complex genomes.

The hidden perils of read mapping as a quality assessment tool in genome sequencing

This article provides a comparative analysis of the various methods of genome sequencing focusing on verification of the assembly quality. The results of a comparative assessment of various de novo assembly tools, as well as sequencing technologies, are presented using a recently completed sequence of the genome of Lactobacillus fermentum 3872. In particular, quality of assemblies is assessed by using CLC Genomics Workbench read mapping and Optical mapping developed by OpGen. Over-extension of contigs without prior knowledge of contig location can lead to misassembled contigs, even when commonly used quality indicators such as read mapping suggest that a contig is well assembled. Precautions must also be undertaken when using long read sequencing technology, which may also lead to misassembled contigs.

Assembly of long error-prone reads using de Bruijn graphs

The recent breakthroughs in assembling long error-prone reads were based on the overlap-layout-consensus (OLC) approach and did not utilize the strengths of the alternative de Bruijn graph approach to genome assembly. Moreover, these studies often assume that applications of the de Bruijn graph approach are limited to short and accurate reads and that the OLC approach is the only practical paradigm for assembling long error-prone reads. We show how to generalize de Bruijn graphs for assembling long error-prone reads and describe the ABruijn assembler, which combines the de Bruijn graph and the OLC approaches and results in accurate genome reconstructions.

TruSPAdes: barcode assembly of TruSeq synthetic long reads

The recently introduced TruSeq synthetic long read (TSLR) technology generates long and accurate virtual reads from an assembly of barcoded pools of short reads. The TSLR method provides an attractive alternative to existing sequencing platforms that generate long but inaccurate reads. We describe the truSPAdes algorithm (http://bioinf.spbau.ru/spades) for TSLR assembly and show that it results in a dramatic improvement in the quality of metagenomics assemblies.