Question: How do I specify single end Read in CIRCquant?
0
gravatar for jaqx008
11 months ago by
jaqx00870
jaqx00870 wrote:

Hey everyone, I am trying to identify circular RNAs in my seq data using CIRCquant. I have read and understood the usage and all the files required I have generated. However, my reads are single end reads and the usage and manual for CIRIquant do not provide any option for single end reads input even though in the published paper they mentioned that the program can take in signle end reads. can any one help ? thanks and see usage below for paired end reads

CIRIquant -t 20 --config ./chr1.yml -1 ./fastq1 -2 ./fastq2 --config ./chr1.yml --no-gene -o ./GCF -p GCFloc
ADD COMMENTlink modified 10 months ago by kevinzjy0 • written 11 months ago by jaqx00870

If this is the published paper (https://genomemedicine.biomedcentral.com/articles/10.1186/s13073-019-0614-1), where do they mention that the program can take in single end reads?

Given that the method is based on (reverse) overlap detection of paired-end reads

RO read detection and verification is designed to detect RO reads from paired-end reads based on their 5′ reverse overlaps

feeding it single-end reads does not really make sense to me...

Edit: citation

ADD REPLYlink modified 11 months ago • written 11 months ago by cschu1812.5k

this is the publication and they mentioned it here

https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4316645/#Sec18title under Optimal choice of parameters in CIRI

ADD REPLYlink modified 11 months ago • written 11 months ago by jaqx00870

Hmm, have you tried to only use the -1 option?

ADD REPLYlink written 11 months ago by cschu1812.5k

Yes it doesnt work. at this point I just assume its not possible. Thanks for trying to help.

ADD REPLYlink written 11 months ago by jaqx00870

Just judging from the first couple lines of code at the CIRIquant github repo, the tool will terminate if it doesn't receive -1 and -2, so yea, not without modifying that (which might be easy enough, though!).

ADD REPLYlink written 11 months ago by cschu1812.5k

Yes I agree. I tested with a paired end it ran just fine.

ADD REPLYlink written 11 months ago by jaqx00870
0
gravatar for kevinzjy
10 months ago by
kevinzjy0
Beijing Institutes of Life Science
kevinzjy0 wrote:

Hi jaqx, you can still use CIRI2 for identification of circRNAs from single end data, CIRIquant is designed for quantification of circRNAs based on alignment result of paired-end reads, and does not support single end reads yet.

By the way, our manuscript is still under review, I think what you mentioned is another tool CIRI-full, which is designed to detect internal splicing patterns of circRNAs.

ADD COMMENTlink written 10 months ago by kevinzjy0

Hi, I am sorry for some quick questions about CIRI2:

  1. As I understand, CIRI2 can only be used for RNA-seq data that did not undergo polyA enrichment, isn't it?

  2. Can I use CIRI2 to construct circRNAs from Small RNA-seq data?

Thank you!

ADD REPLYlink modified 7 months ago • written 7 months ago by pbigbig210
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: 805 users visited in the last hour