Question: I am not getting any plasmid in my mircobial genome data? why is it so ?
0
gravatar for hafiz.talhamalik
9 days ago by
Pakistan
hafiz.talhamalik60 wrote:

I have fastq files of microbial genomics data. and want to get plasmids from it. I have tried plasmidfinder and spadesplasmid as well. but all the time file show up empty. why is it so ? any suggestion ?

ADD COMMENTlink written 9 days ago by hafiz.talhamalik60

How was the DNA prepared?

Do you have a reference sequence for any plasmids you can blast contigs against?

Do you know for certain this organism, under the conditions you tested, contains a plasmid?

ADD REPLYlink written 9 days ago by jrj.healey12k

The DNA was extracted by using whole DNA extraction kit under optimal conditions. Yes, We do have reference sequence of a plasmid.

ADD REPLYlink written 8 days ago by ambrinaakbar10

I think its still possible to lose plasmids (especially if they’re small), when doing preps which purport to be “total DNA preps” (but I’m not convinced this is definitively the answer anyway).

It may be worth taking a step back, and just doing a ‘normal’ de novo assembly, with no special consideration for the plasmids, and then BLAST all the contigs against the reference plasmids sequence. Step one is figuring out whether you even have the corresponding DNA there at all, before plasmidSPAdes etc, to see if they have anything to work with.

ADD REPLYlink written 8 days ago by jrj.healey12k

I tried aligning fastq file with plasmid sequence and found some hits as well but when I annotated those hits with plasidfinder I found nothing there.

ADD REPLYlink written 8 days ago by hafiz.talhamalik60
1

I would assemble before you BLAST - so BLAST contigs instead of reads.

Plasmid assemblers aren't magic however, so assuming you're doing nothing wrong, there's still no guarantee they'll find what you're looking for in a neat way.

ADD REPLYlink written 8 days ago by jrj.healey12k

I have tried plasmidfinder and spadesplasmid as well.

What exactly have you tried? Can you show the commands?

ADD REPLYlink written 9 days ago by h.mon25k

I tried assembly using velvet, spades and unicycler. then I used assembled contigs for plasmid finding for web based plasmidfiner.

for plasmid spades I used this command.

plasmidspades.py -1 fastqfiles/AK-17_S2_L001_R1_001.fastq.gz -2 fastqfiles/AK-17_S2_L001_R2_001.fastq.gz -o spades_output/S2-plasmid
ADD REPLYlink modified 8 days ago by h.mon25k • written 8 days ago by hafiz.talhamalik60
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1634 users visited in the last hour