gprofiler_gost_functional_enrichment_gene_list_analysis.py

1from gprofiler import GProfiler
2
3# Initialize the GProfiler object
4gp = GProfiler(return_dataframe=True)
5
6# Define a list of gene identifiers (e.g., Entrez IDs, Ensembl IDs, or gene symbols)
7genes = ["TP53", "TNF", "APOE", "EGFR", "VEGFA"]
8
9# Perform functional enrichment analysis (g:GOSt)
10results = gp.profile(organism='hsapiens', query=genes)
11
12# Print the top enrichment results
13print(results.head())