Flow cytometry is a high-throughput technology used to quan- tify multiple surface and intracellular markers at the level of a sin- gle cell. This enables us to identify cell sub-types, and to determine their relative proportions. Improvements of this technology allow us to describe millions of individual cells from a blood sample using multiple markers. This results in high-dimensional datasets, whose manual analysis is highly time-consuming and poorly reproducible. While several methods have been developed to perform automatic recognition of cell populations, most of them treat and analyze each sample independently. However, in practice, individual samples are rarely independent, especially in longitudinal studies. Here we ana- lyze new longitudinal flow-cytometry data from the DALIA-1 trial which evaluates a therapeutic vaccine against HIV, by proposing a new Bayesian nonparametric approach with Dirichlet process mix- ture (DPM) of multivariate skew t-distributions to perform model based clustering of flow-cytometry data. DPM models directly esti- mate the number of cell populations from the data, avoiding model selection issues, and skew t-distributions provides robustness to out- liers and non-elliptical shape of cell populations. To accommodate repeated measurements, we propose a sequential strategy relying on a parametric approximation of the posterior. We illustrate the good performance of our method on simulated data and on an experimen- tal benchmark dataset. This sequential strategy outperforms all other methods evaluated on the benchmark dataset, and leads to improved performance on the DALIA-1 data.