I am working on a RNASeq data set of Gastric Cancer.
However after running Deseq2, The gene MUC5AC which is one the most important genes in this malignancy is not present among 23367 genes.
I also checked for MUC5AC synonyms, but I could not find it.
I used hg19 gtf for this analysis.
I will appreciate any comments
It is probably filtered-out by the independent filtering when running
results. Try to turn it off with
independentFiltering=FALSE. Also be sure that it is indeed present in the input count matrix but you probably checked that already.
@ATpoint, do you think it is likely that a key gene in some malignancy is lost after independent filtering?
@nazaninhoseinkhan, I would double check if your gene of interest is really present in your genome annotation file. There are major and minor versions, for example,
MUC2, another mucine gene, is present in
GRCh38.93.gtfbut not in
I used hg19 (GRCH37) gtf file for annotation.
I carried out this analysis for a team working on Gastric Cancer.
They told me that they wanted to check the expression of MUC5AC.
All other Mucins are present in deseq2 output but not for MUC5AC.
Thanks for your comment.
I turned off "independentFiltering=FALSE", but nothing changed. The number of genes was again 23367.
Please do this. Be sure the gene is present in the input.
It is not present in the htseq count files