Question: Detecting bacterial TSS from RNA-Seq data, how?
0
gravatar for R.Blues
3.2 years ago by
R.Blues90
Somewhere
R.Blues90 wrote:

Hello everyone,

Sorry for bothering you, but I am not sure if there is a standard protocol for inferring or detecting transcriptional start sites (TSS) from bacterial RNA-Seq paired-end data, and I would be very grateful if someone could help me here. I have been searching and I am not really sure what to do! :)

The data is already mapped vs the reference genome using the gtf annotation file. I just want to determinate the +1 for the different transcripts! Is that as easy as taking the first nucleotide of each operon? Or is there a better way?

I am quite lost, and I would really appreciate any help. Thank you very much.

tss rna-seq bacteria • 2.5k views
ADD COMMENTlink modified 3.0 years ago by Biostar ♦♦ 20 • written 3.2 years ago by R.Blues90
0
gravatar for Asaf
3.2 years ago by
Asaf6.4k
Israel
Asaf6.4k wrote:

There are methods to enrich for primary transcripts by catching the ppp at the 5' end of the molecule (simply digesting 5'p and 5'OH). Take a look at: http://nar.oxfordjournals.org/content/44/W1/W46.long

ADD COMMENTlink written 3.2 years ago by Asaf6.4k

Yes, it's true! Thanks! However, I have to work with this RNA-Seq data, so I cannot redo the experimento. I will take this into account for the next time (thank you!) :)

ADD REPLYlink written 3.2 years ago by R.Blues90
0
gravatar for geek_y
3.2 years ago by
geek_y10.0k
Barcelona
geek_y10.0k wrote:

You can not detect TSS from RNA-Seq data. You need data from protocols like CAGE-Seq.

In RNA-seq, the mRNA is randomly fragmented and sequenced. so the precise TSS can not be determined as the mRNA is sequenced randomly. In protocols like CAGE, they explicitly sequence the 5' end of the transcript, so we know the exact TSS with 1bp resolution and even we can quantify the dominant TSS.

ADD COMMENTlink modified 3.2 years ago • written 3.2 years ago by geek_y10.0k

Hello, this looks to me like a strong statement, could you please substantiate it?

ADD REPLYlink written 3.2 years ago by Anima Mundi2.5k

I updated the answer

ADD REPLYlink written 3.2 years ago by geek_y10.0k

Thanks. I would also say that even if optimized protocols are by definition better at finding TSSs, TSS information can still be extracted by conventional RNA-Seq pipelines.

ADD REPLYlink written 3.2 years ago by Anima Mundi2.5k

Yes. Approximate TSSs can be estimated from RNA-Seq but not precise TSS at 1 bp resolution can not be estimated. It depends on the use case whether we need precise TSS or approximate.

ADD REPLYlink written 3.2 years ago by geek_y10.0k

I have found alternatives as TSSpredator (http://it.inf.uni-tuebingen.de/?page_id=190). I understand that it will not work as well as a CAGE, but I have to work with this RNA-Seq data.

ADD REPLYlink modified 3.2 years ago • written 3.2 years ago by R.Blues90

I am very interested in your experiences with the TSS predator. Is was usefule for the TSS detection with RNA-Seq data? How did you convert the bam file to .gf file format? The .gf file format is a Wiggle Traack Farmat (.wig)? Thank you for your reply!

ADD REPLYlink written 3.1 years ago by h.botond40
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: 1305 users visited in the last hour