RNA-seq analysis between a wildtype vs a double microRNA kockout mouse brain
1
0
Entering edit mode
faiza • 0
@37cf7afe
Last seen 2.0 years ago
Germany

I recently performed RNA-seq analysis on auditory brainstem of mouse.The experiment involves wildtype mouse brains vs a transgenic mouse where two microRNA,s are knocked out. The idea was to see the difference in expression of genes between wildtype vs KO as microRNAs regulate gene expression. Problems: 1.The data is messy where the RPM values between 3 replicates are too different. The PCA plot is all over the place where samples from same area and same mouse don't group together. 2.The differential gene expression by edgeR reveals a total of 5 genes and none of those are interesting to me as future research candidates. The genes I identified from literature as targets of these microRNAs and would be interesting candidates for my research question are in the total list produced by edgeR but they are not significantly different.

My question is if we perform a gene ontology /enrichment analysis ,would it reveal the genes(as differentialy expressed) I am interested in as compared to an over all result that edgeR produces? For example now I have genes related to blood or calcium uptake as differentially regulated but I am interested in genes specifically related to inhibitory synapses and see if their expression is changed in KO samples compared to wildtype ones.

I hope the question is not too confusing since I don't know much about R and these analysis .The analysis was done by the company who did the RNA-seq. My apologies!

```

microRNA gene RNASeqRData edgeR geneontology • 1.1k views
ADD COMMENT
0
Entering edit mode

You aren't use RPM as input into edgeR, right? Doing so might not be the cause of your problem, but it won't be helping. But if replicates are very different, and the effects of your knockout are subtle, your experiment might not turn up anything.

ADD REPLY
0
Entering edit mode

No, the RPM are not used for edgeR..yeah i think the issue is difference between replicates to an extent that we see no significant differences.

ADD REPLY
1
Entering edit mode
@james-w-macdonald-5106
Last seen 1 day ago
United States

The analysis using edgeR is a by-gene analysis that specifically tests for evidence that a given gene is differentially expressed in one group versus another. All of the various gene set tests test for evidence that a set of genes is differentially expressed. A conventional over-representation analysis like a GO hypergeometric test won't be useful in this context, as it needs significant by-gene results.

You might be able to use something like mroast or fry to do self-contained gene set tests and achieve significant results for some interesting sets of genes even if none of the genes are themselves significant, particularly if you choose a more powerful method than the default.

That said, if the incoming data are quite noisy and you don't have sufficient replication, any results you get may be suspect. There is no substitute for good data.

ADD COMMENT
0
Entering edit mode

Ah,Thankyou. I am glad you understood what exactly was my concern. The company suggested that Go test might help us that's why I was curious.We thought of adding more replicates to get a more robust data at the end.I hope that solves the issue.

ADD REPLY

Login before adding your answer.

Traffic: 686 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