scholarly journals PlasmidSeeker: identification of known plasmids from bacterial whole genome sequencing reads

PeerJ ◽  
2018 ◽  
Vol 6 ◽  
pp. e4588 ◽  
Author(s):  
Märt Roosaare ◽  
Mikk Puustusmaa ◽  
Märt Möls ◽  
Mihkel Vaher ◽  
Maido Remm

BackgroundPlasmids play an important role in the dissemination of antibiotic resistance, making their detection an important task. Using whole genome sequencing (WGS), it is possible to capture both bacterial and plasmid sequence data, but short read lengths make plasmid detection a complex problem.ResultsWe developed a tool named PlasmidSeeker that enables the detection of plasmids from bacterial WGS data without read assembly. The PlasmidSeeker algorithm is based onk-mers and usesk-mer abundance to distinguish between plasmid and bacterial sequences. We tested the performance of PlasmidSeeker on a set of simulated and real bacterial WGS samples, resulting in 100% sensitivity and 99.98% specificity.ConclusionPlasmidSeeker enables quick detection of known plasmids and complements existing tools that assemble plasmids de novo. The PlasmidSeeker source code is stored on GitHub:https://github.com/bioinfo-ut/PlasmidSeeker.

2015 ◽  
Vol 68 (10) ◽  
pp. 835-838 ◽  
Author(s):  
Björn A Espedido ◽  
Borce Dimitrijovski ◽  
Sebastiaan J van Hal ◽  
Slade O Jensen

AimsTo characterise the resistome of a multi-drug resistant Klebsiella pneumoniae (Kp0003) isolated from an Australian traveller who was repatriated to a Sydney Metropolitan Hospital from Myanmar with possible prosthetic aortic valve infective endocarditis.MethodsKp0003 was recovered from a blood culture of the patient and whole genome sequencing was performed. Read mapping and de novo assembly of reads facilitated in silico multi-locus sequence and plasmid replicon typing as well as the characterisation of antibiotic resistance genes and their genetic context. Conjugation experiments were also performed to assess the plasmid (and resistance gene) transferability and the effect on the antibiotic resistance phenotype.ResultsImportantly, and of particular concern, the carbapenem-hydrolysing β-lactamase gene blaNDM-4 was identified on a conjugative IncX3 plasmid (pJEG027). In this respect, the blaNDM-4 genetic context is similar (at least to some extent) to what has previously been identified for blaNDM-1 and blaNDM-4-like variants.ConclusionsThis study highlights the potential role that IncX3 plasmids have played in the emergence and dissemination of blaNDM-4-like variants worldwide and emphasises the importance of resistance gene surveillance.


2017 ◽  
Author(s):  
JR Tyson ◽  
NJ O’Neil ◽  
M Jain ◽  
HE Olsen ◽  
P Hieter ◽  
...  

ABSTRACTAdvances in 3rd generation sequencing have opened new possibilities for ‘benchtop’ whole genome sequencing. The MinION is a portable device that uses nanopore technology and can sequence long DNA molecules. MinION long reads are well suited for sequencing and de novo assembly of complex genomes with large repetitive elements. Long reads also facilitate the identification of complex genomic rearrangements such as those observed in tumor genomes. To assess the feasibility of the de novo assembly of large complex genomes using both MinION and Illumina platforms, we sequenced the genome of a Caenorhabditis elegans strain that contains a complex acetaldehyde-induced rearrangement and a biolistic bombardment-mediated insertion of a GFP containing plasmid. Using ∼5.8 gigabases of MinION sequence data, we were able to assemble a C. elegans genome containing 145 contigs (N50 contig length = 1.22 Mb) that covered >99% of the 100,286,401 bp reference genome. In contrast, using ∼8.04 gigabases of Illumina sequence data, we were able to assemble a C. elegans genome in 38,645 contigs (N50 contig length = ∼26 kb) containing 117 Mb. From the MinION genome assembly we identified the complex structures of both the acetaldehyde-induced mutation and the biolistic-mediated insertion. To date, this is the largest genome to be assembled exclusively from MinION data and is the first demonstration that the long reads of MinION sequencing can be used for whole genome assembly of large (100 Mb) genomes and the elucidation of complex genomic rearrangements.


Author(s):  
M.I. Terekhova ◽  
◽  
E.V. Rogacheva ◽  
I.A. Derevyanchenko ◽  
L.A. Kraeva ◽  
...  

