Scalable differential expression analysis of single cell transcriptomics datasets with complex study designs