I have a set of differentially expressed genes, and want to know what are their GO terms, and pathway that they reside in. I have installed GO.db and KEGG.db, but unsure how to get started. For examples, genes <- c("TP53", "SOX4", "PTEN"), whats next ?