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.

CAR19 monitoring by peripheral blood immunophenotyping reveals histology-specific expansion and toxicity.

ABSTRACT: Chimeric antigen receptor (CAR) T cells directed against CD19 (CAR19) are a revolutionary treatment for B-cell lymphomas (BCLs). CAR19 cell expansion is necessary for CAR19 function but is also associated with toxicity. To define the impact of CAR19 expansion on patient outcomes, we prospectively followed a cohort of 236 patients treated with CAR19 (brexucabtagene autoleucel or axicabtagene ciloleucel) for mantle cell lymphoma (MCL), follicular lymphoma, and large BCL (LBCL) over the course of 5 years and obtained CAR19 expansion data using peripheral blood immunophenotyping for 188 of these patients. CAR19 expansion was higher in patients with MCL than other lymphoma histologic subtypes. Notably, patients with MCL had increased toxicity and required fourfold higher cumulative steroid doses than patients with LBCL. CAR19 expansion was associated with the development of cytokine release syndrome, immune effector cell-associated neurotoxicity syndrome, and the requirement for granulocyte colony-stimulating factor 14 days after infusion. Younger patients and those with elevated lactate dehydrogenase (LDH) had significantly higher CAR19 expansion. In general, no association between CAR19 expansion and LBCL treatment response was observed. However, when controlling for tumor burden, we found that lower CAR19 expansion in conjunction with low LDH was associated with improved outcomes in LBCL. In sum, this study finds CAR19 expansion principally associates with CAR-related toxicity. Additionally, CAR19 expansion as measured by peripheral blood immunophenotyping may be dispensable to favorable outcomes in LBCL.

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.

Clinical presentation and treatment responses in IgM-related AL amyloidosis.

Amyloid light-chain (AL) amyloidosis is a multi-organ disease due to deposition of misfolded monoclonal immunoglobulin light chains. IgM AL amyloidosis is a rare variant, about 6% of AL amyloidosis cases, and more data are needed for treatment guidance. In IgM AL amyloidosis, the clonal cell of origin may be a plasma or lymphoplasmacytic cell, and treatments targeting each are employed. We describe presenting clinical and laboratory features of 95 patients with IgM AL amyloidosis treated at Boston University Amyloidosis Center from 1996 to 2012. The median diagnosis age was 66 years (range: 38-89) with 56% males. Organ involvement rates were: kidney (51%); heart (40%); lymph nodes (25%) and gastrointestinal tract (17%). Treatment responses were analyzed for 46 patients seen after 2003. Five treatment regimens were assigned by bone marrow pathology and patient-specific factors. Overall hematologic response rates and very good partial or complete hematologic response rates, respectively, were: high-dose melphalan/stem cell transplant (HDM/SCT) 100%;80%, Bortezomib 82%;27%, Rituximab 80%;27%, immunomodulatory agents (IMids) 75%;0%, and standard dose alkylating agents (Melphalan or cyclophosphamide) 63%;19%. Overall, 5-year survival rates were significantly higher in patients with a hematological response: 79.2 ± 8.5% versus 41 ± 14.9% in non-responders, which is more favorable than typically expected in AL amyloidosis.

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.

Stromal-derived IL-6 alters the balance of myeloerythroid progenitors during Toxoplasma gondii infection.

Inflammation alters hematopoiesis, often by decreasing erythropoiesis and enhancing myeloid output. The mechanisms behind these changes and how the BM stroma contributes to this process are active areas of research. In this study, we examine these questions in the setting of murine Toxoplasma gondii infection. Our data reveal that infection alters early myeloerythroid differentiation, blocking erythroid development beyond the Pre MegE stage, while expanding the GMP population. IL-6 was found to be a critical mediator of these differences, independent of hepcidin-induced iron restriction. Comparing the BM with the spleen showed that the hematopoietic response was driven by the local microenvironment, and BM chimeras demonstrated that radioresistant cells were the relevant source of IL-6 in vivo. Finally, direct ex vivo sorting revealed that VCAM(+)CD146(lo) BM stromal fibroblasts significantly increase IL-6 secretion after infection. These data suggest that BMSCs regulate the hematopoietic changes during inflammation via IL-6.

Superparamagnetic iron oxide nanoparticles labeling of bone marrow stromal (mesenchymal) cells does not affect their "stemness".

Superparamagnetic iron oxide nanoparticles (SPION) are increasingly used to label human bone marrow stromal cells (BMSCs, also called "mesenchymal stem cells") to monitor their fate by in vivo MRI, and by histology after Prussian blue (PB) staining. SPION-labeling appears to be safe as assessed by in vitro differentiation of BMSCs, however, we chose to resolve the question of the effect of labeling on maintaining the "stemness" of cells within the BMSC population in vivo. Assays performed include colony forming efficiency, CD146 expression, gene expression profiling, and the "gold standard" of evaluating bone and myelosupportive stroma formation in vivo in immuncompromised recipients. SPION-labeling did not alter these assays. Comparable abundant bone with adjoining host hematopoietic cells were seen in cohorts of mice that were implanted with SPION-labeled or unlabeled BMSCs. PB+ adipocytes were noted, demonstrating their donor origin, as well as PB+ pericytes, indicative of self-renewal of the stem cell in the BMSC population. This study confirms that SPION labeling does not alter the differentiation potential of the subset of stem cells within BMSCs.