Question: Obtaining input matrix from Kallisto for later downstream analysis with Deseq2 or EdgeR.
gravatar for obuhovadarina
7 months ago by
obuhovadarina0 wrote:

Dear All!

I am very new to bioinformatics and trying to connect phenotype of certain cells with expression of certain genes (all exhibiting the same function). I downloaded FASTA sequences for all genes of interest from NCBI and RNA-seq data from cells. I trimmed adapters from obtained RNA sequences and run Kallisto quant on it. Now I am trying to obtain the matrix from Kallisto to further downstream input.

I am looking at this tutorial However, I am confused with this part

files <- file.path(dir, "kallisto_boot", samples$run, "abundance.h5")

As the authors use tximportData package, I am not sure what files should "kallisto_boot" include. Transcripts? Abundance.tsv and json file?

edger rna-seq kallisto deseq2 • 290 views
ADD COMMENTlink written 7 months ago by obuhovadarina0

I cannot help with kallisto as I do not use it but the index it expects is an entire transcriptome, not just a collection of selected genes. Download a reference transcriptome, either from NCBI/RefSeq or Gencode, then run kallisto, then tximport, normalize data with a tool of your choice, e.g. edgeR or DESeq2 and then do whatever downstream analysis you plan to do. I assume that this kallisto_boot is a column in the output files of the tool that tximport expects. Simply run the example code and see if it works.

ADD REPLYlink written 7 months ago by ATpoint36k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1597 users visited in the last hour