Hi, I have assembled a new virus genome. So I have contigs and scaffolds fasta file. My question is:
- How can I understand if this new viral genome is circular or linear, starting by contigs or scaffolds file?
- How can I understand if it's single or double stranded?
Are there any tools, software or scripts that do this?
Please, Help me. Thanks Giovanna
2 answers
Do you have paired-end reads? Yes? See how the map. You should have pairs with R1 mapping towards the 5'-end of the assembly with R2 mapping towards the 3'-end (or vice versa). The process is more complicated if you have many contigs (which seems pretty odd for a viral genome assembly). Another tell tale is if the 5'-end of the contig and the 3'-end of the contig (again assuming just one contig assembly) overlap with e.g. self-vs-self blast query..
First map the reads to your assembly. Then filter away all pairs that map near to each other. See what remains. Be sure to allow the reads to map to multiple loci. It's not so easy to answer because you have told us very little about your data. Where did the nucleic acids come from and how were they isolated, sequenced and assembled? Etc. Stuff like this alone could maybe answer questions about ssDNA vs dsDNA..
I had paired fastq file that were mapped with Genieous and I obtained the contigs fasta file. Then I performed SOPRA scaffolding and the statistics of its. Now I would like to know if my genome is single or double stranded and if it is circular or linear. Now I'm testing Circlator, but it is still running. Do you have any other ideas? Thanks all.
I don't think you can easily, if at all, figure out the strandedness of the virus.
As 5hekki points out, in order to decide if it's circular, you can just look to see if the sequences are circularly permuted at each end.
The following tool might be useful:
Log in to answer this question.
-
How to get host and virus UMI count per cell with Bacterial genome data?
written by San •
-
mapping scaffolds(bacteria) to a reference genome
written by noobiestar •
-
Generate single genome consensus fasta file from SPAdes output: viral genome sequencing
written by snow_seq •
-
How to construct draft genome of virus using Trinity de novo assembler?
written by mathavanbioinfo
-
How to assemble viral genomes when my data contains host DNA as well
written by GBC_Zonatos •
-
genome circularity cross validation
written by K.Gee •
-
Tophat problem - R1 and R2 reads mapped to different strandeds of viral genome
written by WZOUXM •
-
grabbing the beginning and end of viral genome in denovo sequensing
written by ebrahimiet •
-
Copy number (repeat count) of assembled contigs
written by misaghb •
-
Circular Genome??
written by Bdv •
The easiest solution will be a wet-lab method.
Thanks, but can you suggest me how?
Off the top of my head, a simple way would be to try and restriction digest your DNA with double-strand specific restriction enzymes. If they work, its probably dsDNA, if not, probably ssDNA.
If you also have to consider RNA viruses things get more complicated, but I'm guessing if that was the case you wouldn't have NGS data for it already.
Spectrophotometry, like nanodrop or qubit, should be able to distinguish ssDNA from dsDNA.