Distinct Hodgkin lymphoma subtypes defined by noninvasive genomic profiling.

The scarcity of malignant Hodgkin and Reed-Sternberg cells hampers tissue-based comprehensive genomic profiling of classic Hodgkin lymphoma (cHL). By contrast, liquid biopsies show promise for molecular profiling of cHL due to relatively high circulating tumour DNA (ctDNA) levels1-4. Here we show that the plasma representation of mutations exceeds the bulk tumour representation in most cases, making cHL particularly amenable to noninvasive profiling. Leveraging single-cell transcriptional profiles of cHL tumours, we demonstrate Hodgkin and Reed-Sternberg ctDNA shedding to be shaped by DNASE1L3, whose increased tumour microenvironment-derived expression drives high ctDNA concentrations. Using this insight, we comprehensively profile 366 patients, revealing two distinct cHL genomic subtypes with characteristic clinical and prognostic correlates, as well as distinct transcriptional and immunological profiles. Furthermore, we identify a novel class of truncating IL4R mutations that are dependent on IL-13 signalling and therapeutically targetable with IL-4Rα-blocking antibodies. Finally, using PhasED-seq5, we demonstrate the clinical value of pretreatment and on-treatment ctDNA levels for longitudinally refining cHL risk prediction and for detection of radiographically occult minimal residual disease. Collectively, these results support the utility of noninvasive strategies for genotyping and dynamic monitoring of cHL, as well as capturing molecularly distinct subtypes with diagnostic, prognostic and therapeutic potential.

Bone marrow skeletal stem/progenitor cell defects in dyskeratosis congenita and telomere biology disorders.

Dyskeratosis congenita (DC) is an inherited multisystem disorder, characterized by oral leukoplakia, nail dystrophy, and abnormal skin pigmentation, as well as high rates of bone marrow (BM) failure, solid tumors, and other medical problems such as osteopenia. DC and telomere biology disorders (collectively referred to as TBD here) are caused by germline mutations in telomere biology genes leading to very short telomeres and limited proliferative potential of hematopoietic stem cells. We found that skeletal stem cells (SSCs) within the BM stromal cell population (BMSCs, also known as BM-derived mesenchymal stem cells), may contribute to the hematologic phenotype. TBD-BMSCs exhibited reduced clonogenicity, spontaneous differentiation into adipocytes and fibrotic cells, and increased senescence in vitro. Upon in vivo transplantation into mice, TBD-BMSCs failed to form bone or support hematopoiesis, unlike normal BMSCs. TERC reduction (a TBD-associated gene) in normal BMSCs by small interfering TERC-RNA (siTERC-RNA) recapitulated the TBD-BMSC phenotype by reducing proliferation and secondary colony-forming efficiency, and by accelerating senescence in vitro. Microarray profiles of control and siTERC-BMSCs showed decreased hematopoietic factors at the messenger RNA level and decreased secretion of factors at the protein level. These findings are consistent with defects in SSCs/BMSCs contributing to BM failure in TBD.

A human lymphoma organoid model for evaluating and targeting the follicular lymphoma tumor immune microenvironment.

Heterogeneity in the tumor microenvironment (TME) of follicular lymphomas (FLs) can affect clinical outcomes. Current immunotherapeutic strategies, including antibody- and cell-based therapies, variably overcome pro-tumorigenic mechanisms for sustained disease control. Modeling the intact FL TME, with its native, syngeneic tumor-infiltrating leukocytes, is a major challenge. Here, we describe an organoid culture method for cultivating patient-derived lymphoma organoids (PDLOs), which include cells from the native FL TME. We define the robustness of this method by successfully culturing cryopreserved FL specimens from diverse patients and demonstrate the stability of TME cellular composition, tumor somatic mutations, gene expression profiles, and B/T cell receptor dynamics over 3 weeks. PDLOs treated with CD3:CD19 and CD3:CD20 therapeutic bispecific antibodies showed B cell killing and T cell activation. This stable system offers a robust platform for advancing precision medicine efforts in FL through patient-specific modeling, high-throughput screening, TME signature identification, and treatment response evaluation.

Integrative genomic analysis identifies unique immune environments associated with immunotherapy response in diffuse large B cell lymphoma.

