The molecular drivers of classic Hodgkin lymphoma (HL) underlying its pathologic and clinical distinctions from non-Hodgkin lymphoma (NHL) have been challenging to uncover. The study of HL has been impeded by the rarity of the Hodgkin-Reed-Sternberg (HRS) cells within a dense tumor microenvironment (TME). Previously, we developed a FACS approach to isolate the HRS cells and profiled their whole exomes and genomes (Reichel et al 2015; Maura et al 2023). These studies did not reveal an HRS-specific mutational signature, but rather mutations that overlap with those of NHL, such as those in SOCS1 and HIST1H1E. These data suggested that non-genetic factors such as the underlying cell state may cooperate with somatic mutations to induce HL. However, we previously lacked the ability to co-map somatic mutations and cell states within the same cancer cells. We thus developed Genotyping of Transcriptomes (GoT; Nam et al 2019) that links single-cell RNA-seq (scRNA-seq) and somatic mutations within the same thousands of individual cells. We applied GoT to primary HL samples (n = 10 HL; 68291 cells, Fig. 1A) after FACS-isolation of HRS cells. We analytically integrated the cells across the samples based on gene expression and identified the expected cell types including B-cells, T-cells, monocytes and HRS cells (Fig. 1B). As expected, the HRS cells displayed low expression of B-cell antigens ( CD19, CD79A) and high levels of TNFRSF8 (CD30), CCL17, CCR7 and PD-L1 (Fig. 1C). GoT data revealed that mutations were restricted to the HRS cells (Fig. 1D). Copy number variation (CNV) analysis further confirmed that CNVs were specific to the HRS cells (Fig. 1E). As scRNA-seq has unveiled marked intratumoral cell state heterogeneity across cancers, we posited that the HRS cells may display heterogeneous cell states. Indeed, we identified distinct cell states (Fig. 1F), including cycling HRS and metabolically-active HRS cells upregulating the translation machinery ( NACA, RACK1). We also observed a cell state with elevated inflammatory signals ( TNF, JUN, FOS, CCL22). One group demonstrated a partially preserved B cell program with high CD19 and CD79A. Unexpectedly, we also identified an HRS state with features of neuronal-glial (NG) transdifferentiation, including genes involved in synapse formation ( NRG2, NRXN3). Consistently, genes upregulated in the NG state were enriched in nervous system and neuronal development pathways (Fig. 1G). We further assessed NG gene expressions across normal cell types and confirmed their specificity to NG lineage cells (Fig. 1H). These findings were consistent with dendritic processes of HRS cells highlighted by CD15 and CD30 (Fig. 1I, red arrows). Areas of CD15/CD30 staining devoid of nucleus (Fig. 1I, blue arrows) corroborated the presence of dendritic processes. These results suggested that the NG differentiation may be a key defining feature of HRS cells, whereby dendritic processes facilitate a tight control of the TME. Consistently, we identified strong HRS-immune complexes that persisted through tissue processing and sorting (Fig. 1F). To determine the degree of plasticity or heritability of the HRS cell states, we integrated VDJ sequencing with GoT (GoT-VDJ), whereby somatic hypermutations in the VDJ served as endogenous barcodes. Phylogenetic reconstructions revealed that HRS cell states were distributed across clades (Fig. 1J), unveiling a plastic interchange of HRS cell states. Consistently, GoT showed that mutant cells were present across the cell states (Fig. 1K), indicating that the cell states were decoupled from the genetic identities. We therefore postulated that the NG program may be encoded epigenetically. We performed single nuclei ATAC-seq on HRS-enriched HL samples (n = 5 samples; Fig. 1L). Compared to admixed B cells, the HRS cells displayed enhanced accessibility for binding motifs of the NG lineage transcription factors (TF), e.g. OLIG1/2 and SOX family (Fig. 1M), revealing the TF networks that govern the lineage infidelity program. We further identified that the NG program is cell intrinsic and at least partially stable through in vitro passaging, as 4 of 5 HL cell lines exhibited the highest NG gene expression compared to NHL cell lines (Fig. 1N). Altogether, these single-cell multi-omics studies indicated that an epigenetically encoded NG program may cooperate with lymphoma driving mutations to drive the development of HL, paving novel avenues of therapy.
Read full abstract