This is a test version of Biostars. For the public version, visit https://www.biostars.org.
RNAseq counts: transcript to gene

Hi All,

I have RNAseq counts that I want to use for a differential expression analysis. However, there are approximately 1500 duplicated gene symbols with counts for each transcript (geneID). Can I simply collapse those to create unique rows for each geneSymbol?

For example,

counts <- counts %>%
  select(-c("geneID", "bioType", "annotationLevel")) %>%
  group_by(geneSymbol) %>%
  summarise(across(everything(), sum, na.rm = TRUE))

Thanks for your advice.

enter image description here

rnaseq
• 819 views • link
updated by Ram Ram • written by MAPK2 • MAPK2

1 answer

No. Restrict yourself to canonical chromosomes and you won't run into this issue - as often at least. You cannot collapse counts that map to different loci to the same "gene" just because HGNC and ENSEMBL name things differently. ENSEMBL is more unique so you should ideally pick the entries you want to keep instead of aggregating anything.

• 0 views • link
written by Ram Ram

Log in to answer this question.

More posts like this