Most diffuse large B-cell lymphoma (DLBCL) patients treated with bispecific antibodies (BsAb) or chimeric antigen receptor (CAR) T cells fail to achieve durable treatment responses, underscoring the need for a deeper understanding of mechanisms that regulate the immune environment and response to treatment. Here, an integrative, multi-omic approach was employed to characterize DLBCL immune environments, which effectively segregated DLBCLs into four quadrants - termed DLBCL-immune quadrants (IQ) - defined by cell-of-origin and immune-related gene set expression scores. Recurrent genomic alterations were enriched in each IQ, suggesting that lymphoma cell-intrinsic alterations contribute to orchestrating unique DLBCL immune environments. In relapsed/refractory DLBCL patients, DLBCL-IQ assignment correlated significantly with clinical benefit with the CD20 x CD3 BsAb, mosunetuzumab, but not with CD19-directed CAR T cells. DLBCL-IQ provides a new framework to conceptualize the DLBCL immune landscape and uncovers the differential impact of the endogenous immune environment on outcomes to BsAb and CAR T cell treatment.

Cell-free DNA in large B-cell lymphoma: MRD and beyond.

Large B-cell lymphomas (LBCLs) are a strikingly diverse set of diseases, including clinical, biological, and molecular heterogeneity. Despite a wealth of information resolving this heterogeneity in the research setting, applying molecular features routinely in the clinic remains challenging. The advent of circulating tumor DNA (ctDNA) liquid biopsies promises to unlock additional molecular information in the clinic, including mutational genotyping, molecular classification, and minimal residual disease detection. Here, we examine the technologies, applications, and studies exploring the utility of ctDNA in LBCLs.

Tracing Founder Mutations in Circulating and Tissue-Resident Follicular Lymphoma Precursors.

Follicular lymphomas (FL) are characterized by BCL2 translocations, often detectable in blood years before FL diagnosis, but also observed in aging healthy individuals, suggesting additional lesions are required for lymphomagenesis. We directly characterized early cooperating mutations by ultradeep sequencing of prediagnostic blood and tissue specimens from 48 subjects who ultimately developed FL. Strikingly, CREBBP lysine acetyltransferase (KAT) domain mutations were the most commonly observed precursor lesions, and largely distinguished patients developing FL (14/48, 29%) from healthy adults with or without detected BCL2 rearrangements (0/13, P = 0.03 and 0/20, P = 0.007, respectively). CREBBP variants were detectable a median of 5.8 years before FL diagnosis, were clonally selected in FL tumors, and appeared restricted to the committed B-cell lineage. These results suggest that mutations affecting the CREBBP KAT domain are common lesions in FL cancer precursor cells (CPC), with the potential for discriminating subjects at risk of developing FL or monitoring residual disease. SIGNIFICANCE: Our study provides direct evidence for recurrent genetic aberrations preceding FL diagnosis, revealing the combination of BCL2 translocation with CREBBP KAT domain mutations as characteristic committed lesions of FL CPCs. Such prediagnostic mutations are detectable years before clinical diagnosis and may help discriminate individuals at risk for lymphoma development. This article is highlighted in the In This Issue feature, p. 1275.

Determinants of resistance to engineered T cell therapies targeting CD19 in large B cell lymphomas.

Most relapsed/refractory large B cell lymphoma (r/rLBCL) patients receiving anti-CD19 chimeric antigen receptor (CAR19) T cells relapse. To characterize determinants of resistance, we profiled over 700 longitudinal specimens from two independent cohorts (n = 65 and n = 73) of r/rLBCL patients treated with axicabtagene ciloleucel. A method for simultaneous profiling of circulating tumor DNA (ctDNA), cell-free CAR19 (cfCAR19) retroviral fragments, and cell-free T cell receptor rearrangements (cfTCR) enabled integration of tumor and both engineered and non-engineered T cell effector-mediated factors for assessing treatment failure and predicting outcomes. Alterations in multiple classes of genes are associated with resistance, including B cell identity (PAX5 and IRF8), immune checkpoints (CD274), and those affecting the microenvironment (TMEM30A). Somatic tumor alterations affect CAR19 therapy at multiple levels, including CAR19 T cell expansion, persistence, and tumor microenvironment. Further, CAR19 T cells play a reciprocal role in shaping tumor genotype and phenotype. We envision these findings will facilitate improved chimeric antigen receptor (CAR) T cells and personalized therapeutic approaches.

Inferring gene expression from cell-free DNA fragmentation profiles.

Profiling of circulating tumor DNA (ctDNA) in the bloodstream shows promise for noninvasive cancer detection. Chromatin fragmentation features have previously been explored to infer gene expression profiles from cell-free DNA (cfDNA), but current fragmentomic methods require high concentrations of tumor-derived DNA and provide limited resolution. Here we describe promoter fragmentation entropy as an epigenomic cfDNA feature that predicts RNA expression levels at individual genes. We developed 'epigenetic expression inference from cell-free DNA-sequencing' (EPIC-seq), a method that uses targeted sequencing of promoters of genes of interest. Profiling 329 blood samples from 201 patients with cancer and 87 healthy adults, we demonstrate classification of subtypes of lung carcinoma and diffuse large B cell lymphoma. Applying EPIC-seq to serial blood samples from patients treated with PD-(L)1 immune-checkpoint inhibitors, we show that gene expression profiles inferred by EPIC-seq are correlated with clinical response. Our results indicate that EPIC-seq could enable noninvasive, high-throughput tissue-of-origin characterization with diagnostic, prognostic and therapeutic potential.

