This is a test version of Biostars. For the public version, visit https://www.biostars.org.
how to calculaute Tissue Specificity of genes

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

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

2 answers

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

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.

Log in to answer this question.