Immune reconstitution after allogeneic stem cell transplantation is a dynamic and complex process depending on the recipient and donor characteristics, on the modalities of transplantation, and on the occurrence of graft-versus-host disease. Multivariate methods widely used for gene expression profiling can simultaneously analyze the patterns of a great number of biological variables on a heterogeneous set of patients. Here we use these methods on flow cytometry assessment of up to 25 lymphocyte populations to analyze the global pattern of long-term immune reconstitution after transplantation. Immune patterns were most distinct from healthy controls at six months, and had not yet fully recovered as long as two years after transplant. The two principal determinants of variability were linked to the balance of B and CD8+ T cells and of natural killer and B cells, respectively. Recipient’s cytomegalovirus serostatus, cytomegalovirus replication, and chronic graft-versus-host disease were the main factors shaping the immune pattern one year after transplant. We identified a complex signature of under-and over-representation of immune populations dictated by recipient’s cytomegalovirus seropositivity. Finally, we identified dimensions of variance in immune patterns as significant predictors of long-term non-relapse mortality, independently of chronic graft-versus-host disease.