Question: How can Investigate if a Viral genome is circular or linear? And if it is single or double stranded?
0
gravatar for giovannaventola3es9
24 months ago by
giovannaventola3es930 wrote:

Hi, I have assembled a new virus genome. So I have contigs and scaffolds fasta file. My question is:

  1. How can I understand if this new viral genome is circular or linear, starting by contigs or scaffolds file?
  2. 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

assembly genome • 874 views
ADD COMMENTlink modified 24 months ago by Joe18k • written 24 months ago by giovannaventola3es930

How can I understand if it's single or double stranded?

The easiest solution will be a wet-lab method.

ADD REPLYlink written 24 months ago by WouterDeCoster44k

Thanks, but can you suggest me how?

ADD REPLYlink written 24 months ago by giovannaventola3es930
1

Spectrophotometry, like nanodrop or qubit, should be able to distinguish ssDNA from dsDNA.

ADD REPLYlink written 24 months ago by WouterDeCoster44k

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.

ADD REPLYlink written 24 months ago by Joe18k
3
gravatar for 5heikki
24 months ago by
5heikki9.1k
Finland
5heikki9.1k wrote:

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..

ADD COMMENTlink modified 24 months ago • written 24 months ago by 5heikki9.1k

Hi, yes I have paired reads and I have many contigs. So how can I do?

ADD REPLYlink written 24 months ago by giovannaventola3es930

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..

ADD REPLYlink modified 24 months ago • written 24 months ago by 5heikki9.1k

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.

ADD REPLYlink written 24 months ago by giovannaventola3es930
2
gravatar for Joe
24 months ago by
Joe18k
United Kingdom
Joe18k wrote:

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:

https://github.com/sanger-pathogens/circlator

ADD COMMENTlink written 24 months ago by Joe18k
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: 1706 users visited in the last hour