low counts from microRNA seq
1
0
Entering edit mode
@abhaykanodia-23398
Last seen 23 days ago
United States

I have 24 samples on which smallRNA seq was performed, most of the genes have 0 counts o0r very low counts. Is this normal for microRNA seq? Attached a screenshot of .csv file with count matrix.

I have three different treatment groups. Due to this, all the treatment groups have same differentially expressed microRNAs but the log2fold change is different. Thank you!

Code should be placed in three backticks as shown below


# include your problematic code here with any corresponding output
# please also include the results of running the following in an R session

sessionInfo( )

DESeq2 microRNA HTSeq • 289 views
1
Entering edit mode
ATpoint ★ 2.1k
@atpoint-13662
Last seen 1 hour ago
Germany

The support site is intended to help with technical problems towards Bioconductor packages. This is not a DESeq2 problem. Better communities would be biostars.org or a suitable sub-reddit. In any case smallRNA-seq sequences small RNA, your table shows some protein-coding genes, these would naturally not be captured here. Check existing pipelines how to process these sorts of data, for example https://nf-co.re/smrnaseq and others to be found via your favourit search engine.

0
Entering edit mode

Hi @ATpoint,

Thanks for your comments. However, when I perform DESeq2, I get the following results even after filtering for miRNA, not sure if its a DESeq2 problem or not. Only 45 miRNAs have some log2foldchange value in it, the rest do not.

0
Entering edit mode

Basemean of zero means these genes have no counts, it's something upstream of DESeq2.

0
Entering edit mode

Thanks so much @atpoint. So that means all the miroRNAs should have had some count value? Its unusual to have a countmatrix with such counts for microRNA sequencing?

0
Entering edit mode

0
Entering edit mode

Thanks for your response. I have posted it on biostars.