Abstract Purpose: Cutaneous T-cell lymphomas (CTCL), encompassing a spectrum of T-cell lymphoproliferative disorders involving the skin, have collectively increased in incidence over the last 40 years. Sézary syndrome is an aggressive form of CTCL characterized by significant presence of malignant cells in both the blood and skin. The guarded prognosis for Sézary syndrome reflects a lack of reliably effective therapy, due, in part, to an incomplete understanding of disease pathogenesis. Experimental Design: Using single-cell sequencing of RNA and the machine-learning reverse graph embedding approach in the Monocle package, we defined a model featuring distinct transcriptomic states within Sézary syndrome. Gene expression used to differentiate the unique transcriptional states were further used to develop a boosted tree classification for early versus late CTCL disease. Results: Our analysis showed the involvement of FOXP3+ malignant T cells during clonal evolution, transitioning from FOXP3+ T cells to GATA3+ or IKZF2+ (HELIOS) tumor cells. Transcriptomic diversities in a clonal tumor can be used to predict disease stage, and we were able to characterize a gene signature that predicts disease stage with close to 80% accuracy. FOXP3 was found to be the most important factor to predict early disease in CTCL, along with another 19 genes used to predict CTCL stage. Conclusions: This work offers insight into the heterogeneity of Sézary syndrome, providing better understanding of the transcriptomic diversities within a clonal tumor. This transcriptional heterogeneity can predict tumor stage and thereby offer guidance for therapy.