Biology Stack Exchange is a question and answer site for biology researchers, academics, and students. Join them; it only takes a minute:

Sign up
Here's how it works:
  1. Anybody can ask a question
  2. Anybody can answer
  3. The best answers are voted up and rise to the top

I want to learn more about the data analysis and statistics on transcriptome sequencing data. I would like to read some important papers of the field and books and maybe some MOOCS, if they are available.

More precisely I have data of differentially expressed genes across different groups of individuals and I want to test, if the genes are more expressed in one group are the genes also more polymorphic?

Any ideas?

share|improve this question

closed as too broad by kmm Nov 13 '14 at 13:52

There are either too many possible answers, or good answers would be too long for this format. Please add details to narrow the answer set or to isolate an issue that can be answered in a few paragraphs.If this question can be reworded to fit the rules in the help center, please edit the question.

there are several tools available for this analysis. Tophat-Cufflinks is a popular combo that people use. For analysis, again, there are several tools. – WYSIWYG Mar 19 '14 at 13:54
Are you studying a species that currently have its genome sequenced? – biojl Mar 19 '14 at 16:20
No, it does not have a published genome sequence, de novo transcriptome assemly was used and cuffdiff was used to get the information about differentially expressed genes between two treatments. Now I want to statistically test if the genes, that are more expressed in one treatment also exhibit more polymorphism. Also, the data is very overdispersed. – Piret Avila Mar 20 '14 at 20:58
@PiretAvila. You would need different sample isolates (biological replicates) to study polymorphism. – WYSIWYG Mar 21 '14 at 5:21
Yes, I have many biological replicates for each treatment. So you see, my problem is not obtaining the data, but analyzing it. I have differential gene expression values for genes that were expressed differently in two treatments (log2(FPKMy/FPKMx)) and I know how many SNPs each of these genes have. I want to know if the genes that were more expressed in one treatment (log2(FPKMy/FPKMx) is positive), are they on average more polymorphic. What would be the best statistical method? I used GLM's (hurdle and zero inflated neg. bin.), since there are a lot of sites with no polymorphism. – Piret Avila Mar 21 '14 at 7:33