Probably you picked up the top differential genes with large fold changes and high expression. I agree that these data are not useless, but one has to remember what the point of differential analysis is. It is not only to pick the top and highly expressed genes which one could probably even pick by examining a MA-plot. It is rather to decide in a data-driven manner what the cutoff is for genes that are variable, moderately- or lowly expressed and have either moderate FCs or high FCs due to small counts, and this is not possible without replicates. The question would be if filtering by a proper logCPM and logFC cutoff would not have resulted in similar findings.
I agree thoigh that my statement that you cannot make any claims was probably too much, as in situations as you decribe it one probably has at least some data-driven basis to make conclusions. It depends what the underlying question is and how dramatic the expression changes are.
I can only shake my head when reading studies like this. First of all they use simulated data, not a single confirmation experiment was done to back up their strategy. Second, they benchmark (among others) against DESeq (deprecated since 2014 upon DESeq2 release) and NOI-seq (now deprecated in the now Bioc version). And third, why would you make a study to develop methods to process inherently underpowered (n=1) studies? If you do not have replicates then you cannot make any claims, simply as that. No way possible you can separate technical noise from biological effects. They should have generated some data themselves (or download experiments with replicates) and then compare their method with standard approaches that use replicates. False discovery rate would probably skyrocket.
Hi- I think I know what you mean and I agree but I think sometimes this statement is exaggerated. There are circumstances where even n=1 can be informative, even if not conclusive, at least to formulate further hypotheses. For example, if you expect low variability to start with (e.g. cell cultures), the effect of the treatment is large and you have some expectation of what genes should change and what should stay the same, even n=1 is something useful. Besides, in some cases n=1 may be all you have so either you do the experiment with that one or you don't do it at all (which may be wiser sometimes but not always, it depends...).
On the other hand, in some cases n=2 or 3 is only marginally better than n=1 and it may give a false sense of security. For example, if you assign people completely blindly to treatment and control, with very small n you have good chances that confounders like sex or age are completely associated to the condition and you get very significant changes which in fact are misleading.
Sorry - it's only that sometimes I have the impression that n=1 is taken as "useless" while n=2 as "ok, now everything's fine".
I agree that replicates are necessary, but I am not sure you cannot make any claims at all. If you have geneA with 10 counts in both conditions and geneB that goes from 10 to 1000 counts, would you say that the probabilities that those genes are differentially expressed are completely equal?
The difference in counts between A and B can be convincing but without randomization of replicates you can't be sure that the difference is due to the treatment applied. The difference may be due to something other than the treatment (e.g. sex, age, handling of samples etc), unless you have some prior belief that tells you that that gene is unlikely to be different for reasons other than the treatment applied (basically I just rephrased my previous comment)
I agree. My main point was that I don't think that the values are completely random in a single-replicate experiment. As you pointed out in the earlier comment, having 2 replicates does not automatically transform complete noise into perfect signal.
Sure, this suggests that this single gene is DEG, but this is nothing you can base a claim on. With claim I mean biological messages like alterations of pathways, groups of genes reacting upon treatment, drug response etc. You are right though that n=1 is not completely useless but for a systematic analysis it is not reliable and I guess most people do a RNA-seq since they want a global picture. For single gene studies one can and should do qPCR or similar low-throughput methods.