I am writing to ask your help in regard to performing differential analysis to analyze miR RNA-seq data. I have miRNA seq data and I have analyzed in the following way-
Initial QC - removed adaptors at 3’end of fastq files
Used MiRDeep2 Mapper process and mapped reads to HG38 genome
Then used MiRDeep2 for the identification of novel and known miRNAs. In this procedure, I used collapsed reads and mapping output (ARF format) from MiRDeep2 Mapper I used human mature and precursor miRNA sequences from miRBase database.
I have also performed MiRDeep2 Quantifier analysis for quantitation of reads mapping to known miRBase precursors. Now, I would like to perform miRNA differential analysis (control vs treatment) using DESeq2. For that, I have mapped the fastq files to HG38 genome using Bowtie2 and then used featureCounts. But, the DESeq2 tool failed to perform. Could you please help me with this issue? I have also tried using MiRDeep2 Quantifier output to perform DESeq2 but failed again.
Mayur Doke, DESeq2, or edgeR use like input raw counts of transcripts (mRNA, miRNA, and other features). So then, you want to guarantee this condition. But how Asgar Hussain wrote, you need to tell us what's error did you found.
Thank you for your reply Wentao Cai, André Nicolau Aquime Gonçalves and Asgar Hussain . I am sorry for the delayed response. Please find below error that DESeq2 gave me-
Mayur Doke You have some miRNAs with more than one rows. Remove the redundant miRNAs or you can set "check.names =False" while reading expression file.
This error isn't from DESeq2 package. You need to check your expression file to remove duplicated gene symbols. Load your file with rownames = FALSE parameter. So, then you can try to run this:
#Colapse symbols that are duplicated by taking the one with higest expression