I was wondering does the verboseScatterplot function in WGCNA allow to specific a few subset of genes? This is a subset overlapping with another analysis of mine and I wanted to highlight those genes in the module mebership vs gene significance plot. Below is the code I am using:

verboseScatterplot(abs(geneModuleMembership[moduleGenes,column]),
                   abs(geneTraitSignificance.chol[moduleGenes,1]),
                   xlab=paste("Module membership in", module,"module"),
                   ylab="Gene significance for total choelsterol",
                   main=paste("Module membership vs gene significance"),
                   cex.main = 1.2,cex.lab = 1.2,cex.axis = 1.2,col="turquoise",
                   abline = TRUE,lmFnc = rlm,ylim=c(-0.1,0.3))


# include your problematic code here with any corresponding output 
# please also include the results of running the following in an R session 

sessionInfo( )