Crowd-sourced benchmarking of single-sample tumor subclonal reconstruction.

Subclonal reconstruction algorithms use bulk DNA sequencing data to quantify parameters of tumor evolution, allowing an assessment of how cancers initiate, progress and respond to selective pressures. We launched the ICGC-TCGA (International Cancer Genome Consortium-The Cancer Genome Atlas) DREAM Somatic Mutation Calling Tumor Heterogeneity and Evolution Challenge to benchmark existing subclonal reconstruction algorithms. This 7-year community effort used cloud computing to benchmark 31 subclonal reconstruction algorithms on 51 simulated tumors. Algorithms were scored on seven independent tasks, leading to 12,061 total runs. Algorithm choice influenced performance substantially more than tumor features but purity-adjusted read depth, copy-number state and read mappability were associated with the performance of most algorithms on most tasks. No single algorithm was a top performer for all seven tasks and existing ensemble strategies were unable to outperform the best individual methods, highlighting a key research need. All containerized methods, evaluation code and datasets are available to support further assessment of the determinants of subclonal reconstruction accuracy and development of improved methods to understand tumor evolution.

Allosteric PI3Kα Inhibition Overcomes On-target Resistance to Orthosteric Inhibitors Mediated by Secondary PIK3CA Mutations.

PIK3CA mutations occur in ∼8% of cancers, including ∼40% of HR-positive breast cancers, where the PI3K-alpha (PI3Kα)-selective inhibitor alpelisib is FDA approved in combination with fulvestrant. Although prior studies have identified resistance mechanisms, such as PTEN loss, clinically acquired resistance to PI3Kα inhibitors remains poorly understood. Through serial liquid biopsies and rapid autopsies in 39 patients with advanced breast cancer developing acquired resistance to PI3Kα inhibitors, we observe that 50% of patients acquire genomic alterations within the PI3K pathway, including PTEN loss and activating AKT1 mutations. Notably, although secondary PIK3CA mutations were previously reported to increase sensitivity to PI3Kα inhibitors, we identified emergent secondary resistance mutations in PIK3CA that alter the inhibitor binding pocket. Some mutations had differential effects on PI3Kα-selective versus pan-PI3K inhibitors, but resistance induced by all mutations could be overcome by the novel allosteric pan-mutant-selective PI3Kα-inhibitor RLY-2608. Together, these findings provide insights to guide strategies to overcome resistance in PIK3CA-mutated cancers. SIGNIFICANCE: In one of the largest patient cohorts analyzed to date, this study defines the clinical landscape of acquired resistance to PI3Kα inhibitors. Genomic alterations within the PI3K pathway represent a major mode of resistance and identify a novel class of secondary PIK3CA resistance mutations that can be overcome by an allosteric PI3Kα inhibitor. See related commentary by Gong and Vanhaesebroeck, p. 204 . See related article by Varkaris et al., p. 240 . This article is featured in Selected Articles from This Issue, p. 201.

Germline mutations and developmental mosaicism underlying EGFR-mutant lung cancer.

While the development of multiple primary tumors in smokers with lung cancer can be attributed to carcinogen-induced field cancerization, the occurrence of multiple primary tumors in individuals with EGFR-mutant lung cancer who lack known environmental exposures remains unexplained. We identified ten patients with early-stage, resectable non-small cell lung cancer who presented with multiple anatomically distinct EGFR-mutant tumors. We analyzed the phylogenetic relationships among multiple tumors from each patient using whole exome sequencing (WES) and hypermutable poly-guanine (poly-G) repeat genotyping, as orthogonal methods for lineage tracing. In two patients, we identified germline EGFR variants, which confer moderately enhanced signaling when modeled in vitro. In four other patients, developmental mosaicism is supported by the poly-G lineage tracing and WES, indicating a common non-germline cell-of-origin. Thus, developmental mosaicism and germline variants define two distinct mechanisms of genetic predisposition to multiple EGFR-mutant primary tumors, with implications for understanding their etiology and clinical management.

Genetic events associated with venetoclax resistance in CLL identified by whole-exome sequencing of patient samples.

