This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to separate plasmid from whole complete genome

Hi, actually i had download complete genome from ncbi and few of genome has plasmid, i don't know how to remove it, if anyone know how to separate plasmid from whole complete genome file by using any script please let me know.

genebank ncbi
• 3,133 views • link
updated by GenoMax GenoMax • written by Neel • Neel

Can you provide more information? For example (but not limited to): What genome did you download? What format is it in? What kind of chromosome identifiers are there?

• 0 views • link
written by Trivas • Trivas

Hi, i had download genebank file for example GCF_000006765.1.fna and let's suppose it has plasmid sequence inside complete genome of this file with some plasmid accession number >cp027137.1. so basically i want to separate bacterial genome and it plasmid genome in different file.

Thank you for your time.

• 0 views • link
written by Neel • Neel

The plasmid sequence would be separate from the chromosomal sequence, it is not embedded in the main genome.

Each separate DNA has a separate file.

• 0 views • link
written by Istvan Albert Istvan Albert

Plasmids should generally be a separate entry in genome file. Can you post an example of the accession you are looking at?

• 0 views • link
written by GenoMax GenoMax

Hi, But i had read somewhere that's complete genome have one or multiple chromosomes as well as plasmid and i also check multiple time my genome has plasmid but i don't know how to separate from bacterial genome.

Thank you for your reply

• 0 views • link
written by Neel • Neel

Based on example accession @istvan posted the genome file contains 4 plasmids in addition to the main genome.

$ zgrep ">" GCF_000009725.1_ASM972v1_genomic.fna.gz 
>NC_000911.1 Synechocystis sp. PCC 6803, complete sequence
>NC_005230.1 Synechocystis sp. PCC 6803 plasmid pSYSA, complete sequence
>NC_005231.1 Synechocystis sp. PCC 6803 plasmid pSYSG, complete sequence
>NC_005229.1 Synechocystis sp. PCC 6803 plasmid pSYSM, complete sequence
>NC_005232.1 Synechocystis sp. PCC 6803 plasmid pSYSX, complete sequence

You can use a program like faSplit from Jim Kent to separate the sequences into separate files. Then choose the one that contains the full genome.

• 0 views • link
written by GenoMax GenoMax

Thank you for your reply, i just give this file as a example GCF_000006765.1.

• 0 views • link
written by Neel • Neel

Ok, then I'm wasting my poor neurons here... thanks for the clarification anyway

• 0 views • link
written by patrickdm • patrickdm

patrickdm I used the example that @istvan had provided to Illustrate how plasmids would appear in a genome file that NCBI provides. If the plasmid is integrated in the genome then you will need to use some sort of search strategy to identify it.

• 0 views • link
written by GenoMax GenoMax

GenoMax , yes I eventually got it :) I had just looked up the accession in previous Neel's comment, and got tricked. Thanks.

• 0 views • link
written by patrickdm • patrickdm

1 answer

The plasmids, if known, ought to be annotated for example like so:

Another way to tell plasmids is to find genes known to be on a plasmid or to align against known plasmid databases:

If all you have are reads or fragmented contigs then you may need to use a tool like:

We developed a tool named PlasmidSeeker that enables the detection of plasmids from bacterial WGS data without read assembly. The PlasmidSeeker algorithm is based on k-mers and uses k-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.

• 0 views • link
written by Istvan Albert Istvan Albert

Hi, i am trying to run plasmidseeker and here is the error-

./plasmidseeker.pl -d /home/bvs/neelam/PlasmidSeeker/db_w20 -i /home/bvs/neelam/genome/fna/ -b P_aeruginosa_PaO1.fna -o output2 Loading database... Converting sample reads to k-mers... gt4_wordmap_new: could not mmap file ./16512089695639tmp_sample_32.list gt4_wordmap_new: could not mmap file ./16512089695639tmp_sample_32.list Error: Creating the wordmap failed! Finding coverage of bacterial isolate... gt4_wordmap_new: could not mmap file ./16512089695639tmp_sample_32.list gt4_wordmap_new: could not mmap file ./16512089695639tmp_closest_32.list Error: Creating the wordmap failed! Died at ./plasmidseeker.pl line 156.

Thank you

• 0 views • link
written by Neel • Neel

Log in to answer this question.

More posts like this