I have a list of genes (n = 120) that are involved in breast cancer. And I have a list of differentially expressed genes (n = 70) that are present in the breast cancer gene list. My question is: Are differentially expressed genes significantly more likely to be breast cancer related vs non-differentially expressed genes?
which statistical test will be suitable to address this question and how I can visualize its results in R
.
-
$\begingroup$ I would suggest using an existing tool to evaluate this issue, such as the GO enrichment analysis toolkit: geneontology.org/docs/go-enrichment-analysis. Some of the details are finicky, there are a lot of people who have thought very hard about the issue of annotating genes and GO is what people have settled on, rather than a bespoke list of BC genes. If you must visualize this (rather than simply report the numbers or have a table of enriched categories, which I'd prefer), I'd use a Venn diagram, which you can google for R. $\endgroup$– Maximilian PressJun 14, 2022 at 15:47
1 Answer
This question can be addressed by gene set enrichment analysis. You should use Gene Ontology and KEGG Pathways to find out which genes from your gene set are enriched in GO terms and KEGG pathway terms. You can use clusterProfiler package in R for this anlysis and for visualization also. Here is the link to the guide: https://yulab-smu.top/biomedical-knowledge-mining-book/enrichment-overview.html