Although BCL2 mutations are reported as later occurring events leading to venetoclax resistance, many other mechanisms of progression have been reported though remain poorly understood. Here, we analyze longitudinal tumor samples from 11 patients with disease progression while receiving venetoclax to characterize the clonal evolution of resistance. All patients tested showed increased in vitro resistance to venetoclax at the posttreatment time point. We found the previously described acquired BCL2-G101V mutation in only 4 of 11 patients, with 2 patients showing a very low variant allele fraction (0.03%-4.68%). Whole-exome sequencing revealed acquired loss(8p) in 4 of 11 patients, of which 2 patients also had gain (1q21.2-21.3) in the same cells affecting the MCL1 gene. In vitro experiments showed that CLL cells from the 4 patients with loss(8p) were more resistant to venetoclax than cells from those without it, with the cells from 2 patients also carrying gain (1q21.2-21.3) showing increased sensitivity to MCL1 inhibition. Progression samples with gain (1q21.2-21.3) were more susceptible to the combination of MCL1 inhibitor and venetoclax. Differential gene expression analysis comparing bulk RNA sequencing data from pretreatment and progression time points of all patients showed upregulation of proliferation, B-cell receptor (BCR), and NF-κB gene sets including MAPK genes. Cells from progression time points demonstrated upregulation of surface immunoglobulin M and higher pERK levels compared with those from the preprogression time point, suggesting an upregulation of BCR signaling that activates the MAPK pathway. Overall, our data suggest several mechanisms of acquired resistance to venetoclax in CLL that could pave the way for rationally designed combination treatments for patients with venetoclax-resistant CLL.

Genomic and transcriptomic analysis of checkpoint blockade response in advanced non-small cell lung cancer.

Anti-PD-1/PD-L1 agents have transformed the treatment landscape of advanced non-small cell lung cancer (NSCLC). To expand our understanding of the molecular features underlying response to checkpoint inhibitors in NSCLC, we describe here the first joint analysis of the Stand Up To Cancer-Mark Foundation cohort, a resource of whole exome and/or RNA sequencing from 393 patients with NSCLC treated with anti-PD-(L)1 therapy, along with matched clinical response annotation. We identify a number of associations between molecular features and outcome, including (1) favorable (for example, ATM altered) and unfavorable (for example, TERT amplified) genomic subgroups, (2) a prominent association between expression of inducible components of the immunoproteasome and response and (3) a dedifferentiated tumor-intrinsic subtype with enhanced response to checkpoint blockade. Taken together, results from this cohort demonstrate the complexity of biological determinants underlying immunotherapy outcomes and reinforce the discovery potential of integrative analysis within large, well-curated, cancer-specific cohorts.

Inferring early genetic progression in cancers with unobtainable premalignant disease.

Analysis of premalignant tissue has identified the typical order of somatic events leading to invasive tumors in several cancer types. For other cancers, premalignant tissue is unobtainable, leaving genetic progression unknown. Here, we demonstrate how to infer progression from exome sequencing of primary tumors. Our computational method, PhylogicNDT, recapitulated the previous experimentally determined genetic progression of human papillomavirus-negative (HPV-) head and neck squamous cell carcinoma (HNSCC). We then evaluated HPV+ HNSCC, which lacks premalignant tissue, and uncovered its previously unknown progression, identifying early drivers. We converted relative timing estimates of driver mutations and HPV integration to years before diagnosis based on a clock-like mutational signature. We associated the timing of transitions to aneuploidy with increased intratumor genetic heterogeneity and shorter overall survival. Our approach can establish previously unknown early genetic progression of cancers with unobtainable premalignant tissue, supporting development of experimental models and methods for early detection, interception and prognostication.

Evolutionary history of transformation from chronic lymphocytic leukemia to Richter syndrome.

