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.

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.

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.

Evaluating upfront high-dose consolidation after R-CHOP for follicular lymphoma by clinical and genetic risk models.

High-dose therapy and autologous stem cell transplantation (HDT/ASCT) is an effective salvage treatment for eligible patients with follicular lymphoma (FL) and early progression of disease (POD). Since the introduction of rituximab, HDT/ASCT is no longer recommended in first remission. We here explored whether consolidative HDT/ASCT improved survival in defined subgroups of previously untreated patients. We report survival analyses of 431 patients who received frontline rituximab, cyclophosphamide, doxorubicin, vincristine, and prednisone (R-CHOP) for advanced FL, and were randomized to receive consolidative HDT/ASCT. We performed targeted genotyping of 157 diagnostic biopsies, and calculated genotype-based risk scores. HDT/ASCT improved failure-free survival (FFS; hazard ratio [HR], 0.8, P = .07; as-treated: HR, 0.7, P = .04), but not overall survival (OS; HR, 1.3, P = .27; as-treated: HR, 1.4, P = .13). High-risk cohorts identified by FL International Prognostic Index (FLIPI), and the clinicogenetic risk models m7-FLIPI and POD within 24 months-prognostic index (POD24-PI) comprised 27%, 18%, and 22% of patients. HDT/ASCT did not significantly prolong FFS in high-risk patients as defined by FLIPI (HR, 0.9; P = .56), m7-FLIPI (HR, 0.9; P = .91), and POD24-PI (HR, 0.8; P = .60). Similarly, OS was not significantly improved. Finally, we used a machine-learning approach to predict benefit from HDT/ASCT by genotypes. Patients predicted to benefit from HDT/ASCT had longer FFS with HDT/ASCT (HR, 0.4; P = .03), but OS did not reach statistical significance. Thus, consolidative HDT/ASCT after frontline R-CHOP did not improve OS in unselected FL patients and subgroups selected by genotype-based risk models.

An experimental design framework for Markovian gene regulatory networks under stationary control policy.

BACKGROUND: A fundamental problem for translational genomics is to find optimal therapies based on gene regulatory intervention. Dynamic intervention involves a control policy that optimally reduces a cost function based on phenotype by externally altering the state of the network over time. When a gene regulatory network (GRN) model is fully known, the problem is addressed using classical dynamic programming based on the Markov chain associated with the network. When the network is uncertain, a Bayesian framework can be applied, where policy optimality is with respect to both the dynamical objective and the uncertainty, as characterized by a prior distribution. In the presence of uncertainty, it is of great practical interest to develop an experimental design strategy and thereby select experiments that optimally reduce a measure of uncertainty. RESULTS: In this paper, we employ mean objective cost of uncertainty (MOCU), which quantifies uncertainty based on the degree to which uncertainty degrades the operational objective, that being the cost owing to undesirable phenotypes. We assume that a number of conditional probabilities characterizing regulatory relationships among genes are unknown in the Markovian GRN. In sum, there is a prior distribution which can be updated to a posterior distribution by observing a regulatory trajectory, and an optimal control policy, known as an "intrinsically Bayesian robust" (IBR) policy. To obtain a better IBR policy, we select an experiment that minimizes the MOCU remaining after applying its output to the network. At this point, we can either stop and find the resulting IBR policy or proceed to determine more unknown conditional probabilities via regulatory observation and find the IBR policy from the resulting posterior distribution. For sequential experimental design this entire process is iterated. Owing to the computational complexity of experimental design, which requires computation of many potential IBR policies, we implement an approximate method utilizing mean first passage times (MFPTs) - but only in experimental design, the final policy being an IBR policy. CONCLUSIONS: Comprehensive performance analysis based on extensive simulations on synthetic and real GRNs demonstrate the efficacy of the proposed method, including the accuracy and computational advantage of the approximate MFPT-based design.

Effect of separate sampling on classification accuracy.

MOTIVATION: Measurements are commonly taken from two phenotypes to build a classifier, where the number of data points from each class is predetermined, not random. In this 'separate sampling' scenario, the data cannot be used to estimate the class prior probabilities. Moreover, predetermined class sizes can severely degrade classifier performance, even for large samples. RESULTS: We employ simulations using both synthetic and real data to show the detrimental effect of separate sampling on a variety of classification rules. We establish propositions related to the effect on the expected classifier error owing to a sampling ratio different from the population class ratio. From these we derive a sample-based minimax sampling ratio and provide an algorithm for approximating it from the data. We also extend to arbitrary distributions the classical population-based Anderson linear discriminant analysis minimax sampling ratio derived from the discriminant form of the Bayes classifier. AVAILABILITY: All the codes for synthetic data and real data examples are written in MATLAB. A function called mmratio, whose output is an approximation of the minimax sampling ratio of a given dataset, is also written in MATLAB. All the codes are available at: http://gsp.tamu.edu/Publications/supplementary/shahrokh13b.