Entering edit mode
3.6 years ago
wathunya.h
•
0
Now, I'm comparing RNA expressions that have RNA-Seq and HTSeq count How can I interpret it together with different unit or Can I convert HTSeq count equivalent RNA-Seq? or if you have other suggestions, please feel free to suggest me.
It's not very clear what you are trying to achieve. Generally, if you are comparing expression values (in an RNA-seq dataset) across different sample groups you want to normalize by using for example median of ratios (default DESeq2 normalization method for DE analysis). RPKM/FPKM and TPM are more suitable for gene to gene comparison within the same sample.
Thanks for the information. Now, I have datasets from different databases such as cBioportal, GEO database. Dataset in cBioportal give mRNA expression as RNA-Seq v2 but some datasets in GEO database give RNA-Seq-htseq count because of analysis with high-throughput methods. For my questions and doubts, I want to compare different datasets in DE analysis. Can I convert RNA count from htseq-count equivalence to same RNA-seq dataset or should not convert and compare direct (head-to-head)?
As far as I understood, I would suggest that in all datasets where you have retrieved raw data you process them until generating a counts matrix with
htseq count. This way, at the end, you have counts matrices generated byhtseq countfor all datasets. However, keep in mind 2 important factors:htseq countdoes not normalize the read counts. Normalization happens in downstream analysis, for example while performing DE analysis in DESeq2 which normalizes the counts by median of ratios as default.Hope I managed to help somehow!
Thanks a lot for the useful suggestion. You are so helpful.
You're welcome! Enjoy your day :)
Already answered here