Question: how do Differential gene expression analysis by Salmon output
gravatar for lkianmehr
6 days ago by
lkianmehr10 wrote:

Hello, I did transcript quantification by Salmon. now I have an output like this table;

Name    Length  EffectiveLength TPM NumReads                                
ENSMUST00000193812.1    ENSMUSG00000102693.1    OTTMUSG00000049935.1    OTTMUST00000127109.1    RP23-271O17.1-001   RP23-271O17.1   1070    TEC     1070    885.998 12.758544   75.864598
ENSMUST00000082908.1    ENSMUSG00000064842.1    -   -   Gm26206-201 Gm26206 110 snRNA       110 27.847  0   0
ENSMUST00000162897.1    ENSMUSG00000051951.5    OTTMUSG00000026353.2    OTTMUST00000086625.1    AC157543.1-003  Xkr4    4153    processed_transcript        4153    3968.998    0.384028    10.229377
ENSMUST00000159265.1    ENSMUSG00000051951.5    OTTMUSG00000026353.2    OTTMUST00000086624.1    AC157543.1-002  Xkr4    2989    processed_transcript        2989    2804.998    4.503059    84.770623
ENSMUST00000070533.4    ENSMUSG00000051951.5    OTTMUSG00000026353.2    OTTMUST00000065166.1    AC157543.1-001  Xkr4    3634    protein_coding      3634    3449.998    0   0
ENSMUST00000192857.1    ENSMUSG00000102851.1    OTTMUSG00000049958.1    OTTMUST00000127143.1    RP23-317L18.1-001   RP23-317L18.1   480 processed_pseudogene        480 314.826 0   0
ENSMUST00000195335.1    ENSMUSG00000103377.1    OTTMUSG00000049960.1    OTTMUST00000127145.1    RP23-317L18.4-001   RP23-317L18.4   2819    TEC     2819    2634.998    1.92262 34
ENSMUST00000192336.1    ENSMUSG00000104017.1    OTTMUSG00000049961.1    OTTMUST00000127146.1    RP23-317L18.3-001   RP23-317L18.3   2233    TEC     2233    2048.998    0   0

I want to do differential gene expression analysis by DESeq2 on this output and on the basis of Gene_id and NumReads columns, I am not sure actually if aggregated gene name based on gene_id is correct? do you have any suggestion about that? how do DGE analysis on?

thanks in advance

ADD COMMENTlink modified 6 days ago by genomax57k • written 6 days ago by lkianmehr10

Take a look at the tximport package

ADD REPLYlink written 6 days ago by WouterDeCoster32k

Take a look at section 2.1 in RNASeqGene workflow for DESeq2. Detailed information in this tximport vignette.

ADD REPLYlink modified 6 days ago • written 6 days ago by genomax57k

thanks, don't you have any suggestion except tximport, because I am new in R and I can't use it, can I use aggregate option to sum genes ?"

ADD REPLYlink written 6 days ago by lkianmehr10

because I am new in R and I can't use it,

Sure you can, it's well documented

ADD REPLYlink written 6 days ago by WouterDeCoster32k

Please embrace the opportunity to learn a new skill, lkianmehr. In my vocabulary, at least, the phrase "I can't" does not exist.

Make a start from Here.

ADD REPLYlink written 4 days ago by Kevin Blighe30k
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: 1827 users visited in the last hour