Hello,
After using Minimap2 to align my mRNA/(cDNA) transcripts to the GENCODE human genome and counting with featureCounts to the GENCODE gtf file, I have a ridiculously large count of MALAT1 genes for every sample (up to 5% of all reads are this gene, ~5-25K counts/gene non-normalised, sequenced with Nanopore minION).
My sample is bulk rna from resected brain tissue. This gene is apparently up regulated in cancer, so it makes just enough sense to see a lot of counts. But this many counts? I do not know.
Do you all have any experience in dealing with something like this and can give me some tips on how to see if I should exclude this gene or see if this is an error somewhere in my pipeline.
Many Thanks
•
written
by
camppatrick •
1 answer
MALAT1 is a known super-expressed gene. You see tons of it in single cell RNA-seq as well as in bulk.
It's not unusual to get genes like this, depending on the tissue. Sometimes you get lots of mito transcripts (e.g. in muscle), sometimes it's something else. Blood is dominated with globins (there are special kits for globin depletion).
I don't think there's error in your pipeline, you should be able to analyse the dataset in a normal way.
Log in to answer this question.
-
Analysis of Nanopore cDNA sequencing data
written by Nico80 •
-
IGV Sashimi plot from ONT long read sequencing of cDNA amplicons
written by WRUZIC1 •
-
featureCounts cDNA vs DNA files
written by camppatrick •
-
Using normalised counts for pairwise gene correlation (pearson) across samples
written by amjass •
-
RNA-Seq featureCounts Results Annotation File
written by santos48 •
-
Question regarding minimap2 presets
written by Rogerio Ribeiro •
-
Snakemake input error with files from inconsistent naming scheme
written by skbrimer
-
Plenty of gaps with my sam alignment file, why is that?
written by daniel.castanedamogo •
-
question about identifying differential expressed genes in TCGA
written by tujuchuanli
-
Where are my reads in IGV?
written by francois •