Secreted frizzled related-protein 2 (Sfrp2) deficiency decreases adult skeletal stem cell function in mice.

In a previous transcriptomic study of human bone marrow stromal cells (BMSCs, also known as bone marrow-derived "mesenchymal stem cells"), SFRP2 was highly over-represented in a subset of multipotent BMSCs (skeletal stem cells, SSCs), which recreate a bone/marrow organ in an in vivo ectopic bone formation assay. SFRPs modulate WNT signaling, which is essential to maintain skeletal homeostasis, but the specific role of SFRP2 in BMSCs/SSCs is unclear. Here, we evaluated Sfrp2 deficiency on BMSC/SSC function in models of skeletal organogenesis and regeneration. The skeleton of Sfrp2-deficient (KO) mice is overtly normal; but their BMSCs/SSCs exhibit reduced colony-forming efficiency, reflecting low SSC self-renewal/abundancy. Sfrp2 KO BMSCs/SSCs formed less trabecular bone than those from WT littermates in the ectopic bone formation assay. Moreover, regeneration of a cortical drilled hole defect was dramatically impaired in Sfrp2 KO mice. Sfrp2-deficient BMSCs/SSCs exhibited poor in vitro osteogenic differentiation as measured by Runx2 and Osterix expression and calcium accumulation. Interestingly, activation of the Wnt co-receptor, Lrp6, and expression of Wnt target genes, Axin2, C-myc and Cyclin D1, were reduced in Sfrp2-deficient BMSCs/SSCs. Addition of recombinant Sfrp2 restored most of these activities, suggesting that Sfrp2 acts as a Wnt agonist. We demonstrate that Sfrp2 plays a role in self-renewal of SSCs and in the recruitment and differentiation of adult SSCs during bone healing. SFRP2 is also a useful marker of BMSC/SSC multipotency, and a factor to potentially improve the quality of ex vivo expanded BMSC/SSC products.

The landscape of tumor cell states and ecosystems in diffuse large B cell lymphoma.

