par(bg = "gray30")
plot(1:25, pch = 1:25, col = hcl.colors(25, palette = "viridis"), bg = "red", lwd = 2, cex = 1.5, axes = FALSE)
axis(2, at = seq(0, 25, by = 5), col = "white", col.axis = "white", las = 1)
axis(1, at = seq(0, 25, by = 5), col = "white", col.axis = "white")
viridis palette. A few symbols, those shown with red fill, accept a background fill and outline color.