Multiple comparisions of distribution in single-cell data

I have a single-cell RNAseq data set in which I identified several cell populations. The data set consists of samples derived from different sexes and time points.
What will be the best way to test whether the distribution of a gene expression differs between males and females in each of the identified populations? Do I have to correct for multiple testing? What if in addition do that, I would like to consider the changes in gene expression in two separate days?
I attach the plot demonstrating the data. Please help!