Bulk gene expression measurements cloud individual cell differences in heterogeneous cell populations. Here, we describe a protocol for how single-cell gene expression analysis and index sorting by Florescence Activated Cell Sorting (FACS) can be combined to delineate heterogeneity and immunophenotypically characterize molecularly distinct cell populations.