small dimensions in the matrix of counts from smallRNA using DESEq2
1
0
Entering edit mode
andreia ▴ 10
@andreia-23745
Last seen 2.6 years ago
Portugal

Hi there,

I am doing smallRNA data analysis and I have a matrix of counts 17 miRNA and 16 samples I wonder if its correct to perform the results() from DESeq2 to get the DEM's.

Is there a minimum number of miRNA that is correct to perform the analysis with DESEq2 package?

Please, can anyone help me about this doubt?

Thanks in advance

Andreia

matrix DESeq2 SmallRNA dimensions small • 1.6k views
ADD COMMENT
1
Entering edit mode

What kind of assay is this? With only 17 miRNAs it is probably not NGS, is it? What are the values, are these counts?

ADD REPLY
0
Entering edit mode

Yes they are raw counts.

I know that is very small data! I did in past 2 data smallRNAseq analysis (not alot experience but i have a steady workflow...) but nothing with this quality. I already told the client about this.

This is from semen samples and the wetlab protocol is hard to get a good quality and quatity of RNA... :/

ADD REPLY
0
Entering edit mode

And yes is NGS :)

ADD REPLY
1
Entering edit mode
@mikelove
Last seen 22 hours ago
United States

With 17 features, I would first doubt the accuracy of the size factors. The dispersion estimates will be decent because of the 16 samples, but if you don't know what is "housekeeping" then there's no hope for DE analysis with so few features.

ADD COMMENT
0
Entering edit mode

What method do you recommend to use in this type of dimensions data smallrna for differential expression?

Thanks

Andreia

ADD REPLY
1
Entering edit mode

I don't have any suggestions here, I'd be worried how to normalize the data.

ADD REPLY
0
Entering edit mode

I tested TMM and quantile normalization like some papers about smallRNA tested, actually i think there is no THE BEST normalization method :/ for this type of the analysis, the same thing with the aligner (i tested bwa and bowtie for this dataset with similar results).

In my particular case its harder, the data that came out its not helping i'm afraid.

ADD REPLY

Login before adding your answer.

Traffic: 963 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6