The increasing number of antibiotic-resistant isolates of L. monocytogenes is required to establish a genotypic resistance profile to ensure appropriate antibiotic therapy of listeriosis. In this study, whole-genome sequencing and de novo assembly was performed on L. monocytogenes strains from St. Petersburg and the Vologda region. We obtained the MLST ST, phylogenetic lineage and PCR-serogroups in silico for isolates under the study, revealed genes and mutations associated with antibiotic resistance. In general, the genetic composition was similar between the strains from different regions and included a wide range of antibiotic resistance mechanisms. Listeria strains possessed genes that code for resistance to β-lactam antibiotics, fluoroquinolones, tetracyclines and macrolides, — classes that are commonly used in the treatment of listeria infection. The present study is important in the sanitary and epidemiological surveillance of listeriosis in Russia.


2018 ◽  
Author(s):  
Alba Sanchis-Juan ◽  
Jonathan Stephens ◽  
Courtney E French ◽  
Nicholas Gleadall ◽  
Karyn Mégy ◽  
...  

AbstractComplex structural variants (cxSVs) are genomic rearrangements comprising multiple structural variants, typically involving three or more breakpoint junctions. They contribute to human genomic variation and can cause Mendelian disease, however they are not typically considered during genetic testing. Here, we investigate the role of cxSVs in Mendelian disease using short-read whole genome sequencing (WGS) data from 1,324 individuals with neurodevelopmental or retinal disorders from the NIHR BioResource project. We present four cases of individuals with a cxSV affecting Mendelian disease-associated genes. Three of the cxSVs are pathogenic: a de novo duplication-inversion-inversion-deletion affecting ARID1B in an individual with Coffin-Siris syndrome, a deletion-inversion-duplication affecting HNRNPU in an individual with intellectual disability and seizures, and a homozygous deletion-inversion-deletion affecting CEP78 in an individual with cone-rod dystrophy. Additionally, we identified a de novo duplication-inversion-duplication overlapping CDKL5 in an individual with neonatal hypoxic-ischaemic encephalopathy. Long-read sequencing technology used to resolve the breakpoints demonstrated the presence of both a disrupted and an intact copy of CDKL5 on the same allele; therefore, it was classified as a variant of uncertain significance. Analysis of sequence flanking all breakpoint junctions in all the cxSVs revealed both microhomology and longer repetitive sequences, suggesting both replication and homology based processes. Accurate resolution of cxSVs is essential for clinical interpretation, and here we demonstrate that long-read WGS is a powerful technology by which to achieve this. Our results show cxSVs are an important although rare cause of Mendelian disease, and we therefore recommend their consideration during research and clinical investigations.


2021 ◽  
Author(s):  
Marco Toffoli ◽  
Xiao Chen ◽  
Fritz J Sedlazeck ◽  
Chiao-Yin Lee ◽  
Stephen Mullin ◽  
...  

GBA variants cause the autosomal recessive Gaucher disease, and carriers are at increased risk of Parkinson disease (PD) and Lewy body dementia (LBD). The presence of a highly homologous nearby pseudogene (GBAP1) predisposes to a range of structural variants arising from either gene conversion or reciprocal recombination, the latter resulting in copy number gains or losses, complicating genetic testing and analysis. To date, short-read sequencing has not been able to fully resolve these or other variants in the key homology region, and targeted long-read sequencing has not previously resolved reciprocal recombinants. We present and validate two independent methods to resolve recombinant alleles and other variants in GBA: Gauchian, a novel bioinformatics tool for short-read, whole-genome sequencing data analysis, and Oxford Nanopore long-read sequencing after enrichment with appropriate PCR. The methods were concordant for 42 samples including 30 with a range of recombinants and GBAP1-related mutations, and Gauchian outperforms the GATK Best Practices pipeline. Applying Gauchian to Illumina sequencing of over 10,000 individuals from publicly available cohorts shows that copy number variants (CNVs) spanning GBAP1 are relatively common in Africans. CNV frequencies in PD and LBD are similar to controls, but gains may coexist with other mutations in patients, and a modifying effect cannot be excluded. Gauchian detects a higher frequency of GBA variants in LBD than PD, especially severe ones. These findings highlight the importance of accurate GBA mutation detection in these patients, which is possible by either Gauchian analysis of short-read whole genome sequencing, or targeted long-read sequencing.


2020 ◽  
Vol 58 (4) ◽  
Author(s):  
Ellen N. Kersh ◽  
Cau D. Pham ◽  
John R. Papp ◽  
Robert Myers ◽  
Richard Steece ◽  
...  