Richter syndrome (RS) arising from chronic lymphocytic leukemia (CLL) exemplifies an aggressive malignancy that develops from an indolent neoplasm. To decipher the genetics underlying this transformation, we computationally deconvoluted admixtures of CLL and RS cells from 52 patients with RS, evaluating paired CLL-RS whole-exome sequencing data. We discovered RS-specific somatic driver mutations (including IRF2BP2, SRSF1, B2M, DNMT3A and CCND3), recurrent copy-number alterations beyond del(9p21)(CDKN2A/B), whole-genome duplication and chromothripsis, which were confirmed in 45 independent RS cases and in an external set of RS whole genomes. Through unsupervised clustering, clonally related RS was largely distinct from diffuse large B cell lymphoma. We distinguished pathways that were dysregulated in RS versus CLL, and detected clonal evolution of transformation at single-cell resolution, identifying intermediate cell states. Our study defines distinct molecular subtypes of RS and highlights cell-free DNA analysis as a potential tool for early diagnosis and monitoring.

Pirtobrutinib targets BTK C481S in ibrutinib-resistant CLL but second-site BTK mutations lead to resistance.

Covalent inhibitors of Bruton tyrosine kinase (BTK) have transformed the therapy of chronic lymphocytic leukemia (CLL), but continuous therapy has been complicated by the development of resistance. The most common resistance mechanism in patients whose disease progresses on covalent BTK inhibitors (BTKis) is a mutation in the BTK 481 cysteine residue to which the inhibitors bind covalently. Pirtobrutinib is a highly selective, noncovalent BTKi with substantial clinical activity in patients whose disease has progressed on covalent BTKi, regardless of BTK mutation status. Using in vitro ibrutinib-resistant models and cells from patients with CLL, we show that pirtobrutinib potently inhibits BTK-mediated functions including B-cell receptor (BCR) signaling, cell viability, and CCL3/CCL4 chemokine production in both BTK wild-type and C481S mutant CLL cells. We demonstrate that primary CLL cells from responding patients on the pirtobrutinib trial show reduced BCR signaling, cell survival, and CCL3/CCL4 chemokine secretion. At time of progression, these primary CLL cells show increasing resistance to pirtobrutinib in signaling inhibition, cell viability, and cytokine production. We employed longitudinal whole-exome sequencing on 2 patients whose disease progressed on pirtobrutinib and identified selection of alternative-site BTK mutations, providing clinical evidence that secondary BTK mutations lead to resistance to noncovalent BTKis.

osr1 Maintains Renal Progenitors and Regulates Podocyte Development by Promoting wnt2ba via the Antagonism of hand2.

Knowledge about the genetic pathways that control nephron development is essential for better understanding the basis of congenital malformations of the kidney. The transcription factors Osr1 and Hand2 are known to exert antagonistic influences to balance kidney specification. Here, we performed a forward genetic screen to identify nephrogenesis regulators, where whole genome sequencing identified an osr1 lesion in the novel oceanside (ocn) mutant. The characterization of the mutant revealed that osr1 is needed to specify not renal progenitors but rather their maintenance. Additionally, osr1 promotes the expression of wnt2ba in the intermediate mesoderm (IM) and later the podocyte lineage. wnt2ba deficiency reduced podocytes, where overexpression of wnt2ba was sufficient to rescue podocytes and osr1 deficiency. Antagonism between osr1 and hand2 mediates podocyte development specifically by controlling wnt2ba expression. These studies reveal new insights about the roles of Osr1 in promoting renal progenitor survival and lineage choice.

Clinical Efficacy and Whole-Exome Sequencing of Liquid Biopsies in a Phase IB/II Study of Bazedoxifene and Palbociclib in Advanced Hormone Receptor-Positive Breast Cancer.

