How do I specify single end Read in CIRCquant?
1
0
Entering edit mode
2.3 years ago
jaqx008 ▴ 110

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
bwa CircularRNA rnaseq CIRIquant • 1.1k views
ADD COMMENT
0
Entering edit mode

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 REPLY
0
Entering edit mode

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 REPLY
0
Entering edit mode

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

ADD REPLY
0
Entering edit mode

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

ADD REPLY
0
Entering edit mode

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 REPLY
0
Entering edit mode

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

ADD REPLY
0
Entering edit mode
2.2 years ago
kevinzjy • 0

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 COMMENT
0
Entering edit mode

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 REPLY

Login before adding your answer.

Traffic: 1983 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6