After a search in Gene Expression Omnibus (GEO), series with accession number GSE2457 comprising of 5 normal and 5 diabetic-ED rats were selected. 
Gene Expression Omnibus (GEO)
Investigation of gene expression and genetic simultaneous control associated with erectile dysfunction and diabetes.Diabetes can cause some diseases or abnormalities.  ED is sexual dysfunction characterized by the inability to establish or maintain an erect penis during sexual activity and is a common problem of men with chronic type 2 diabetes.  After a search in Gene Expression Omnibus (GEO), series with accession number GSE2457 comprising of 5 normal and 5 diabetic-ED rats were selected.  Raw CEL files of these samples were normalized with robust multi-array average (RMA) expression measure method by using the linear models for microarray data (LIMMA) R package.  Then, differentially expressed genes (DEGs) were identified between control and normal penile mucosa by employing the LIMMA R package. 
RMA
Then, differentially expressed genes (DEGs) were identified between control and normal penile mucosa by employing the LIMMA R package. 
LIMMA R package
DEGs were classified by utilizing KEGG to underlying pathways by Enrichr. 
KEGG
DEGs were classified by utilizing KEGG to underlying pathways by Enrichr. 
Enrichr
The expression values of DEGs were used to construct a gene regulatory network (GRN), by the GENEI3 R package. 
GENEI3 R package
Genes with higher betweenness centrality scores were then identified, through the CytoNCA.  We then took the commonality of DEGs genes and high-top ranking genes from CytoNCA via a predicted interaction network using GeneMANIA as the most likely important genes in erectile dysfunction. 
CytoNCA
We then took the commonality of DEGs genes and high-top ranking genes from CytoNCA via a predicted interaction network using GeneMANIA as the most likely important genes in erectile dysfunction. 
GeneMANIA