Differential expression between ACE2 + TMPRSS2 + cells and other cells of the matched cell subtype (e.g., Secretory Cells) were performed using the “bimod” likelihood-ratio test within each cell subtype and filtered on Benjamini-Hochberg-corrected p value < 0.05.