I’m facing an issue while running edgeR on one of my datasets (endometrial samples). The same pipeline and settings worked correctly for my colorectal dataset, but for this one I’m consistently getting flat results (0 upregulated and 0 downregulated genes).
I’ve rechecked the input count matrix and rerun the analysis multiple times, but the issue persists. Could you please help me check if there might be a problem with how the data is being interpreted or processed on Galaxy?
I don’t see an obvious problem with the technical steps we can see, but to explore this scientifically, you could compare to a pipeline that performs the mapping and counting as a distinct steps! This may help to tune the parameters more with RNA-Star. Below are some Galaxy resources we can suggest for the exploration!
This workflow from the IWC workflow library is one choice. It runs QA through counting. These counts can be input into any of the Bioconductor differential expression tools!
Then, as an example, this workflow consumes the output of the first and generates statistics with DESeq2.
If you didn’t perform QA on the reads first, that is one area to back up to try to improve the mapping/feature capture results with your current pipeline. Removing untrimmed adaptor may be more important for some samples?
Or, you can try a template workflow like this one (or, see the fastp step included in the first counting workflow above!).
I would also suggest confirming the strand for these reads! This exploration can inform about upstream data preparation steps, too. The workflows from this intro tutorial can run some numbers for you.
And this step in another covers reviewing mapping in IGV.
Then this step in another explains more about how to review results from tools like Infer Experiment.
Overall, a workflow would likely be a benefit to ensure you are running the exact same steps and parameters between your two groups of samples. This could be extracted from one of the histories, then tuned up to ensure it has the correct Inputs and connections, then run on both as a batch?
