Shrinkage estimation of dispersion in Negative Binomial models for RNA-seq experiments with small sample size