PURPOSE: Sensitivity to endocrine therapy (ET) is critical for the clinical benefit from the combination of palbociclib plus ET in hormone receptor-positive/HER2-negative (HR+/HER2-) advanced breast cancer. Bazedoxifene is a third-generation selective estrogen receptor (ER) modulator and selective ER degrader with activity in preclinical models of endocrine-resistant breast cancer, including models harboring ESR1 mutations. Clinical trials in healthy women showed that bazedoxifene is well tolerated. PATIENTS AND METHODS: We conducted a phase Ib/II study of bazedoxifene plus palbociclib in patients with HR+/HER2- advanced breast cancer who progressed on prior ET (N = 36; NCT02448771). RESULTS: The study met its primary endpoint, with a clinical benefit rate of 33.3%, and the safety profile was consistent with what has previously been seen with palbociclib monotherapy. The median progression-free survival (PFS) was 3.6 months [95% confidence interval (CI), 2.0-7.2]. An activating PIK3CA mutation at baseline was associated with a shorter PFS (HR = 4.4; 95% CI, 1.5-13; P = 0.0026), but activating ESR1 mutations did not impact the PFS. Longitudinal plasma circulating tumor DNA whole-exome sequencing (WES; N = 68 plasma samples) provided an overview of the tumor heterogeneity and the subclonal genetic evolution, and identified actionable mutations acquired during treatment. CONCLUSIONS: The combination of palbociclib and bazedoxifene has clinical efficacy and an acceptable safety profile in a heavily pretreated patient population with advanced HR+/HER2- breast cancer. These results merit continued investigation of bazedoxifene in breast cancer.

Combined tumor and immune signals from genomes or transcriptomes predict outcomes of checkpoint inhibition in melanoma.

Immune checkpoint blockade (CPB) improves melanoma outcomes, but many patients still do not respond. Tumor mutational burden (TMB) and tumor-infiltrating T cells are associated with response, and integrative models improve survival prediction. However, integrating immune/tumor-intrinsic features using data from a single assay (DNA/RNA) remains underexplored. Here, we analyze whole-exome and bulk RNA sequencing of tumors from new and published cohorts of 189 and 178 patients with melanoma receiving CPB, respectively. Using DNA, we calculate T cell and B cell burdens (TCB/BCB) from rearranged TCR/Ig sequences and find that patients with TMBhigh and TCBhigh or BCBhigh have improved outcomes compared to other patients. By combining pairs of immune- and tumor-expressed genes, we identify three gene pairs associated with response and survival, which validate in independent cohorts. The top model includes lymphocyte-expressed MAP4K1 and tumor-expressed TBX3. Overall, RNA or DNA-based models combining immune and tumor measures improve predictions of melanoma CPB outcomes.

The RNA helicase Ddx21 controls Vegfc-driven developmental lymphangiogenesis by balancing endothelial cell ribosome biogenesis and p53 function.

The development of a functional vasculature requires the coordinated control of cell fate, lineage differentiation and network growth. Cellular proliferation is spatiotemporally regulated in developing vessels, but how this is orchestrated in different lineages is unknown. Here, using a zebrafish genetic screen for lymphatic-deficient mutants, we uncover a mutant for the RNA helicase Ddx21. Ddx21 cell-autonomously regulates lymphatic vessel development. An established regulator of ribosomal RNA synthesis and ribosome biogenesis, Ddx21 is enriched in sprouting venous endothelial cells in response to Vegfc-Flt4 signalling. Ddx21 function is essential for Vegfc-Flt4-driven endothelial cell proliferation. In the absence of Ddx21, endothelial cells show reduced ribosome biogenesis, p53 and p21 upregulation and cell cycle arrest that blocks lymphangiogenesis. Thus, Ddx21 coordinates the lymphatic endothelial cell response to Vegfc-Flt4 signalling by balancing ribosome biogenesis and p53 function. This mechanism may be targetable in diseases of excessive lymphangiogenesis such as cancer metastasis or lymphatic malformation.

Molecular features of exceptional response to neoadjuvant anti-androgen therapy in high-risk localized prostate cancer.

High-risk localized prostate cancer (HRLPC) is associated with a substantial risk of recurrence and disease mortality. Recent clinical trials have shown that intensifying anti-androgen therapies administered before prostatectomy can induce pathologic complete responses or minimal residual disease, called exceptional response, although the molecular determinants of these clinical outcomes are largely unknown. Here, we perform whole-exome and transcriptome sequencing on pre-treatment multi-regional tumor biopsies from exceptional responders (ERs) and non-responders (NRs, pathologic T3 or lymph node-positive disease) to intensive neoadjuvant anti-androgen therapies. Clonal SPOP mutation and SPOPL copy-number loss are exclusively observed in ERs, while clonal TP53 mutation and PTEN copy-number loss are exclusively observed in NRs. Transcriptional programs involving androgen signaling and TGF-ß signaling are enriched in ERs and NRs, respectively. These findings may guide prospective validation studies of these molecular features in large HRLPC clinical cohorts treated with neoadjuvant anti-androgens to improve patient stratification.

