Differential expression analysis between ACE2 + TMPRSS2 + and negative epithelial cells was performed in Seurat using a Wilcoxon test and Bonferroni p value correction.