Search

OBJECTIVES: In this study, we characterize a concurrent disseminated infection with a virulent hypermucoviscous (HMV) Klebsiella pneumoniae and an OXA-181-producing XDR K. pneumoniae from a patient with recent hospitalization in India. During exposure to meropenem therapy, the highly susceptible HMV K. pneumoniae became resistant to carbapenems, consistent with the acquisition of bla(OXA-181). METHODS: Twelve K. pneumoniae isolates were recovered from the patient and the hospital room environment over a 3 month hospitalization. Phenotypic and molecular studies were completed to characterize the isolates. Oxford Nanopore and Illumina MiSeq WGS were performed to study phylogeny (MLST and SNPs), plasmids and virulence genes and demonstrate changes in the organism’s resistome that occurred over time. RESULTS: WGS revealed that the HMV K. pneumoniae belonged to ST23 and harboured an IncH1B virulence plasmid, while the XDR K. pneumoniae belonged to ST147 and possessed two MDR plasmids (IncR and IncFII), the bla(OXA-181)-bearing ColKP3 plasmid and chromosomal mutations conferring the XDR phenotype. Sequential isolates demonstrated plasmid diversification (fusion of the IncR and IncFII plasmids), mobilization of resistance elements (ompK35 inactivation by ISEcp1-bla(CTX-M-15) mobilization, varying numbers of resistance genes on plasmid scaffolds) and chromosomal mutations (mutations in mgrB) leading to further antibiotic resistance that coincided with antibiotic pressure. Importantly, the HMV strain in this study was unable to preserve the carbapenem-resistant phenotype without the selective pressure of meropenem. CONCLUSIONS: To the best of our knowledge, we are the first to report a carbapenem-resistant HMV K. pneumoniae strain in the USA. Ultimately, this case demonstrates the role of antibiotic pressure in the acquisition and loss of important genetic elements.

SummaryGenomeScope is an open-source web tool to rapidly estimate the overall characteristics of a genome, including genome size, heterozygosity rate, and repeat content from unprocessed short reads. These features are essential for studying genome evolution, and help to choose parameters for downstream analysis. We demonstrate its accuracy on 324 simulated and 16 real datasets with a wide range in genome sizes, heterozygosity levels, and error rates.Availability and Implementationhttp://genomescope.org,https://github.com/schatzlab/genomescope.gitContactmschatz@jhu.edu.Supplementary informationSupplementary data are available atBioinformaticsonline.

Third-generation long-range DNA sequencing and mapping technologies are creating a renaissance in high-quality genome sequencing. Unlike second-generation sequencing, which produces short reads a few hundred base-pairs long, third-generation single-molecule technologies generate over 10,000 bp reads or map over 100,000 bp molecules. We analyze how increased read lengths can be used to address longstanding problems in de novo genome assembly, structural variation analysis and haplotype phasing.

The free-living flatworm, Macrostomum lignano has an impressive regenerative capacity. Following injury, it can regenerate almost an entirely new organism because of the presence of an abundant somatic stem cell population, the neoblasts. This set of unique properties makes many flatworms attractive organisms for studying the evolution of pathways involved in tissue self-renewal, cell-fate specification, and regeneration. The use of these organisms as models, however, is hampered by the lack of a well-assembled and annotated genome sequences, fundamental to modern genetic and molecular studies. Here we report the genomic sequence of M. lignano and an accompanying characterization of its transcriptome. The genome structure of M. lignano is remarkably complex, with ∼75% of its sequence being comprised of simple repeats and transposon sequences. This has made high-quality assembly from Illumina reads alone impossible (N50 = 222 bp). We therefore generated 130× coverage by long sequencing reads from the Pacific Biosciences platform to create a substantially improved assembly with an N50 of 64 Kbp. We complemented the reference genome with an assembled and annotated transcriptome, and used both of these datasets in combination to probe gene-expression patterns during regeneration, examining pathways important to stem cell function.

Background Datasets generated on deep-sequencing platforms have been deposited in various public repositories such as the Gene Expression Omnibus (GEO), Sequence Read Archive (SRA) hosted by the NCBI, or the DNA Data Bank of Japan (ddbj). Despite being rich data sources, they have not been used much due to the difficulty in locating and analyzing datasets of interest. Results Geoseq http://geoseq.mssm.edu provides a new method of analyzing short reads from deep sequencing experiments. Instead of mapping the reads to reference genomes or sequences, Geoseq maps a reference sequence against the sequencing data. It is web-based, and holds pre-computed data from public libraries. The analysis reduces the input sequence to tiles and measures the coverage of each tile in a sequence library through the use of suffix arrays. The user can upload custom target sequences or use gene/miRNA names for the search and get back results as plots and spreadsheet files. Geoseq organizes the public sequencing data using a controlled vocabulary, allowing identification of relevant libraries by organism, tissue and type of experiment. Conclusions Analysis of small sets of sequences against deep-sequencing datasets, as well as identification of public datasets of interest, is simplified by Geoseq. We applied Geoseq to, a) identify differential isoform expression in mRNA-seq datasets, b) identify miRNAs (microRNAs) in libraries, and identify mature and star sequences in miRNAS and c) to identify potentially mis-annotated miRNAs. The ease of using Geoseq for these analyses suggests its utility and uniqueness as an analysis tool.

Background The use of high throughput genome-sequencing technologies has uncovered a large extent of structural variation in eukaryotic genomes that makes important contributions to genomic diversity and phenotypic variation. When the genomes of different strains of a given organism are compared, whole genome resequencing data are typically aligned to an established reference sequence. However, when the reference differs in significant structural ways from the individuals under study, the analysis is often incomplete or inaccurate. Results Here, we use rice as a model to demonstrate how improvements in sequencing and assembly technology allow rapid and inexpensive de novo assembly of next generation sequence data into high-quality assemblies that can be directly compared using whole genome alignment to provide an unbiased assessment. Using this approach, we are able to accurately assess the ‘pan-genome’ of three divergent rice varieties and document several megabases of each genome absent in the other two. Conclusions Many of the genome-specific loci are annotated to contain genes, reflecting the potential for new biological properties that would be missed by standard reference-mapping approaches. We further provide a detailed analysis of several loci associated with agriculturally important traits, including the S5 hybrid sterility locus, the Sub1 submergence tolerance locus, the LRK gene cluster associated with improved yield, and the Pup1 cluster associated with phosphorus deficiency, illustrating the utility of our approach for biological discovery. All of the data and software are openly available to support further breeding and functional studies of rice and other species. Electronic supplementary material The online version of this article (doi:10.1186/s13059-014-0506-z) contains supplementary material, which is available to authorized users.