Differential gene expression analysis based on the negative binomial distribution