how to calculaute Tissue Specificity of genes
2
0
Entering edit mode
5.3 years ago
Bioiris • 0

salut,

I have a count table , and I need to calculate the tissue specificity of genes

table(A) ————————

gene_id sample1 sample2 sample3 . ……..sample20 g1 23 0 1 3 g2 35 1 4 8 g3 4 6 70 0 g4 98 12 7 6 g5 60 2 8 12 . . . g1907 13 0 0 0

is there a tool that calculates tissue specificity ??

Thanks.

RNA-Seq • 1.7k views
ADD COMMENT
1
Entering edit mode

There are multiple methods to estimate tissue specificity. Here's an article that compares some approaches: http://bib.oxfordjournals.org/content/early/2016/02/17/bib.bbw008.full

ADD REPLY
0
Entering edit mode

I am doing the same thing. I read that paper last night. Very fast read ...recommends the TAU method. But also points to some pitfalls to avoid.

Other papers from this year:

https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4932457/

https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4968975/

https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4916594/

What tissue are you working on?

ADD REPLY
0
Entering edit mode

Try Jaccard index, I have used them with RNAseq data before.

ADD REPLY
0
Entering edit mode
5.3 years ago
Satyajeet Khare ★ 1.6k

You can use Specificity matrix in cummeRbund. S=1 if the gene is expressed exclusively in that sample.

ADD COMMENT
0
Entering edit mode
5.3 years ago
Charles Plessy ★ 2.8k

I like the richness index that is implemented in the "rarefy" function of the "vegan" package for R. For example, if one calculates the richness index on a scale (subsample size) of 10 for a given gene, it is basically asking the question of how many samples would seen be expessing that gene if the total number of counts for that gene were 10. The lower the richness, the higher the tissue specificity. To normalise the expression table before calculating the richness, I downsample all the libraries at the same level. Note as well that a meaningful richness index can not be calculated on a scale that is lower than the total number of counts for a given gene. Thus, in the example above, genes with less than 10 counts in total over all the sample should be excluded.

ADD COMMENT

Login before adding your answer.

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