Biological heterogeneity in diffuse large B cell lymphoma (DLBCL) is partly driven by cell-of-origin subtypes and associated genomic lesions, but also by diverse cell types and cell states in the tumor microenvironment (TME). However, dissecting these cell states and their clinical relevance at scale remains challenging. Here, we implemented EcoTyper, a machine-learning framework integrating transcriptome deconvolution and single-cell RNA sequencing, to characterize clinically relevant DLBCL cell states and ecosystems. Using this approach, we identified five cell states of malignant B cells that vary in prognostic associations and differentiation status. We also identified striking variation in cell states for 12 other lineages comprising the TME and forming cell state interactions in stereotyped ecosystems. While cell-of-origin subtypes have distinct TME composition, DLBCL ecosystems capture clinical heterogeneity within existing subtypes and extend beyond cell-of-origin and genotypic classes. These results resolve the DLBCL microenvironment at systems-level resolution and identify opportunities for therapeutic targeting (https://ecotyper.stanford.edu/lymphoma).

Short Diagnosis-to-Treatment Interval Is Associated With Higher Circulating Tumor DNA Levels in Diffuse Large B-Cell Lymphoma.

PURPOSE: Patients with Diffuse Large B-cell Lymphoma (DLBCL) in need of immediate therapy are largely under-represented in clinical trials. The diagnosis-to-treatment interval (DTI) has recently been described as a metric to quantify such patient selection bias, with short DTI being associated with adverse risk factors and inferior outcomes. Here, we characterized the relationships between DTI, circulating tumor DNA (ctDNA), conventional risk factors, and clinical outcomes, with the goal of defining objective disease metrics contributing to selection bias. PATIENTS AND METHODS: We evaluated pretreatment ctDNA levels in 267 patients with DLBCL treated across multiple centers in Europe and the United States using Cancer Personalized Profiling by Deep Sequencing. Pretreatment ctDNA levels were correlated with DTI, total metabolic tumor volumes (TMTVs), the International Prognostic Index (IPI), and outcome. RESULTS: Short DTI was associated with advanced-stage disease (P < .001) and higher IPI (P < .001). We also found an inverse correlation between DTI and TMTV (RS = -0.37; P < .001). Similarly, pretreatment ctDNA levels were significantly associated with stage, IPI, and TMTV (all P < .001), demonstrating that both DTI and ctDNA reflect disease burden. Notably, patients with shorter DTI had higher pretreatment ctDNA levels (P < .001). Pretreatment ctDNA levels predicted short DTI independent of the IPI (P < .001). Although each risk factor was significantly associated with event-free survival in univariable analysis, ctDNA level was prognostic of event-free survival independent of DTI and IPI in multivariable Cox regression (ctDNA: hazard ratio, 1.5; 95% CI [1.2 to 2.0]; IPI: 1.1 [0.9 to 1.3]; -DTI: 1.1 [1.0 to 1.2]). CONCLUSION: Short DTI largely reflects baseline tumor burden, which can be objectively measured using pretreatment ctDNA levels. Pretreatment ctDNA levels therefore have utility for quantifying and guarding against selection biases in prospective DLBCL clinical trials.

Enhanced detection of minimal residual disease by targeted sequencing of phased variants in circulating tumor DNA.

Circulating tumor-derived DNA (ctDNA) is an emerging biomarker for many cancers, but the limited sensitivity of current detection methods reduces its utility for diagnosing minimal residual disease. Here we describe phased variant enrichment and detection sequencing (PhasED-seq), a method that uses multiple somatic mutations in individual DNA fragments to improve the sensitivity of ctDNA detection. Leveraging whole-genome sequences from 2,538 tumors, we identify phased variants and their associations with mutational signatures. We show that even without molecular barcodes, the limits of detection of PhasED-seq outperform prior methods, including duplex barcoding, allowing ctDNA detection in the ppm range in participant samples. We profiled 678 specimens from 213 participants with B cell lymphomas, including serial cell-free DNA samples before and during therapy for diffuse large B cell lymphoma. In participants with undetectable ctDNA after two cycles of therapy using a next-generation sequencing-based approach termed cancer personalized profiling by deep sequencing, an additional 25% have ctDNA detectable by PhasED-seq and have worse outcomes. Finally, we demonstrate the application of PhasED-seq to solid tumors.

Predicting HLA class II antigen presentation through integrated deep learning.

Accurate prediction of antigen presentation by human leukocyte antigen (HLA) class II molecules would be valuable for vaccine development and cancer immunotherapies. Current computational methods trained on in vitro binding data are limited by insufficient training data and algorithmic constraints. Here we describe MARIA (major histocompatibility complex analysis with recurrent integrated architecture; https://maria.stanford.edu/ ), a multimodal recurrent neural network for predicting the likelihood of antigen presentation from a gene of interest in the context of specific HLA class II alleles. In addition to in vitro binding measurements, MARIA is trained on peptide HLA ligand sequences identified by mass spectrometry, expression levels of antigen genes and protease cleavage signatures. Because it leverages these diverse training data and our improved machine learning framework, MARIA (area under the curve = 0.89-0.92) outperformed existing methods in validation datasets. Across independent cancer neoantigen studies, peptides with high MARIA scores are more likely to elicit strong CD4+ T cell responses. MARIA allows identification of immunogenic epitopes in diverse cancers and autoimmune disease.

Molecular profile of clonal strains of human skeletal stem/progenitor cells with different potencies.

Bone marrow stromal cells (BMSCs, also known as bone marrow-derived mesenchymal stem cells) are fibroblastic reticular cells, a subset of which is composed of multipotent skeletal stem cells (SSCs). SSCs/BMSCs are able to recreate a bone/marrow organ in vivo. To determine differences between clonogenic multipotent SSCs and similarly clonogenic but non-multipotent BMSCs, we established single colony-derived strains (SCDSs, initiated by individual Colony Forming Unit-Fibroblasts) and determined their differentiation capacity by vivo transplantation. In this series of human SCDSs (N=24), 20.8% formed fibrous tissue (F), 66.7% formed bone (B), and 12.5% formed a bone/marrow organ, and thus were multipotent (M). RNA isolated from 12 SCDSs just prior to transplantation was analyzed by microarray. Although highly similar, there was variability from one SCDS to another, and SCDSs did not strictly segregate into the three functional groups (F, B or M) by unsupervised hierarchical clustering. We then compared 3 F-SCDSs to 3 M-SCDSs that did segregate. Genes associated with skeletogenesis, osteoblastogeneis, hematopoiesis, and extracellular matrix were over-represented in M-SCDSs compared with F-SCDSs. These results highlight the heterogeneity of SSCs/BMSCs, even between functionally similar SCDSs, but also indicate that differences can be detected that may shed light on the character of the SSC.