VST normalization and presentation of gene expression as VST-Counts_sample,gene_ - mean_gene_ has been suggested as a better alternative to z-scaling https://support.bioconductor.org/p/123969/. We could add this as a parameter, and leave the current method as the default.