output file conversion from FeatureCount to DESeq2 - no duplicate for samples
1
0
Entering edit mode
Jiayi • 0
@70976543
Last seen 12 days ago
Hong Kong

Hello, I did a bulkRNA-seq and now have an output gene count file from: featureCounts -s 0 -p -P -d 0 -D 1000 -B --primary -t exon -g gene_name -a gtf -T 6 -o output bam1 bam2 bam3 (I did it via hisat2 then samtools sort then featurecounts using linux command line)

The three bam files belong to 3 cell lines and I want to do a differential analysis on their RNA gene expression, see which cell line expresses higher level of what genes.

The problem is, I did not do any duplicates, so I only have one sample per each cell line, and when I tried doing dds1 <- DESeq(dds1) it tells me:

The design matrix has the same number of samples and coefficients to fit,
so estimation of dispersion is not possible. Treating samples
as replicates was deprecated in v1.20 and no longer supported since v1.22.


What should I do if I want to compare them and get a result on which gene is expressed higher in one cell line compared to another?

Meanwhile, my data looks like below:

            H1_LIM9 NiPSC_LIM9 RUESC_LIM9
DDX11L1     0         0            0
WASH7P      217       209          116
MIR6859-1   1         0            0
MIR1302-2HG 0         0            2


Currently I'm importing them into dds via splitting the data matrix into 3 files each with two cell lines, so that they get to be compared with only one other cell line. Is there other better ways?

Thank you!

DESeq2 FeatureCount • 124 views
0
Entering edit mode
0
Entering edit mode
Assa Yeroslaviz ★ 1.5k
@assa-yeroslaviz-1597
Last seen 9 days ago
Germany

see here:

In the last one, Michael Love also mentions, that there's no real support for DESeq w.o. replicates.