ABSTRACT U.S. gonorrhea rates are rising, and antibiotic-resistant Neisseria gonorrhoeae (AR-Ng) is an urgent public health threat. Since implementation of nucleic acid amplification tests for N. gonorrhoeae identification, the capacity for culturing N. gonorrhoeae in the United States has declined, along with the ability to perform culture-based antimicrobial susceptibility testing (AST). Yet AST is critical for detecting and monitoring AR-Ng. In 2016, the CDC established the Antibiotic Resistance Laboratory Network (AR Lab Network) to shore up the national capacity for detecting several resistance threats including N. gonorrhoeae. AR-Ng testing, a subactivity of the CDC’s AR Lab Network, is performed in a tiered network of approximately 35 local laboratories, four regional laboratories (state public health laboratories in Maryland, Tennessee, Texas, and Washington), and the CDC’s national reference laboratory. Local laboratories receive specimens from approximately 60 clinics associated with the Gonococcal Isolate Surveillance Project (GISP), enhanced GISP (eGISP), and the program Strengthening the U.S. Response to Resistant Gonorrhea (SURRG). They isolate and ship up to 20,000 isolates to regional laboratories for culture-based agar dilution AST with seven antibiotics and for whole-genome sequencing of up to 5,000 isolates. The CDC further examines concerning isolates and monitors genetic AR markers. During 2017 and 2018, the network tested 8,214 and 8,628 N. gonorrhoeae isolates, respectively, and the CDC received 531 and 646 concerning isolates and 605 and 3,159 sequences, respectively. In summary, the AR Lab Network supported the laboratory capacity for N. gonorrhoeae AST and associated genetic marker detection, expanding preexisting notification and analysis systems for resistance detection. Continued, robust AST and genomic capacity can help inform national public health monitoring and intervention.


2021 ◽  
Vol 6 (1) ◽  
Author(s):  
Brent S. Pedersen ◽  
Joe M. Brown ◽  
Harriet Dashnow ◽  
Amelia D. Wallace ◽  
Matt Velinder ◽  
...  

AbstractIn studies of families with rare disease, it is common to screen for de novo mutations, as well as recessive or dominant variants that explain the phenotype. However, the filtering strategies and software used to prioritize high-confidence variants vary from study to study. In an effort to establish recommendations for rare disease research, we explore effective guidelines for variant (SNP and INDEL) filtering and report the expected number of candidates for de novo dominant, recessive, and autosomal dominant modes of inheritance. We derived these guidelines using two large family-based cohorts that underwent whole-genome sequencing, as well as two family cohorts with whole-exome sequencing. The filters are applied to common attributes, including genotype-quality, sequencing depth, allele balance, and population allele frequency. The resulting guidelines yield ~10 candidate SNP and INDEL variants per exome, and 18 per genome for recessive and de novo dominant modes of inheritance, with substantially more candidates for autosomal dominant inheritance. For family-based, whole-genome sequencing studies, this number includes an average of three de novo, ten compound heterozygous, one autosomal recessive, four X-linked variants, and roughly 100 candidate variants following autosomal dominant inheritance. The slivar software we developed to establish and rapidly apply these filters to VCF files is available at https://github.com/brentp/slivar under an MIT license, and includes documentation and recommendations for best practices for rare disease analysis.


PLoS ONE ◽  
2021 ◽  
Vol 16 (6) ◽  
pp. e0253440
Author(s):  
Samantha Gunasekera ◽  
Sam Abraham ◽  
Marc Stegger ◽  
Stanley Pang ◽  
Penghao Wang ◽  
...  

Whole-genome sequencing is essential to many facets of infectious disease research. However, technical limitations such as bias in coverage and tagmentation, and difficulties characterising genomic regions with extreme GC content have created significant obstacles in its use. Illumina has claimed that the recently released DNA Prep library preparation kit, formerly known as Nextera Flex, overcomes some of these limitations. This study aimed to assess bias in coverage, tagmentation, GC content, average fragment size distribution, and de novo assembly quality using both the Nextera XT and DNA Prep kits from Illumina. When performing whole-genome sequencing on Escherichia coli and where coverage bias is the main concern, the DNA Prep kit may provide higher quality results; though de novo assembly quality, tagmentation bias and GC content related bias are unlikely to improve. Based on these results, laboratories with existing workflows based on Nextera XT would see minor benefits in transitioning to the DNA Prep kit if they were primarily studying organisms with neutral GC content.


Sign in / Sign up

Export Citation Format

Share Document