Rna-Seq Analysis With Bacteria, Rrna Reduced, Strand Specific
3
3
Entering edit mode
11.1 years ago
epistatic ▴ 190

I do mostly eukaryotic RNA-seq using Tophat2 and then cufflinks, DEseq, edgeR, etc.

I have some bacterial rRNA- (Ribo-zero) and TruSeq directional (strand-specific) samples, PE-50 on my HiSeq 2500. I am not sure the best way to go about the alignment since these have no introns, do I use bowtie2 or bwa directly or should I still use a GTF and Tophat2?

For GTF files, is there a standard location to retrieve these? If I have a GFF can I use it or do I need to format into a GTF?

Thank you!

rna-seq tophat2 • 4.3k views
ADD COMMENT
2
Entering edit mode
11.1 years ago
epistatic ▴ 190

EDGE-pro, Estimated Degree of Gene Expression in PROkaryots -> DESeq works great for Bacterial RNA-seq. It uses bowtie for the alignment and includes a script to take the FPKM output and create a count table that you can immediately upload into R. I also had a bacterial genome with two plasmids and it handled it all extremely well.

Input is genome fasta, ptt, rnt, and read fastq

http://ccb.jhu.edu/software/EDGE-pro/

ADD COMMENT
0
Entering edit mode
11.1 years ago
Mathew Bunj ▴ 40

Bowtie 2 should be fine, I dont think you need TopHat2, BTW what was the purpose of Ribo zero library?

ADD COMMENT
0
Entering edit mode

The RiboZero library will remove the ribosomal RNA, which typically makes up >95% of the total RNA. Unlike with eukaryotes, we can't enrich for mRNA using polyA tails, so we deplete rRNA instead.

ADD REPLY
0
Entering edit mode
11.1 years ago
epistatic ▴ 190

Response to different anti-microbial agents.

ADD COMMENT

Login before adding your answer.

Traffic: 2235 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