Parallel Genomic Alterations of Antigen and Payload Targets Mediate Polyclonal Acquired Clinical Resistance to Sacituzumab Govitecan in Triple-Negative Breast Cancer.

Sacituzumab govitecan (SG), the first antibody-drug conjugate (ADC) approved for triple-negative breast cancer, incorporates the anti-TROP2 antibody hRS7 conjugated to a topoisomerase-1 (TOP1) inhibitor payload. We sought to identify mechanisms of SG resistance through RNA and whole-exome sequencing of pretreatment and postprogression specimens. One patient exhibiting de novo progression lacked TROP2 expression, in contrast to robust TROP2 expression and focal genomic amplification of TACSTD2/TROP2 observed in a patient with a deep, prolonged response to SG. Analysis of acquired genomic resistance in this case revealed one phylogenetic branch harboring a canonical TOP1 E418K resistance mutation and subsequent frameshift TOP1 mutation, whereas a distinct branch exhibited a novel TACSTD2/TROP2 T256R missense mutation. Reconstitution experiments demonstrated that TROP2T256R confers SG resistance via defective plasma membrane localization and reduced cell-surface binding by hRS7. These findings highlight parallel genomic alterations in both antibody and payload targets associated with resistance to SG. SIGNIFICANCE: These findings underscore TROP2 as a response determinant and reveal acquired SG resistance mechanisms involving the direct antibody and drug payload targets in distinct metastatic subclones of an individual patient. This study highlights the specificity of SG and illustrates how such mechanisms will inform therapeutic strategies to overcome ADC resistance.This article is highlighted in the In This Issue feature, p. 2355.

Longitudinal Single-Cell Dynamics of Chromatin Accessibility and Mitochondrial Mutations in Chronic Lymphocytic Leukemia Mirror Disease History.

While cancers evolve during disease progression and in response to therapy, temporal dynamics remain difficult to study in humans due to the lack of consistent barcodes marking individual clones in vivo. We employ mitochondrial single-cell assay for transposase-accessible chromatin with sequencing to profile 163,279 cells from 9 patients with chronic lymphocytic leukemia (CLL) collected across disease course and utilize mitochondrial DNA (mtDNA) mutations as natural genetic markers of cancer clones. We observe stable propagation of mtDNA mutations over years in the absence of strong selective pressure, indicating clonal persistence, but dramatic changes following tight bottlenecks, including disease transformation and relapse posttherapy, paralleled by acquisition of copy-number variants and changes in chromatin accessibility and gene expression. Furthermore, we link CLL subclones to distinct chromatin states, providing insight into nongenetic sources of relapse. mtDNA mutations thus mirror disease history and provide naturally occurring genetic barcodes to enable patient-specific study of cancer subclonal dynamics. SIGNIFICANCE: Single-cell multi-omic profiling of CLL reveals the utility of somatic mtDNA mutations as in vivo barcodes, which mark subclones that can evolve over time along with changes in accessible chromatin and gene expression profiles to capture dynamics of disease evolution. See related commentary by Hilton and Scott, p. 2965. This article is highlighted in the In This Issue feature, p. 2945.

Characterizing genetic intra-tumor heterogeneity across 2,658 human cancer genomes.

Intra-tumor heterogeneity (ITH) is a mechanism of therapeutic resistance and therefore an important clinical challenge. However, the extent, origin, and drivers of ITH across cancer types are poorly understood. To address this, we extensively characterize ITH across whole-genome sequences of 2,658 cancer samples spanning 38 cancer types. Nearly all informative samples (95.1%) contain evidence of distinct subclonal expansions with frequent branching relationships between subclones. We observe positive selection of subclonal driver mutations across most cancer types and identify cancer type-specific subclonal patterns of driver gene mutations, fusions, structural variants, and copy number alterations as well as dynamic changes in mutational processes between subclonal expansions. Our results underline the importance of ITH and its drivers in tumor evolution and provide a pan-cancer resource of comprehensively annotated subclonal events from whole-genome sequencing data.

