Performing enrichment analysis using extracted GOs from topGO

Question

[KennedyOmufwoko]

Hello everyone, does anyone have an idea on performing vigorous enrichment analysis in non model species? I have managed to extract the GO using topGO and I have a list my 'top genes of interest' by adj p value and fold change. From my understanding it would be much easier to use packages like clusterprofiler, GSEA or msigdb but my study species hasn't been explored that much. Here's the data

Hoping to do something of this sort Clustering of DAVID gene enrichment results from gene expression studies