how to build an interaction network in R

I’m investigating the role of my gene of interest (GOI) in cancer, using gene co-expression analysis (guilty by association). Performing gene co-expression analysis, I’ve detected my GOI associated with 163 other genes in a module.

I’ve already performed an enrichment analysis with the genes of this module and I know the potential pathways associated with this gene.

Now I would like to build an interaction network with the genes of this module to better understand how these genes are related and detect the major hubs.

Can you suggest me a package in R to do that? I would like to do it all in R, instead of using softwares as Ingenuity or similar approaches.