Activation of the MAPK pathway mediates resistance to PI3K inhibitors in chronic lymphocytic leukemia.

Inhibitors of Bruton tyrosine kinase (BTK) and phosphatidylinositol 3-kinase δ (PI3Kδ) that target the B-cell receptor (BCR) signaling pathway have revolutionized the treatment of chronic lymphocytic leukemia (CLL). Mutations associated with resistance to BTK inhibitors have been identified, but limited data are available on mechanisms of resistance to PI3Kδ inhibitors. Here we present findings from longitudinal whole-exome sequencing of cells from patients with multiply relapsed CLL (N = 28) enrolled in trials of PI3K inhibitors. The nonresponder subgroup was characterized by baseline activating mutations in MAP2K1, BRAF, and KRAS genes in 60% of patients. PI3Kδ inhibition failed to inhibit ERK phosphorylation (pERK) in nonresponder CLL cells with and without mutations, whereas treatment with a MEK inhibitor rescued ERK inhibition. Overexpression of MAP2K1 mutants in vitro led to increased basal and inducible pERK and resistance to idelalisib. These data demonstrate that MAPK/ERK activation plays a key role in resistance to PI3Kδ inhibitors in CLL and provide a rationale for therapy with a combination of PI3Kδ and ERK inhibitors.

Distinct evolutionary paths in chronic lymphocytic leukemia during resistance to the graft-versus-leukemia effect.

Leukemic relapse remains a major barrier to successful allogeneic hematopoietic stem cell transplantation (allo-HSCT) for aggressive hematologic malignancies. The basis for relapse of advanced lymphoid malignancies remains incompletely understood and may involve escape from the graft-versus-leukemia (GvL) effect. We hypothesized that for patients with chronic lymphocytic leukemia (CLL) treated with allo-HSCT, leukemic cell-intrinsic features influence transplant outcomes by directing the evolutionary trajectories of CLL cells. Integrated genetic, transcriptomic, and epigenetic analyses of CLL cells from 10 patients revealed that the clinical kinetics of post-HSCT relapse are shaped by distinct molecular dynamics. Early relapses after allo-HSCT exhibited notable genetic stability; single CLL cell transcriptional analysis demonstrated a cellular heterogeneity that was static over time. In contrast, CLL cells relapsing late after allo-HSCT displayed notable genetic evolution and evidence of neoantigen depletion, consistent with marked single-cell transcriptional shifts that were unique to each patient. We observed a greater rate of epigenetic change for late relapses not seen in early relapses or relapses after chemotherapy alone, suggesting that the selection pressures of the GvL bottleneck are unlike those imposed by chemotherapy. No selective advantage for human leukocyte antigen (HLA) loss was observed, even when present in pretransplant subpopulations. Gain of stem cell modules was a common signature associated with leukemia relapse regardless of posttransplant relapse kinetics. These data elucidate the biological pathways that underlie GvL resistance and posttransplant relapse.

The cationic amino acid exporter Slc7a7 is induced and vital in zebrafish tissue macrophages with sustained efferocytic activity.

Most tissues harbor a substantial population of resident macrophages. Here, we elucidate a functional link between the Slc7a7 cationic amino acid transporter and tissue macrophages. We identified a mutant zebrafish devoid of microglia due to a mutation in the slc7a7 gene. We found that in Slc7a7-deficient larvae, macrophages do enter the retina and brain to become microglia, but then die during the developmental wave of neuronal apoptosis, which triggers intense efferocytic work from them. A similar macrophage demise occurs in other tissues, at stages where macrophages have to engulf many cell corpses, whether due to developmental or experimentally triggered cell death. We found that Slc7a7 is the main cationic amino acid transporter expressed in macrophages of zebrafish larvae, and that its expression is induced in tissue macrophages within 1-2 h upon efferocytosis. Our data indicate that Slc7a7 is vital not only for microglia but also for any steadily efferocytic tissue macrophages, and that slc7a7 gene induction is one of the adaptive responses that allow them to cope with the catabolism of numerous dead cells without compromising their own viability.