how to analyse gene expression using RNA-seq data
2
0
Entering edit mode
7.8 years ago
ashish ▴ 680

I want know how do we calculate gene expression using RNA-seq data available on NCBI. For example I want analyze gene expression of a family of transcription factors, I've the sequences of all these transcription factors and then i download the RNA-seq data from NCBI's SRA database. What do i do next in order to get the FPKM value for creating a heatmap. I am very new to Bioinformatics and learning on my own and also do not know any language except some Biopython and R. Thank you

RNA-Seq gene expression python • 5.0k views
ADD COMMENT
0
Entering edit mode

Check out RNA-seqlopedia, lots of information on understanding the process, not just the technical bioinformatics http://rnaseq.uoregon.edu

ADD REPLY
1
Entering edit mode
7.8 years ago
Vivek ★ 2.7k

You could visit a few RNA seq tutorials. They are by far the easiest to find in bioinformatics. Alternatively you can read through the manuals for the cufflinks suite.

http://cole-trapnell-lab.github.io/cufflinks/manual/

The basic workflow is Reads -> Mapping to genome (STAR/Tophat etc are good splice aware aligners for this job) -> Perform a guided/unguided transcriptome assembly (Cufflinks) -> Create a merged assembly from your samples (Cuffmerge) -> Quantify against a known reference (Cuffcompare)

If you found the sequences from SRA may be you could also find the associated BAM files and reduce your analysis time.

ADD COMMENT
0
Entering edit mode

Thanks Vivek. That is exactly what I was looking for. Can you also explain what is guided/unguided transcriptome assembly?

ADD REPLY
0
Entering edit mode

In guided transcriptome assembly you give the assembler a GTF file of known transcripts/isoforms for the organism you are interested in to aid as a guiding reference. In unguided assembly, you don't which can lead the assembler into assembling more de novo isoforms, if that is your objective.

ADD REPLY
0
Entering edit mode

Alright I got it. That means I'll be doing a guided assembly.

ADD REPLY
2
Entering edit mode
7.8 years ago

Take a look at DESeq2. Here's a very nice tutorial from their authors : https://f1000research.com/articles/4-1070/v2

ADD COMMENT

Login before adding your answer.

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