Cutaneous T-cell lymphoma (CTCL) is a rare cancer of skin-homing T cells. A subgroup of patients develops large cell transformation with rapid progression to an aggressive lymphoma. Here, we investigated the transformed CTCL (tCTCL) tumor ecosystem using integrative multiomics spanning whole-exome sequencing (WES), single-cell RNA sequencing, and immune profiling in a unique cohort of 56 patients. WES of 70 skin biopsies showed high tumor mutation burden, UV signatures that are prognostic for survival, exome-based driver events, and most recurrently mutated pathways in tCTCL. Single-cell profiling of 16 tCTCL skin biopsies identified a core oncogenic program with metabolic reprogramming toward oxidative phosphorylation (OXPHOS), cellular plasticity, upregulation of MYC and E2F activities, and downregulation of MHC I suggestive of immune escape. Pharmacologic perturbation using OXPHOS and MYC inhibitors demonstrated potent antitumor activities, whereas immune profiling provided in situ evidence of intercellular communications between malignant T cells expressing macrophage migration inhibitory factor and macrophages and B cells expressing CD74.
Our study contributes a key resource to the community with the largest collection of tCTCL biopsies that are difficult to obtain. The multiomics data herein provide the first comprehensive compendium of genomic alterations in tCTCL and identify potential prognostic signatures and novel therapeutic targets for an incurable T-cell lymphoma.
This article is highlighted in the In This Issue feature, p. 1171
INTRODUCTION
Cutaneous T-cell lymphomas (CTCL) are a group of non-Hodgkin T-cell lymphomas. They affect males more often than females, and Black/African American (AA) patients show a more aggressive clinical course and inferior survival (1, 2). The most common form of CTCL, mycosis fungoides (MF), is characterized by cutaneous patch/plaque (PP) lesions or tumors, whereas Sezary syndrome (SS) is the leukemic variant of CTCL with circulating malignant T cells in the peripheral blood (3–6). Large cell transformation (LCT) occurs in a subset of patients with MF and SS when the lymphoma cells undergo histopathologic transition from neoplastic small- to-medium-sized lymphocytes to large, blast-like T cells. The primary site of detectable MF/SS transformation is the skin (7, 8). Similar to Richter transformation in chronic lymphocytic leukemia, transformation in CTCL heralds immediate transition to aggressive clinical behavior, especially for those who transform within two years of MF diagnosis (8), rapid decline in survival (median survival 19–36 months; refs. 7–10) and resistance to multiple forms of therapy.
Major advances have been made in the treatment of advanced-stage CTCL with the FDA approval of brentuximab vedotin (anti-CD30) in 2017 (ALCANZA trial; ref. 11) and mogamulizumab (anti-CCR4) as a breakthrough therapy in 2018 (MAVORIC trial; ref. 12). Although brentuximab shows superior efficacy in the skin compartment and in CD30+ MF with LCT (13) and mogamulizumab shows more reliable activity in the blood compartment (14, 15), the responses are not universal or durable. The median progression-free survival was 15.9 months for brentuximab and 7.7 months for mogamulizumab, and notably, LCT was one of the exclusion criteria in the MAVORIC trial. These trials thus highlight a critical unmet need to identify novel therapeutic targets for advanced-stage CTCL, particularly following disease transformation.
Although cancer death rates have significantly declined for many common cancers in the past decade, there is a sobering underrepresentation of this success in rare cancers, and particularly in vulnerable racial and ethnic minority groups (16, 17). The lack of research tissue samples and clinical trials further compounds this inequality. Several research groups have attempted to profile the genomic landscape of SS/MF by whole-exome sequencing (WES) and whole-genome sequencing (18–29). Although we have gained a wealth of information from these multi-institutional sequencing efforts, there is a lack of genome-level investigation of CTCL at disease transformation, with only seven transformed MF WES reported to date (20, 23). Similarly, recent advances in single-cell technologies have provided a high-resolution window into malignant and benign T-cell transcriptomics in a small number of SS and MF samples and revealed evidence of intertumoral and intratumoral heterogeneity (30–32). Malignant T cells in the skin showed upregulation of T-cell activation, T-cell receptor (TCR) ligation, and cell-cycle progression transcripts compared with those in the blood compartment (32). Despite these efforts, the paucity of investigation of CTCL tumor and immune microenvironment (TIME) at LCT or its relationship to precursor PP lesions contributes to inadequate knowledge concerning potential therapeutic targets for this deadly disease state.
To address this critical unmet need, we tackled the challenging transformed CTCL (tCTCL) TIME by applying integrative multiomics and multiplex immune profiling of skin biopsies from a rare cohort of 56 patients with tCTCL. We comprehensively characterized the genomic landscape of tCTCL using tissue resource of similar size to other rare cancers from The Cancer Genome Atlas (TCGA) and established tCTCL as a high tumor mutation burden (TMB) cancer dominated by UV signatures that are prognostic for enhanced survival. Notably, Black/AA patients in our cohort show significantly lower contribution of UV signatures compared with the white patients. We identified predicted driver genes and recurrently altered pathways in Hippo, RAS/RTK, and Notch and showed that tCTCL in skin exhibits a distinct genomic chromosomal copy-number variation (CNV) profile from SS/leukemic CTCL, which has important therapeutic implications. Using a combination of single-cell RNA sequencing (scRNA-seq), scV(D)J-seq, and CNV inference analytic strategies, we identified a unique malignant T-cell program with enrichment for oxidative phosphorylation (OXPHOS), cellular plasticity, upregulation of MYC and E2F activities, and downregulation of MHC I at transformation suggestive of immune escape. Furthermore, in vitro pharmacologic studies using novel small-molecule inhibitors of OXPHOS (IACS-010759; ref. 33) and MYC (MYCi975; ref. 34) demonstrated potent antitumor activity. Immune profiling further revealed receptor–ligand interaction between macrophage migration inhibitory factor (MIF) in malignant T cells and CD74 in antigen-presenting cells (APC). Notably, malignant T cells in tCTCL demonstrated intratumoral genetic and transcriptional heterogeneity, with upregulation of ribosomal protein subunit gene expression in dominant subclones in patients with the poorest clinical outcomes. Collectively, our study provided the first comprehensive compendium of genomic alterations in CTCL at disease transformation, identified a tCTCL oncogenic program that exploits metabolic reprogramming, cellular plasticity, and proliferation, and highlighted potential therapeutic vulnerabilities for this incurable rare cancer.
RESULTS
tCTCL Is a High TMB Cancer and Presence of UV Signatures Is Prognostic for Survival at LCT
We collected a unique cohort of 56 patients with biopsy-proven tCTCL (53 transformed MF, 3 SS/overlap MF-SS with transformed tumors in skin) who were seen and/or treated between 2014 and 2020. Of the 56 patients, there were 38 males (68%), 18 females (32%), and 17 patients belonging to racial/ethnic minority groups (30%). Clinical restaging was performed at the time of LCT. The median time from initial MF diagnosis to LCT was 25.3 months (range, 0–252.5 months), median time from LCT to death or last follow-up was 20.4 months (range, 1.0–61.3 months), and 26 patients were deceased (46.4% of the cohort; Supplementary Table S1). Available biopsies from these patients were carefully annotated with clinical data elements and processed for multiomics profiling, including WES, 5′ scRNA-seq, scV(D)J-seq, and multiplex immunofluorescence (mIF) immune profiling (Fig. 1A; Supplementary Table S2A–S2D).
We first sought to explore the genomic landscape of tCTCL and collected 70 skin biopsies from 54 patients with confirmed LCT in skin for WES [Fig. 1A; 45 transformed tumors (TT), 9 PP with LCT, and 16 concurrent or precursor PPs to TT]. Of the 54 patients sequenced, 41 had matched germline normal samples and 17 had paired PP and TT (Supplementary Table S2A; Methods). Somatic single-nucleotide variants (SSNV) and small indels were called by comparing tumor and PP to matched normal germline or panel of normal (n = 13 patients). We identified a median of 290 mutations per sample with a dominance of missense mutations (97.1%) and C>T transitions (65.8%) that is comparable to previously described in SS and MF (Supplementary Fig. S1A and S1B; refs. 19, 35). We measured the TMB in all samples with matched germline and compared this with 33 other cancers profiled in TCGA and two independent cohorts of SS (Wang and Choi cohorts; refs. 26, 36). Analyses of mutation burdens of SSNVs and small indels in tCTCL demonstrated a TMB that is exceeded only by lung squamous cell carcinoma and cutaneous melanoma and significantly higher than SS (Fig. 1B; tumor types ordered by median TMB; Supplementary Table S3A and S3B).
To comprehensively interrogate the operative mutational processes in tCTCL, we next cataloged the repertoire of mutational signatures with reference to the COSMIC mutational signatures (v3.2; ref. 37). At the sample level, we used deconstructSigs (38) to analyze the weights of mutation signatures and showed that UV signatures SBS7a and SBS7b carried the maximum weight (Fig. 1C; Supplementary Table S4). Similarly, we implemented MutationalPatterns (39) and demonstrated highest cosine similarity to signatures SBS7a and SBS7b, with no significant difference between PPs and tumors (Supplementary Fig. S2A). At the cohort level, novel mutational signatures derived by nonnegative matrix factorization (NMF) methods showed best match to SBS7b and SBS6 (Supplementary Fig. S2B). The predominance of UV signatures has been reported in SS/MF and other sun-protected cancers (22, 35, 40, 41), and Signature 7 has been shown to be prognostic for favorable outcome in cutaneous melanoma (42). We performed survival analysis and showed that the sum of weighted signatures SBS7a, SBS7b, SBS7c, and SBS7d are also prognostic for favorable survival from the time of LCT (Fig. 1D; Supplementary Table S4). Importantly, Black/AA patients in our cohort show significantly lower contribution of UV signatures compared with the white patients (Fig. 1E), with enrichment of SBS6, SBS1, and SBS87 in patients with low UV signatures (Fig. 1F, Supplementary Fig. S3A). Previous large-scale epidemiology studies showed that Black/AA patients with MF demonstrated poorer survival than white patients (1, 2). We also observed a survival gap between Black/AA and white patients in this cohort (Supplementary Fig. S3B).
Exome-Based Sequencing Identifies Driver Events and Key Oncogenic Pathways in tCTCL
We next sought to detect exome-based driver events in tCTCL utilizing two mutation-based computational tools, dNdScv (43) and MutSigCV (44), to enhance the coverage of our analysis. dNdScv is a statistical model based on refined dN/dS while factoring variation of mutation rate, and MutSigCV assesses mutation significance as a function of gene size, trinucleotide context, and background mutation frequency for highly recurrent mutations. When both dNdScv and MutsigCV were performed on the same data set, a total of 20 potential driver genes were identified (Supplementary Table S5A–S5E), including CCR4 (chemokine receptor), FRG1 (chromatin modifying), TEKT4 (cell motility), CDC27 and ESX1 (cell cycle), MTRNR2L2 (antiapoptosis), and other novel genes previously not implicated in CTCL (Fig. 1G; Supplementary Fig. S4A and S4B). Notably, we observed CCR4 mutations concentrated at the C-terminus, similar to other gain-of-function CCR4 C-terminus mutations seen in adult T-cell leukemia/lymphoma (ATLL) and SS (26, 45). We subsequently performed analyses of recurrently mutated pathways and showed top-ranking mutations in the Hippo (FAT1 and FAT4), Notch (NCOR, SPEN, and CREBBP), and RAS (ROS1, MEK, and MAPK3) pathways and p53 (Fig. 1H; Supplementary Table S6A–S6C). Importantly, Hippo and Notch are key signaling pathways known to be essential for organ size control in development, stemness, and tumor suppression, and deletion of FAT1 has been shown to promote stemness and metastasis through epithelial-to-mesenchymal transition (EMT) in cancer (46). We also stratified the WES samples by CD30 expression (≤1%, >1% to <10%, 10%–50%, >50%; Supplementary Methods; Supplementary Table S7) but did not detect significant differences in overall survival, TMB, UV signatures, driver or top-ranking pathway genes (Supplementary Fig. S5A–S5E).
Cutaneous tCTCL Exhibits Distinct Genomic Gains and Losses from Those of SS/Leukemic CTCL
Previous work in SS reported distinct chromosomal CNV patterns, with characteristic recurrent alteration in genes such as ARID1A, CDKN2A, CDKN1B, ZEB1, DNMT3A, PLCG1, TP53, PDCD1, and CARD11 (21, 22). We next explored the tCTCL genome for candidate somatic copy-number alterations (SCNA). Using GISTIC 2.0, an algorithm that detects genes targeted by SCNAs that drive cancer growth (47), we mapped 42 recurrently deleted loci (33 with Q-values <0.005) and 28 recurrently amplified loci (9 with Q-value <0.005) across the genome (Q-threshold 0.25; Fig. 2A; Supplementary Table S8A–S8C). Recurrent deletions were seen in the PRAME family of genes (cancer testis antigen/stemness), clustered Protocadherin family, Defensins (DEFB109P1 and DEFB130), RECQL4 (genomic stability and chromatin modifying), CDKN2A (cell cycle), SOC6 and SOC7 (cytokine signaling, JAK–STAT), and POLE (DNA repair), whereas recurrent amplifications were seen in TWIST1 (cellular plasticity/stemness), NDUFB2 and NDUFA7 (OXPHOS), RPLs and RPSs (ribosomal protein subunits), LGALS3 (Galectin3/RAS), AHR (T-cell differentiation), CARD11 (T-cell activation), and SYCP1 (cancer testis antigen). We also performed GISTIC analysis on available SS/leukemic CTCL samples (ref. 36; Fig. 2A; Supplementary Table S8D–S8F). Although tCTCL and SS/leukemic CTCL shared conserved alterations in genes such as CARD11 and CDKN2A, the global genomic profiles between the two entities are strikingly different, with more prominent G-scores seen in tCTCL (Fig. 2A).
Recent WES meta-analysis of predominantly patients with SS described 55 putative CTCL driver genes involving pathways that affect chromatin remodeling, immune surveillance, MAPK, NFκB, and PI3K signaling (19). We therefore investigated these candidate genes in tCTCL and compared with a publicly available SS/leukemic CTCL data set from this meta-analysis study (Choi cohort; ref. 36). Interestingly, although both cohorts show similarly recurrent mutations in genes such as STAT5B, POT1, TET2, and CARD11, tCTCL demonstrated drastically more recurrent mutations in the JAK–STAT pathway genes (JAK1, JAK3, and STAT3), key chromatin-modifying genes (KMT2C, KMT2D, ARID1A, and CREBBP), TP53, and CDKN2A (Fig. 2B), suggesting more global dysregulation in genome-wide topology and gene expression in tCTCL (Supplementary Table S9A–S9D). Recurrent mutations previously described as hotspot mutations in T-cell and/or natural killer (NK) T-cell lymphomas such as CD28F51l, JAK3A573V, and STAT5BN642H were conserved in tCTCL, with the latter showing potential therapeutic vulnerability to JAK–STAT inhibitors (48, 49). An unbiased comparison of the tCTCL and SS cohorts identified 116 significantly differentially mutated genes, including novel genes such as FSIP2, TEKT4, and Protocadherin 15, which are mutated only in tCTCL but not in SS (P < 0.05, Q < 0.25; Supplementary Table S9E). Dissimilarity in the genomic landscape of transformed T cells in skin versus leukemic T cells in blood provides opportunities to exploit differential or synergistic therapeutic vulnerabilities in the two body compartments at advanced-stage disease.
Dissecting the Transformed CTCL TIME at Single-Cell Resolution
To interrogate the tCTCL tumor ecosystem, we next sought to profile the tCTCL TIME at single-cell resolution. We collected 16 fresh skin biopsies from 8 patients with tCTCL, each with paired PP and TT, and profiled 34,912 cells using complementary 5′ scRNA-seq and scV(D)J-seq strategies (10× chromium; Fig. 1A, Fig. 3A). Single-cell V(D)J-seq identifies the precise TCR α- and β-chain clonotype combination in individual T cells and allows linkage of clonotype information to the whole transcriptome by complementary 5′ scRNA-seq. Despite the advantage of scV(D)J-seq in identifying clonal TCR, it is still imprecise due to “dropout” read, a well-known phenomenon in both scRNA-seq and scV(D)J-seq data sets. To tackle this issue, we sought to separate malignant from benign T cells by two complementary profiling strategies. First, we defined true malignant and true benign T cells by scV(D)J-seq (Fig. 3A; Supplementary Table S10). Next, we distinguished malignant cells by patterns of inferred SCNAs using the 5′ scRNA-seq data (inferCNV; Fig. 3B; Methods). In 7 of 8 patients, we observed a dominant TCR clonotype of paired α- and β-chain sequences, consistent with TCR monoclonality in cutaneous tCTCL (Supplementary Table S10). Concurrent PP and TT lesions from the same patient also shared the same TCR clonotype. Interestingly, even though some lesions appear to have multiple “dominant” clonotypes (e.g., PT35T, clonotypes 1–7), careful inspection reveals the same TCR α-chain (TRA: CAVGDAGGTSYGKLTF) and β-chain sequences (TRB: CASSRQGGSGNTIYF or CASSRLGGSGNTIYF) across all seven clonotypes. Two TRB amino acid sequences that differed by single amino acid were observed in the most dominant clonotype (TRA: CAVGDAGGTSYGKLTF; TRB: CASSRQGGSGNTIYF; TRB:CASSRLGGSGNTIYF), perhaps due to somatic hypermutation. Cells exhibiting unpaired TRA and TRB (e.g., PT35 clonotypes 4, 6, and 7) were present, compatible with “dropout” reads in either TCR α- or β-chain expression, respectively. Due to the observed TCR read dropout, a second malignancy detection method using inferCNV was implemented.
An important basis for inferCNV is the requirement for reference “normal” cells, ideally of the same cell type, and a separate group of “observation” cells for testing. Here, we selected cells harboring the nondominant/polyclonal TCR clonotypes from scV(D)J profiling as “true normal.” We further split the “true normal” cells and input two thirds of these cells into the inferCNV reference cell group (Fig. 3A and B). For the observation group, we “spiked in” the remaining “true normal” cells and input all cells with dominant TCR clonotype (i.e., “true malignant”), as well as all CD3+ T cells with TCR dropout as “malignant suspect.” We hypothesized that malignant T cells would demonstrate malignant CNV patterns and form distinct clusters from normal T cells with neutral CNV patterns. Indeed, we observed a striking dichotomy in the chromosomal large-scale CNV pattern between malignant versus normal cells, with increased patient-specific CNVs observed in malignant T cells (Fig. 3B). Uniform manifold approximation and projection (UMAP) clustering of all cells generated 17 clusters, with distinct separation between malignant T cells harboring positive CNV patterns and benign T cells with neutral CNV patterns (Fig. 3C). As in other cancer types, tCTCL demonstrates significant interpatient heterogeneity by gene expression, with malignant T cells clustering by patient of origin, and normal T cells from different patients clustering together, suggesting consistency of the background immune infiltrate.
Notably, of the three samples that lacked dominant TCR clonotypes (PT53 TT, PT55 PP, an dPT55 TT), PT53 TT showed CNV-neutral patterns in all cells, consistent with lack of viable tumor cells in that sample, whereas PT55 PP and PT55 TT showed cell populations with striking malignant CNV patterns, consistent with the presence of scV(D)J-seq dropout. Therefore, our approach demonstrates that the combination of scV(D)J-seq and CNV evaluation provides a more robust methodology for distinguishing malignant from benign reactive T cells in T-cell lymphoma single-cell studies. After separating malignant T cells from the tCTCL TIME, the remaining CD45/PTPRC+ benign immune cell populations were annotated by known marker genes (Fig. 3C; Methods; Supplementary Fig. S6; Supplementary Table S11A–S11C) and SingleR, a computational algorithm that performs unbiased cell type recognition by referencing transcriptomic data sets of pure cell types (Monaco immune database; ref. 50; Methods). Using these two approaches, we identified the benign immune cell clusters as B cells, NK cells, macrophages/monocytes, dendritic cells, and various T-cell subtypes. For the nonimmune cells in the TIME, fibroblasts were identified by COL1A2 and endothelial cells by VWF expression.
tCTCL Exploits OXPHOS Metabolic Reprogramming, Cellular Plasticity, and MYC/E2F Activities at Transformation
To define a malignant T-cell oncogenic program in tCTCL, we first performed differentially expressed gene (DEG) and gene set enrichment pathway analyses comparing malignant T cells in TT and benign CD4+ T cells in the TIME (Fig. 4A). Computational overlap of genes ranked by expression fold change with the Molecular Signatures Database (MSigDB) hallmark gene sets showed significant enrichment of genes in OXPHOS, MYC, EMT, and E2F target pathways and downregulation of IFNα, TNFα, and IFNγ (Fig. 4A; Supplementary Table S12A). Importantly, significant upregulation of OXPHOS, MYC, and E2F target pathways and downregulation of IFNα, TNFα, and IFNγ pathways were also observed in disease evolution from PP to TT (Fig. 4B; Supplementary Table S12B). Further analysis of the DEGs between malignant T cells in TT versus benign CD4+ T cells highlights a 55-gene malignant T-cell signature with upregulation of genes involved in OXPHOS (SLC25A5 and NDUFB2), EMT/cellular plasticity and stemness (TWIST1, EPCAM, and CLDN7), MYC signaling (NME2, SRM, and PIM3), chemokines (CCR7), cytokines (MIF), cell migration and motility (CD9/tetraspanin, LAIR2, and RAB25), as well as downregulation of HLA-A, HLA-B, and HLA-F (MHC I and IFNγ pathway) that is suggestive of immune escape in tCTCL (Fig. 4C; Supplementary Table S13). Of the upregulated genes, TWIST1 (cellular plasticity/stemness), NDUFB2 (OXPHOS), and LGALS3 (galectin 3/RAS) are also predicted drivers within the GISTIC2.0 amplification loci and are therefore considered high confidence drivers in tCTCL.
We next examined the benign immune cell types in the tCTCL TIME. Although there was no quantitative difference in benign CD4+ T cells, CD8+ T cells, macrophages, NK cells, and dendritic cells between concurrent PP and TTs, B cells were significantly enriched in TTs (Supplementary Fig. S7) and may play a protumorigenic role during disease evolution.
Although scRNA-seq data offer a high-dimensional view of cell types, cell states, and malignancy at single-cell resolution, these measurements ultimately represent a snapshot in time of a complex ecosystem with heterogeneous cell substrates in a diverse range of evolutionary stages. We therefore sought to order the trajectory of tCTCL constituents in pseudotime using Monocle 3, a machine learning algorithm that constructs the trajectory of cells between one of several possible end states along geodestic distance from a root node and can learn trajectories that have loops or points of convergence (51). We first arbitrarily assigned benign T cells as the starting root node, here CD4+ naïve T cells, and demonstrated a monodirectional trajectory from benign to malignant CD4+ T cells in PP and TT (Fig. 4D, top). We next tracked the kinetics of the 55-gene malignant T-cell program along pseudotime and observed coordinated upregulation of CXCL13 (chemokine), SLC25A5 (OXPHOS), and NME2 (MYC) from benign to malignant PP and accentuation of cellular plasticity genes such as EPCAM, TWIST1, and PTHLH at the end of the trajectory in TT (Fig. 4D bottom; Supplementary Fig. S8; Supplementary Table S14). Notably, downregulation of MHC I genes (HLA-A, HLA-B, and HLA-F) occurred early at the bifurcation from benign to malignant T cells in PP (Fig. 4D), suggesting immune evasion as an early event in CTCL lymphomagenesis, and is exacerbated from benign CD4+ to malignant T cells in PP to malignant T cells in TT (Fig. 4E).
Small-Molecule Inhibitors of OXPHOS and MYC Suppress CTCL Cells In Vitro
Our scRNA-seq analyses identified OXPHOS and MYC as the most enriched pathways in tCTCL and potential therapeutic vulnerabilities. High OXPHOS is a hallmark of multiple hematopoietic malignancies and solid tumors (33, 52–54). In preclinical models of acute myeloid leukemia and mantle cell lymphoma, inhibition of OXPHOS by IACS-010759, a potent small-molecule inhibitor of mitochondrial electron transport chain complex I, results in inhibition of tumor cell proliferation and induction of apoptosis (33, 54). Likewise, MYCi975, a small-molecule MYC inhibitor, is shown to have potent preclinical in vivo activity against MYC (34). To determine the sensitivity of CTCL and other T-cell lymphomas to OXPHOS and MYC inhibition, we performed pharmacologic perturbation assays using IACS-010759 and MYCi975 in MF and other TCL cell lines (Fig. 4F; Supplementary Table S15A and S15B). Using a group of established T-cell lymphoma cell lines, including Myla (MF), HH (MF, leukemic phase), and Hu78 (SS), we found low nanomolar sensitivity to IACS-010759 as assessed by apoptosis, particularly for Myla, and low micromolar sensitivity to MYCi975, as assessed by cell proliferation, cell viability, and cytotoxicity (MTS assay). Together, these findings provide functional support for the two most significantly enriched malignant T-cell programs in tCTCL and potential use of OXPHOS and MYC inhibitors in treating aggressive T-cell lymphomas.
Analysis of Cellular Cross-talk between Malignant T Cells and the TIME Highlights MIF–CD74 Interactions
To systematically disentangle the complex intercellular network within the tCTCL ecosystem, we next investigated our scRNA-seq data for potential receptor–ligand interactions between malignant T cells and the TIME. We integrated our scRNA-seq data with CellPhoneDB, a computational tool that interrogates scRNA-seq data for protein subunit structures and receptor–ligand interactions, and explored the cross-talk between malignant T cells and the macrophages/monocytes, B cells, dendritic cells, NK cells, fibroblasts, and endothelial cells (Fig. 5A; Supplementary Table S16A–S16C). We showed top-ranking receptor–ligand pairs between MIF (ligand) in malignant T cells and CD74, also known as HLA-DR antigen-associated invariant chain, in macrophages/monocytes, B cells, and dendritic cells (Fig. 5A). Importantly, our scRNA-seq 55-gene malignant T-cell program also highlighted MIF as one of the significant differentially upregulated genes in malignant T cells (Fig. 4C).
To demonstrate in situ evidence for MIF–CD74 interactions, we built a tissue microarray (TMA; 80 tissue cores) using skin biopsies from 21 patients with tCTCL (16 with matched PP-tumor LCT; 64 cores; Fig. 1A) and 9 patients with PP without any history of transformation (PP; 16 cores). We performed immune profiling by Vectra mIF, and regions of interest (ROI) in each biopsy were selected from all lymphoid-dense areas in the biopsy. We stained this TMA with two panels of antibodies (anti-MIF, CD74, CD3, CD8, Ki-67, and CD68 antibodies and DAPI; anti-MIF, CD74, CD3, CD8, Ki-67, PAX5 antibodies and DAPI; Supplementary Table S17A and S17B) and imaged on Vectra 3. We analyzed a total of 488,399 cells by HALO software (Indica Labs) and demonstrated colocalization between MIF in malignant CD4+ T cells (CD3+ CD8− Ki-67+ cells) and CD74 in macrophages (CD68+; Fig. 5B) and B cells (PAX5+; Fig. 5C), providing in situ validation for the CellPhoneDB-predicted receptor–ligand interactions. These findings together raise the possibility that MIF could be an important therapeutic target in tCTCL through inhibiting MIF–CD74 interactions between malignant T cells and APCs. Consistent with our observation in scRNA-seq data set, we also found moderate to prominent B-cell infiltrate in TT, whereas the PP lesions showed either minimal or no B-cell infiltrate (Fig. 5D; Supplementary Fig. S7). This TMA was also evaluated for potential differences in the TIME by CD30 expression (≤1%, >1% to <10%, 10%–50%, >50%) via IHC studies. Here, we did not detect quantitative differences in macrophage, B-cell, and CD8+ T-cell infiltrates, or MIF–CD74 interactions between malignant T cells and macrophages or B cells (Supplementary Fig. S9A–S9F).
Dominant Subclones in tCTCL Show Upregulation of Genes Encoding RPSs
We have observed interpatient tumoral heterogeneity in tCTCL (Fig. 3C), with some patients revealing more than one malignant T-cell cluster. We therefore sought to interrogate intratumoral heterogeneity (ITH) at the genetic and transcriptional levels and determine if ITH correlates with patient outcomes. As CNV inference has been used to distinguish genetically distinct tumor subclones (32, 55–57), we focused our attention on malignant T cells from the scRNA-seq data set and performed CNV inference at the level of tumor subclusters using the inferCNV “subclusters mode” (qnorm method; Fig. 6A). We identified genetically distinct subclones of malignant T cells within TTs from patients 11, 35, and 50 despite TCR monoclonality. Importantly, of the 8 patients in the scRNA-seq cohort, PT11 and PT35 showed highest tumor burden (PT11 TT clinical image; Fig. 1A) and the worst clinical outcomes. UMAP clustering further demonstrated subclonal transcriptional heterogeneity in TTs from the same patients (PT11 clusters 1, 11, 12; PT35 clusters 0, 3, 9; PT50 clusters 4, 8; Fig. 6B and C; Supplementary Table S18A and S18B), with dramatically elevated ribosomal gene expression (large and small RPSs) in the preponderant malignant T-cell subclones (PT11 and PT35). A similar ribosomal signature has been observed in circulating tumor cells from highly aggressive breast cancer (58). Our findings suggest potential therapeutic utility in targeting the ribosomal synthesis and global eukaryotic translational machineries for patients with aggressive, transformed CTCL.
Cutaneous tCTCL and SS Have Distinct Oncogenic Malignant T-Cell Programs
Finally, as we observed distinct patterns of CNA in tCTCL and SS (Fig. 2), we sought to investigate the differences in T-cell oncogenic signatures between these two advanced stages of CTCL at single-cell resolution. A recently reported SS cohort with matched single-cell transcriptomic [expanded CRISPR-compatible cellular indexing of transcriptomes and epitopes by sequencing (ECCITE-seq)] and single-cell TCR clonotype (Herrera cohort; ref. 32) was downloaded for this analysis. Malignant T cells from this SS cohort were identified by the dominant TCR clonotype and benign CD4+ T cells by CD4+ polyclonal T cells (Supplementary Methods), and UMAP clustering showed clear separation between malignant T cells versus benign CD4+ T cells (Supplementary Fig. S10A). Interestingly, malignant T cells in SS do not show upregulation of OXPHOS and MYC activities (Supplementary Fig. S10B) or downregulation of MHC I genes compared with benign CD4+ T cells (Supplementary Fig. S10C; Supplementary Table S19A). We next compared the gene-expression profile between malignant T cells in tCTCL (PP and TT) with malignant T cells in SS. Here, the malignant T cells cluster by patient of origin (Fig. 7A), with tCTCL patient clusters clearly separating from the SS patient clusters (Fig. 7B), and with benign T cells across different patients clustering together (Supplementary Fig. S11). Importantly, malignant T cells in tCTCL (PP and TT) show significant upregulation of OXPHOS, MYC, EMT/stemness, E2F target, and TNFα activities compared with malignant T cells in SS, as well as downregulation of IFNα pathway genes (Fig. 7C; Supplementary Table S19B). Violin plots of MHC I gene expression in malignant T cells further revealed significantly lower HLA-A, B, and C expression in tCTCL TT and PP compared with SS (Fig. 7D), suggesting immunosurveillance escape through loss of MHC I as a more prominent feature in tCTCL.
DISCUSSION
Transformed CTCL is an aggressive large cell lymphoma that is resistant to multiple forms of systemic therapy. There is a critical unmet need for identifying novel therapeutic targets for this incurable cancer, yet its disease biology is poorly understood due to the multitude of challenges associated with tissue-based research in rare cancers. Genomic- and transcriptomic-level information about CTCL at disease transformation is limited, and a better understanding of tCTCL disease biology has broad implications for therapy. Our study contains the largest collection of tCTCL clinical specimens that are difficult to obtain and provides an important resource for the study of tCTCL biology and the identification of novel therapeutic vulnerabilities. We generated a WES data set of similar size to other TCGA rare cancers, the first scRNA-seq atlas of tCTCL, and immune profiled the TIME by scV(D)J-seq and mIF using serial biopsies from patients with tCTCL. Together, this multiomics study uncovered key driver events, prognostic mutational signatures, diverse oncogenic programs, and potential receptor–ligand interactions that malignant T cells exploit at LCT.
Defining the T-cell lymphoma ecosystem at single-cell resolution is particularly challenging, as the TIME is comprised of malignant CD4+ T cells, benign T cells, and other diverse immune cells. Gene expression–based clustering is imprecise in defining malignancy, and dropout read is a well-known phenomenon in scRNA-seq and scVDJseq data sets. Here, we implemented a robust two-tiered algorithm using simultaneous single-cell whole transcriptome-V(D)J profiling and chromosomal copy-number inference to identify malignant T cells from the TIME, and we believe this methodology can be applied to the study of other T- and B-cell leukemia/lymphomas. Although CTCL is clinically viewed as a cancer of monoclonal T cells, a recent study showed TCR clonotypic diversity in CTCL by bulk TCR-seq and suggested that T-cell tumorigenesis occurs prior to TCR gene rearrangement (59). In tCTCL, we observed TCR monoclonality, with monoclonal T cells showing malignant CNV patterns, and background polyclonal T cells in the TIME showing neutral CNV patterns. The difference in these two studies could be due to differential resolution of TCR detection using scV(D)J-seq, with the ability to assign precise TCR α- and β-sequence pairs to individual T cells. Nevertheless, our data suggest that TCR monoclonality is a key feature of transformation. Future single-cell TCR studies involving CTCL lesions at different stages of disease can potentially elucidate the origin of CTCL with respect to the timing of TCR gene rearrangement.
A central finding from our single-cell study is the diverse oncogenic programs that malignant T cells exploit to confer aggressive behavior and survival advantage at transformation, with upregulation of OXPHOS and MYC as the top enriched pathways, which are progressively upregulated in disease evolution, followed by EMT/stemness and E2F target genes (Fig. 4A and B). Notably, similar OXPHOS and MYC upregulation was not observed in our single-cell analysis of an independent SS cohort (Supplementary Fig. S10B; Fig. 7C; ref. 32). Metabolic reprogramming toward OXPHOS, and hence increased ATP production, has been reported in other high-grade leukemias, lymphomas, and solid tumors (33, 52–54). In tCTCL, this metabolic shift can be viewed as advantageous for supplying the energy demand of rapidly growing, high-grade lymphoma cells. As MYC activation has been shown to drive cancer cell metabolism toward glutaminolysis and mitochondrial biogenesis (60, 61), the shift toward OXPHOS can be related to MYC upregulation. Interestingly, quiescent T cells (naïve or memory T cells) are known to use catabolic metabolism by OXPHOS to fuel cell survival (62), raising the possibility that the upregulation of OXPHOS in tCTCL can also be related to T-cell state at transformation. Although future studies are needed to chart the mechanistic action of these pathway components, it is encouraging that novel small-molecule inhibitors of OXPHOS (IACS-010759) and MYC (MYCi975) provided functional support of OXPHOS and MYC as potential therapeutic vulnerabilities. Lastly, in addition to OXPHOS and MYC, our scRNA-seq data set also demonstrated significant upregulation of MIF in malignant T cells and interactions between malignant T cells expressing MIF and macrophages and B cells expressing CD74. These findings raise the possibility of MIF as potential therapeutic targets in tCTCL.
A recurrent theme observed in this study is the upregulation of genes involved in cellular plasticity and stemness in tCTCL. In our WES data set, GISTIC analyses aimed at detecting genomic loci targeted by SCNAs revealed recurrent amplifications in TWIST1, a master transcription factor of EMT/stemness, whereas the scRNA-seq 55-gene malignant T-cell signature demonstrated upregulation of TWIST1, CLND7, and EPCAM, all of which are genes involved in cellular plasticity and stemness. We also observed recurrent mutations in Hippo pathway genes such as FAT1 (∼30%), loss of which has been implicated in driving tumor cell stemness and metastasis in cutaneous squamous cell carcinoma (46). Tumor plasticity is a key feature of aggressive cancer transition, which can certainly account for the transition of CTCL to an aggressive large cell lymphoma with blast-like morphology at transformation. In large-scale gene-expression studies across 21 solid tumors, activation of stemness programs was shown to positively correlate with ITH, drive clonal evolution, and limit antitumor immune responses (63). Indeed, we also observed ITH in tCTCL, with prominent subclonal upregulation of ribosomal subunit gene expression in patients with the worst clinical outcomes.
Besides the upregulated signatures, one intriguing observation from our single-cell data set is the downregulation of MHC I in malignant T cells in tCTCL. Comparison of MHC I gene-expression levels (HLA-A, C, E, F) showed progressive and significant decreases from benign CD4+ T cells to PP to TT, and our pseudotime analysis corroborated these changes in expression across disease evolution. Although high TMB cancers are theoretically more immunogenic and responsive to T cell–based immunotherapies, loss of MHC I antigen presentation can compromise the visibility of tumor cell antigens to CD8+ T cells and response to checkpoint blockade agents. In such MHC I–negative or low tumors with defective antigen presentation machinery, NK cell–based therapies and epigenetic drugs can be considered.
The role of B cells in the tumor microenvironment is recently under intense investigation in various cancer types, though it is still incompletely explored compared with other benign immune substrates of the TIME. The results are also variable, with some studies showing positive correlation between B-cell infiltrate and patient outcomes or response, whereas other studies suggest a protumorigenic role of B cells (64). In tCTCL, our scRNA-seq data show significant enrichment of B-cell infiltrate mostly in TT, a finding that is supported by multiplex IF. Altogether, our results suggest a potential protumorigenic role of B cells in tCTCL, possibly through interactions with malignant T cells via MIF–CD74 signaling.
Lastly, CTCL is a rare cancer with well-known racial disparity. Black/AA patients have a higher incidence rate, younger age of onset, higher disease burden, and inferior survival compared with whites, even after accounting for disease characteristics, socioeconomic factors, and treatments received (1, 2). Nonetheless, potential biological factors underlying these racial disparities are poorly understood. Here, we provide the first attempt at identifying potential genomic correlates for the survival gap between Black/AA and white patients. Although the sample size is admittedly small, we observed a significantly lower contribution of UV signatures that are prognostic for favorable survival and enrichment of other signatures that may drive worse outcomes in Black/AA patients. Future studies involving larger sample sizes from the vulnerable population and research into their TIME, genomic correlates, and potential therapeutic targets will hopefully help reduce racial disparity in CTCL.
Rare cancers represent one of the greatest inequalities in cancer research, with the lack of well-curated tissue specimens to study disease biology and resources for pharmaceutical development. Transformed CTCL exemplifies such challenge where the lack of genomic- and transcriptomic-level information contributes to the paucity of drug-discovery efforts. In this study, we have identified potential genetic driver events and pathways in tCTCL and elucidated a malignant T-cell oncogenic program with potential novel therapeutic vulnerabilities. Although further validation in larger cohorts and preclinical models is needed, our investigation provides a key resource with the largest collection of tCTCL samples studied to date. We anticipate results from this study can be extrapolated to other T-cell lymphomas and will help usher in novel immunotherapeutic strategies to combat this currently incurable cancer.
METHODS
A full description of all methods for WES (DNA extraction and library prep, data processing and somatic mutation filtering, driver gene detection, and mutational signature analysis), scRNA-seq (data processing and benign cell cluster annotation), Vectra mIF studies, and comparison of single-cell transcriptomic profile to SS (Herrera cohort) can be found in the Supplementary Methods.
Patient Cohort and Clinical Annotation
A total of 56 patients with tCTCL were included in this study. The studies were conducted in accordance with the Declaration of Helsinki and approved by the Moffitt Cancer Center (MCC) Institutional Review Board (IRB; MCC19672). Electronic medical charts were reviewed independently by two investigators to document the clinical parameters, including age, gender, time of first rash, time of initial MF diagnosis, time of LCT, time of death/last follow-up, clinical staging at time of initial MF diagnosis or presentation to MCC, clinical restaging at time of LCT, and treatment history (Supplementary Table S1).
WES
Clinical Tissue Samples
Of the 56 patients with tCTCL, 54 had available tissues for WES (n = 70 skin formalin-fixed, paraffin-embedded (FFPE) biopsies; Fig. 1A). For WES and immunofluorescence studies, standard-of-care skin biopsies (FFPE) from 54 patients with tCTCL were analyzed under IRB-approved protocol (MCC19672) and approved waiver of consent where applicable. All samples were assessed by two pathologists for adequacy of lymphoid infiltrate, including 45 TT and 25 PP lesions (9 with LCT). Seventeen patients had paired PP and TT (11 precursor PP, 6 concurrent PP to the LCT lesion). Matched germline samples were available from 41 patients, while remaining patient samples used panel of normal as control. Tumor samples were microdissected before being submitted for WES.
Cancer Driver Gene Detection
Next-generation sequencing data were processed using the established in-house bioinformatics pipeline (65). To enhance the coverage of our analysis, two mutation-based computational algorithms were applied, MutSigCV (version 1.41; ref. 44) and dNdScv (RRID: SCR_017093; ref. 43). MutSigCV was run with default parameters using the MATLAB (RRID:SCR_001622) Compiler Runtime (MCR). A global q-value threshold of 0.1 was applied for selecting driver genes as suggested by the MutSigCV developers. dNdScv was run in R (version 3.6.3) using the default parameters. A global q-value (qallsubs_cv) of <0.05 was used as a cutoff for selecting predicted driver genes.
CNV Detection and Analysis
CNV analyses were performed in samples with matched germline. CNVs were identified using R package ExomeLyzer (54). Copy-number change was determined as the log2 ratio of tumor versus normal reads, and a circular binary segmentation algorithm was applied for CNV segmentation (66). To identify genomic regions significantly enriched with CNVs, GISTIC2.0 (RRID:SCR_000151; ref. 47) was used with default settings. Genomic regions with FDR < 0.25 were considered as significantly enriched and highlighted in the chromosome landscape plot.
TMB
TMB analysis was performed in all samples with matched germline/normal tissue using Maftools (v2.6.0; ref. 67). Results were compared with 33 TCGA cancer types (MC3 data; ref. 68). TMB was calculated by the following formula:
Results were compared with 33 TCGA cancer types (MC3 data; ref. 68) and two SS cohorts (Wang cohort and Choi cohort; refs. 26, 36). The Choi cohort raw sequencing data were downloaded and processed using the same computational pipeline as our tCTCL cohort. The Wang cohort TMB was calculated using processed mutations from the original publication. The sequencing capture size is 39 Mb for the study cohort and 35.8 Mb for the TCGA cancer cohorts after data harmonization.
Mutational Signature Analysis
The repertoire of mutational signatures was cataloged at both the sample and cohort levels. At the sample level, contributions of known signatures (COSMIC mutation signatures version 3) in each tissue were independently calculated using two algorithms, MutationalPatterns (39) and deconstructSigs (38). Both tools were run using the default setting. To decipher the dominant mutagenic processes at the cohort level, de novo representative signatures of the cohort were extracted using Maftools (67). A mutation context matrix was first built, and the top three mutational signatures that optimally represent the mutation profile were extracted using the NMF method. Finally, the constructed signatures were compared against COSMIC version 3 signatures (37).
Survival Analysis
For Kaplan–Meier survival analysis, survival R package (version 3.2-7) was used. Survival analysis was performed from time of initial LCT to the date of death (event time) or the date of last follow-up (censoring time). The patients were grouped into two subgroups according to the median of genomic features, such as the weighted contribution of SBS7 signature. For patients with paired PP and TT tissues, data from the TT sample were used for survival analysis. Survival plot was generated using the survminer R package.
Analysis of Highly Mutated Pathways
A variety of cancer-associated pathways were assessed for mutational enrichment (Supplementary Table S6), including eight oncogenic signaling pathways summarized by the TCGA project and implemented in Maftools (69). Genes with nonsynonymous mutations observed in ≥3 tissues and recurrent in >10% of the samples were plotted in oncoplot using Maftools.
scRNA-seq and V(D)J-seq
Patient Sample Collection
Fresh tissue specimens were collected from eight patients with tCTCL treated and managed at the H. Lee Moffitt Cancer Center Cutaneous Oncology and Malignant Hematology clinics following written informed consent under the IRB-approved Total Cancer Care Protocol (MCC14690) and processed under MCC19672. Two samples were collected per patient (concurrent PP and TT lesions).
Single-Cell V(D)J-seq for TCR Clonotype
TCR annotation was performed using Cell Ranger “vdj” function for sequence assembly and paired clonotype calling. Paired TCR α/β reads attributed to each cell barcode were grouped and assembled into a single contig to determine the combined TCR α/β clonotype, and contigs that were predicted as productive were selected for the downstream analyses.
Single-Cell CNV Inference
InferCNV v1.3.5 (https://github.com/broadinstitute/inferCNV, Trinity CTAT Project) was used to extract chromosomal CNV patterns. Two thirds of nondominant TCR clonotype T cells (CD4+ and CD8+) were used as reference normal cell group for denoise control. For the observation cell group, (i) the remaining nondominant TCR clonotype cells, (ii) all dominant TCR clonotype (“true malignant”) cells, and (iii) all CD3+ T cells with TCR dropout (“malignant suspect”) were included. InferCNV was run using “denoise” mode, and a “cluster by group” parameter was turned off such that observation cells can cluster regardless of tissue of origin.
Differential Expression Testing and Gene Set Enrichment Analysis
Differential expression testing was performed using the Wilcoxon rank sum method with two Seurat functions: FindMakers for pairwise comparisons and FindAllMarkers for comparisons across multiple groups. Both functions were run on the “RNA” data assay with no prefilters. Significant DEGs were filtered with a q-value < 0.05 and an absolute value of fold change (FC) >2 or <0.5. To reveal pathways in which DEGs are enriched, gene set enrichment analysis was performed on all of the detected genes ranked by average logFC. The analysis was performed using the clusterProfiler R package (RRID:SCR_016884; ref. 70) against predefined hallmark gene sets of the MSigDB database.
Pseudotime Trajectory Analysis
T cells, both malignant and benign, were first constructed in a Seurat object and processed and clustered as described above. Trajectory analysis was performed using Monocle 3 (51, 71–74). To meet the format requirement of Monocle 3, Seurat object and clustering information was converted as a cell data set using SeuratWrappers R package. Trajectory was then learned using the default Monocle 3 parameter “learn_graph” function. Naïve CD4+ T cells in the graphical interface was manually selected as the root node and cells were ordered in pseudotime. Next, a branch of the trajectory graph originating from naïve CD4+ T cells to malignant T cells was selected using the “choose_graph_segments” function. The kinetics of 55-gene malignant T-cell program were plotted along the pseudotime of the selected branch using the “plot_genes_in_pseudotime” function of Monocle 3.
CellPhoneDB Intercellular Communication Analysis
We integrated our scRNA-seq data set with CellPhoneDB (RRID:SCR_017054; ref. 75) and explored intercellular receptor–ligand interactions between malignant T cells and macrophages/monocytes, B cells, dendritic cells, NK cells, fibroblasts, and endothelial cells. The raw count matrix of single-cell expression profile was normalized, and the count value of each gene in each cell was divided by the total counts of corresponding cell and multiplied by 10,000. CellPhoneDB was run with the “statistical” method in the virtualenv environment using default parameters. Receptor–ligand pairs with P value < 0.05 were filtered as significantly enriched pairs. To visualize the results, we selected the top three ligand–receptor pairs for each cell type pair (in log2 mean expression value) and summarized the results in a dot plot using R package ggplot2 (RRID:SCR_014601).
Analysis of Intratumoral Heterogeneity
InferCVN subcluster analysis for ITH at the genetic-level was performed on malignant T cells (identified based on scV(D)J-seq and copy-number inference method as described above). Patient-specific clustering was restricted by setting the parameter “cluster_by_groups” as “True.” “Subclusters” mode was enabled with the qnorm partition method to reveal ITH patterns. Malignant cells were clustered based on similarity of the extracted CNV patterns. Gene expression–based clustering of malignant T cells by UMAP was performed to evaluate for ITH at the transcriptional level.
Independent SS Cohort Analysis (Herrera Cohort)
To perform single-cell transcriptomic analysis of SS and malignant T-cell comparison between tCTCL and SS, ECCITE-seq scRNA-seq gene-expression counts and TCR α/β clonotype profiles from six SS blood samples (SS1–SS6) were downloaded from the NCBI Gene Expression Omnibus database (GSE171811; ref. 32; Supplementary Methods).
Pharmacologic Perturbation Assays
Established T-cell lymphoma lines were used for drug assays, including Myla (MF; Sigma-Aldrich; cat. #95051032), HH (MF—leukemic phase; ATCC; cat. #CRL-2105), and Hu78 (SS; ATCC; cat. #TIB-161), Jurkat (ATLL), and MJ (ATLL; ATCC; cat. #8294). For the apoptosis assay, indicated cell lines were seeded with a density of 1 × 105 cells/well in 96-well plate and treated with 8 nmol/L of IACS-010759 (SelleckChem #S8731). At day 5, cells were harvested and stained with Annexin V and PI following the manufacturer's protocol (BioLegend; cat. #640914). Annexin V+PI+ population was gated on singlet population using FlowJo 10 software (RRID:SCR_008520). For cell proliferation assay, indicated cell lines were seeded with a density of 1 × 104 cells/well in 96-well plate and treated with different dose of MYCi975 (MedChemExpress; cat. #HY-129601) for five days. At day 5, cells were incubated with MTS reagent following the manufacture's protocol (Promega; cat. #G3580). Absorbance at OD 490 nm was recorded, and the percentage of growth was normalized to vehicle control. IC50 was calculated based on the curve fitting result using nonlinear regression function of GraphPad Prism 8 (RRID:SCR_002798).
Vectra Multiplex Immunofluorescence
An 80-core TMA was built using skin biopsies from 21 patients with tCTCL (64 tissue cores, including 16 patients with matched PP and TT) and additional 9 patients with PP lesions and no history of LCT (16 tissue cores). ROIs were selected on the basis of the most lymphoid-dense/representative area in the tissue. Quantitative image analysis was performed using the HALO Image Analysis Platform (Indica Labs).
Statistical Analysis
Statistical analysis was performed using R programming language (version > 3.6). In survival test, a two-tailed log-rank test was used to determine the difference between patient subgroups. The statistical significance between category groups was determined using the Wilcoxon rank sum test. ns, P > 0.05; *, P < 0.05; **, P < 0.01; ***, P < 0.001; ****, P < 0.0001. For comparison of survival between patient groups (Fig. 1D; Supplementary Fig. S3B), data from the TT sample were used when a patient had paired PP and TT available.
Data Availability Statement
The sequencing data, including WES, scRNA-seq, and scVDJ-seq, have been submitted to the NCBI BioProject database PRJNA754592. The leukemic/SS WES data analyzed in this study were obtained from the NCBI BioProject database PRJNA285408.
Authors’ Disclosures
L. Seminario-Vidal reports grants and personal fees from Eli Lilly, Helsinn, Boehringer Ingelheim, and Kyowa Kirin, grants from Soligenix, Eisai, Novartis, AbbVie, Bristol Myers Squibb, Celgene, Glenmark, Amgen, AnaptysBio, and Innate Pharma, and personal fees from Regeneron and Blueprint outside the submitted work. S. Whiddon reports personal fees from Kyowa Kirin outside the submitted work. J.R. Conejo-Garcia reports personal fees from Alloy Therapeutics and grants and personal fees from Anixa Biosciences outside the submitted work, as well as a patent on targeting CD277 for Compass Therapeutics pending to The Wistar Institute and a patent on targeting FSHR for Anixa Biosciences issued to The Wistar Institute. L. Sokol reports grants from the NIH/NCI during the conduct of the study, as well as personal fees from Kyowa Hakko Kirin Co., Ltd., Dren Bio, and Daiichi Sankyo outside the submitted work. P. Chen reports a patent for methods and kits for the detection of malignancies pending. No disclosures were reported by the other authors.
Authors’ Contributions
X. Song: Conceptualization, resources, data curation, software, formal analysis, validation, investigation, visualization, methodology, writing–original draft, writing–review and editing. S. Chang: Data curation, formal analysis, validation, investigation, visualization, methodology, writing–original draft, writing–review and editing. L. Seminario-Vidal: Resources, data curation, writing–review and editing. A. de Mingo Pulido: Data curation, validation, investigation, writing–review and editing. L. Tordesillas: Data curation, validation, investigation, writing–review and editing. X. Song: Data curation, software, investigation. R.A. Reed: Resources. A. Harkins: Resources. S. Whiddon: Resources. J.V. Nguyen: Data curation, software, formal analysis, validation, investigation, visualization, methodology, writing–original draft, writing–review and editing. C. Moran Segura: Data curation, investigation. C. Zhang: Data curation. S. Yoder: Data curation, supervision. Z. Sayegh: Resources, data curation. Y. Zhao: Resources, data curation. J.L. Messina: Resources, writing–review and editing. C.M. Harro: Resources. X. Zhang: Resources. J.R. Conejo-Garcia: Resources, writing–review and editing. A.E. Berglund: Writing–review and editing. L. Sokol: Resources, writing–review and editing. J. Zhang: Resources, data curation, software, supervision, project administration. P.C. Rodriguez: Resources, data curation, formal analysis, supervision, investigation, visualization, methodology, writing–review and editing. J.J. Mule: Conceptualization, resources, supervision, funding acquisition, writing–review and editing. A.P. Futreal: Resources, software, supervision, writing–review and editing. K.Y. Tsai: Conceptualization, resources, supervision, writing–review and editing. P.-L. Chen: Conceptualization, resources, data curation, formal analysis, supervision, funding acquisition, validation, investigation, visualization, methodology, writing–original draft, project administration, writing–review and editing.
Acknowledgments
We are grateful to the patients and their families for contributing to this work. This research was made possible through the Total Cancer Care Research Protocol at the H. Lee Moffitt Cancer Center and Research Institute and supported by the Tissue Core, Molecular Genomics Core, and Advanced Analytical and Digital Laboratory. We thank Dr. Jennifer Wargo for her very helpful reading of the manuscript, Christophe Georgescu for his suggestions on implementing the inferCNV analytics, and Dr. Whijae Roh for his thoughtful comments on the genomics data. This research was supported by generous philanthropic contributions to H. Lee Moffitt Cancer Center Moffitt Foundation, Moffitt Clinical Science Fund, Miles for Moffitt fund, Department of Pathology, and Department of Cutaneous Oncology research support fund to P.-L. Chen. Support for shared resources was provided by a Cancer Center Support Grant (CCSG) P30-CA076292 to H. Lee Moffitt Cancer Center. K.Y. Tsai acknowledges support of Donald A. Adam Melanoma and Skin Cancer Center of Excellence and Mr. and Mrs. Winston Weber. J.R. Conejo-Garcia and L. Sokol acknowledge support of this work by NIH R01CA240434.
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.