Lineage-specific canonical and non-canonical activity of EZH2 in advanced prostate cancer subtypes.

Enhancer of zeste homolog 2 (EZH2) is a histone methyltransferase and emerging therapeutic target that is overexpressed in most castration-resistant prostate cancers and implicated as a driver of disease progression and resistance to hormonal therapies. Here we define the lineage-specific action and differential activity of EZH2 in both prostate adenocarcinoma and neuroendocrine prostate cancer (NEPC) subtypes of advanced prostate cancer to better understand the role of EZH2 in modulating differentiation, lineage plasticity, and to identify mediators of response and resistance to EZH2 inhibitor therapy. Mechanistically, EZH2 modulates bivalent genes that results in upregulation of NEPC-associated transcriptional drivers (e.g., ASCL1) and neuronal gene programs in NEPC, and leads to forward differentiation after targeting EZH2 in NEPC. Subtype-specific downstream effects of EZH2 inhibition on cell cycle genes support the potential rationale for co-targeting cyclin/CDK to overcome resistance to EZH2 inhibition.

Parsing digital or analog TCR performance through piconewton forces.

αß T cell receptors (TCRs) principally recognize aberrant peptides bound to major histocompatibility complex molecules (pMHCs) on unhealthy cells, amplifying specificity and sensitivity through physical load placed on the TCR-pMHC bond during immunosurveillance. To understand this mechanobiology, TCRs stimulated by abundantly and sparsely arrayed epitopes (NP366-374/Db and PA224-233/Db, respectively) following in vivo influenza A virus infection were studied with optical tweezers. While certain NP repertoire CD8 T lymphocytes require many ligands for activation, others are digital, needing just few. Conversely, all PA TCRs perform digitally, exhibiting pronounced bond lifetime increases through sustained, energizing volleys of structural transitioning. Optimal digital performance is superior in vivo, correlating with ERK phosphorylation, CD3 loss, and activation marker up-regulation in vitro. Given neoantigen array paucity, digital TCRs are likely critical for immunotherapies.

Lineage-specific canonical and non-canonical activity of EZH2 in advanced prostate cancer subtypes.

Enhancer of zeste homolog 2 (EZH2) is a histone methyltransferase and emerging therapeutic target that is overexpressed in most castration-resistant prostate cancers and implicated as a driver of disease progression and resistance to hormonal therapies. Here we define the lineage-specific action and differential activity of EZH2 in both prostate adenocarcinoma (PRAD) and neuroendocrine prostate cancer (NEPC) subtypes of advanced prostate cancer to better understand the role of EZH2 in modulating differentiation, lineage plasticity, and to identify mediators of response and resistance to EZH2 inhibitor therapy. Mechanistically, EZH2 modulates bivalent genes that results in upregulation of NEPC-associated transcriptional drivers (e.g., ASCL1) and neuronal gene programs, and leads to forward differentiation after targeting EZH2 in NEPC. Subtype-specific downstream effects of EZH2 inhibition on cell cycle genes support the potential rationale for co-targeting cyclin/CDK to overcome resistance to EZH2 inhibition.

KLRG1 Cell Depletion as a Novel Therapeutic Strategy in Patients with Mature T-Cell Lymphoma Subtypes.

PURPOSE: Develop a novel therapeutic strategy for patients with subtypes of mature T-cell and NK-cell neoplasms. EXPERIMENTAL DESIGN: Primary specimens, cell lines, patient-derived xenograft models, commercially available, and proprietary anti-KLRG1 antibodies were used for screening, target, and functional validation. RESULTS: Here we demonstrate that surface KLRG1 is highly expressed on tumor cells in subsets of patients with extranodal NK/T-cell lymphoma (ENKTCL), T-prolymphocytic leukemia (T-PLL), and gamma/delta T-cell lymphoma (G/D TCL). The majority of the CD8+/CD57+ or CD3-/CD56+ leukemic cells derived from patients with T- and NK-large granular lymphocytic leukemia (T-LGLL and NK-LGLL), respectively, expressed surface KLRG1. The humanized afucosylated anti-KLRG1 monoclonal antibody (mAb208) optimized for mouse in vivo use depleted KLRG1+ TCL cells by mechanisms of ADCC, ADCP, and CDC rather than apoptosis. mAb208 induced ADCC and ADCP of T-LGLL patient-derived CD8+/CD57+ cells ex vivo. mAb208 effected ADCC of subsets of healthy donor-derived KLRG1+ NK, CD4+, CD8+ Tem, and TemRA cells while sparing KLRG1- naïve and CD8+ Tcm cells. Treatment of cell line and TCL patient-derived xenografts with mAb208 or anti-CD47 mAb alone and in combination with the PI3K-δ/γ inhibitor duvelisib extended survival. The depletion of macrophages in vivo antagonized mAb208 efficacy. CONCLUSIONS: Our findings suggest the potential benefit of a broader treatment strategy combining therapeutic antibodies with PI3Ki for the treatment of patients with mature T-cell and NK-cell neoplasms. See related commentary by Varma and Diefenbach, p. 2300.

Parsing digital or analogue TCR performance through piconewton forces.

αß T-cell receptors (TCRs) recognize aberrant peptides bound to major histocompatibility complex molecules (pMHCs) on unhealthy cells, amplifying specificity and sensitivity through physical load placed on the TCR-pMHC bond during immunosurveillance. To understand this mechanobiology, TCRs stimulated by abundantly and sparsely arrayed epitopes (NP 366-374 /D b and PA 224-233 /D b , respectively) following in vivo influenza A virus infection were studied with optical tweezers. While certain NP repertoire CD8 T lymphocytes require many ligands for activation, others are digital, needing just few. Conversely, all PA TCRs perform digitally, exhibiting pronounced bond lifetime increases through sustained, energizing volleys of structural transitioning. Optimal digital performance is superior in vivo, correlating with ERK phosphorylation, CD3 loss, and activation marker upregulation in vitro . Given neoantigen array paucity, digital TCRs are likely critical for immunotherapies. One Sentence Summary: Quality of ligand recognition in a T-cell repertoire is revealed through application of physical load on clonal T-cell receptor (TCR)-pMHC bonds.

KDM6A epigenetically regulates subtype plasticity in small cell lung cancer.

Small cell lung cancer (SCLC) exists broadly in four molecular subtypes: ASCL1, NEUROD1, POU2F3 and Inflammatory. Initially, SCLC subtypes were thought to be mutually exclusive, but recent evidence shows intra-tumoural subtype heterogeneity and plasticity between subtypes. Here, using a CRISPR-based autochthonous SCLC genetically engineered mouse model to study the consequences of KDM6A/UTX inactivation, we show that KDM6A inactivation induced plasticity from ASCL1 to NEUROD1 resulting in SCLC tumours that express both ASCL1 and NEUROD1. Mechanistically, KDM6A normally maintains an active chromatin state that favours the ASCL1 subtype with its loss decreasing H3K4me1 and increasing H3K27me3 at enhancers of neuroendocrine genes leading to a cell state that is primed for ASCL1-to-NEUROD1 subtype switching. This work identifies KDM6A as an epigenetic regulator that controls ASCL1 to NEUROD1 subtype plasticity and provides an autochthonous SCLC genetically engineered mouse model to model ASCL1 and NEUROD1 subtype heterogeneity and plasticity, which is found in 35-40% of human SCLCs.

Oncogenic Drivers and Therapeutic Vulnerabilities in KRAS Wild-Type Pancreatic Cancer.

PURPOSE: Approximately 8% to 10% of pancreatic ductal adenocarcinomas (PDAC) do not harbor mutations in KRAS. Understanding the unique molecular and clinical features of this subset of pancreatic cancer is important to guide patient stratification for clinical trials of molecularly targeted agents. EXPERIMENTAL DESIGN: We analyzed a single-institution cohort of 795 exocrine pancreatic cancer cases (including 785 PDAC cases) with a targeted multigene sequencing panel and identified 73 patients (9.2%) with KRAS wild-type (WT) pancreatic cancer. RESULTS: Overall, 43.8% (32/73) of KRAS WT cases had evidence of an alternative driver of the MAPK pathway, including BRAF mutations and in-frame deletions and receptor tyrosine kinase fusions. Conversely, 56.2% of cases did not harbor a clear MAPK driver alteration, but 29.3% of these MAPK-negative KRAS WT cases (12/41) demonstrated activating alterations in other oncogenic drivers, such as GNAS, MYC, PIK3CA, and CTNNB1. We demonstrate potent efficacy of pan-RAF and MEK inhibition in patient-derived organoid models carrying BRAF in-frame deletions. Moreover, we demonstrate durable clinical benefit of targeted therapy in a patient harboring a KRAS WT tumor with a ROS1 fusion. Clinically, patients with KRAS WT tumors were significantly younger in age of onset (median age: 62.6 vs. 65.7 years; P = 0.037). SMAD4 mutations were associated with a particularly poor prognosis in KRAS WT cases. CONCLUSIONS: This study defines the genomic underpinnings of KRAS WT pancreatic cancer and highlights potential therapeutic avenues for future investigation in molecularly directed clinical trials. See related commentary by Kato et al., p. 4527.

Single-cell molecular profiling using ex vivo functional readouts fuels precision oncology in glioblastoma.

BACKGROUND: Functional profiling of freshly isolated glioblastoma (GBM) cells is being evaluated as a next-generation method for precision oncology. While promising, its success largely depends on the method to evaluate treatment activity which requires sufficient resolution and specificity. METHODS: Here, we describe the 'precision oncology by single-cell profiling using ex vivo readouts of functionality' (PROSPERO) assay to evaluate the intrinsic susceptibility of high-grade brain tumor cells to respond to therapy. Different from other assays, PROSPERO extends beyond life/death screening by rapidly evaluating acute molecular drug responses at single-cell resolution. RESULTS: The PROSPERO assay was developed by correlating short-term single-cell molecular signatures using mass cytometry by time-of-flight (CyTOF) to long-term cytotoxicity readouts in representative patient-derived glioblastoma cell cultures (n = 14) that were exposed to radiotherapy and the small-molecule p53/MDM2 inhibitor AMG232. The predictive model was subsequently projected to evaluate drug activity in freshly resected GBM samples from patients (n = 34). Here, PROSPERO revealed an overall limited capacity of tumor cells to respond to therapy, as reflected by the inability to induce key molecular markers upon ex vivo treatment exposure, while retaining proliferative capacity, insights that were validated in patient-derived xenograft (PDX) models. This approach also allowed the investigation of cellular plasticity, which in PDCLs highlighted therapy-induced proneural-to-mesenchymal (PMT) transitions, while in patients' samples this was more heterogeneous. CONCLUSION: PROSPERO provides a precise way to evaluate therapy efficacy by measuring molecular drug responses using specific biomarker changes in freshly resected brain tumor samples, in addition to providing key functional insights in cellular behavior, which may ultimately complement standard, clinical biomarker evaluations.

PGE2 alters chromatin through H2A.Z-variant enhancer nucleosome modification to promote hematopoietic stem cell fate.

Prostaglandin E2 (PGE2) and 16,16-dimethyl-PGE2 (dmPGE2) are important regulators of hematopoietic stem and progenitor cell (HSPC) fate and offer potential to enhance stem cell therapies [C. Cutler et al. Blood 122, 3074-3081(2013); W. Goessling et al. Cell Stem Cell 8, 445-458 (2011); W. Goessling et al. Cell 136, 1136-1147 (2009)]. Here, we report that PGE2-induced changes in chromatin at enhancer regions through histone-variant H2A.Z permit acute inflammatory gene induction to promote HSPC fate. We found that dmPGE2-inducible enhancers retain MNase-accessible, H2A.Z-variant nucleosomes permissive of CREB transcription factor (TF) binding. CREB binding to enhancer nucleosomes following dmPGE2 stimulation is concomitant with deposition of histone acetyltransferases p300 and Tip60 on chromatin. Subsequent H2A.Z acetylation improves chromatin accessibility at stimuli-responsive enhancers. Our findings support a model where histone-variant nucleosomes retained within inducible enhancers facilitate TF binding. Histone-variant acetylation by TF-associated nucleosome remodelers creates the accessible nucleosome landscape required for immediate enhancer activation and gene induction. Our work provides a mechanism through which inflammatory mediators, such as dmPGE2, lead to acute transcriptional changes and modify HSPC behavior to improve stem cell transplantation.

Cancer-associated fibroblasts are the main contributors to epithelial-to-mesenchymal signatures in the tumor microenvironment.

Epithelial-to-mesenchymal transition (EMT) is associated with tumor initiation, metastasis, and drug resistance. However, the mechanisms underlying these associations are largely unknown. We studied several tumor types to identify the source of EMT gene expression signals and a potential mechanism of resistance to immuno-oncology treatment. Across tumor types, EMT-related gene expression was strongly associated with expression of stroma-related genes. Based on RNA sequencing of multiple patient-derived xenograft models, EMT-related gene expression was enriched in the stroma versus parenchyma. EMT-related markers were predominantly expressed by cancer-associated fibroblasts (CAFs), cells of mesenchymal origin which produce a variety of matrix proteins and growth factors. Scores derived from a 3-gene CAF transcriptional signature (COL1A1, COL1A2, COL3A1) were sufficient to reproduce association between EMT-related markers and disease prognosis. Our results suggest that CAFs are the primary source of EMT signaling and have potential roles as biomarkers and targets for immuno-oncology therapies.

ALK Amplification and Rearrangements Are Recurrent Targetable Events in Congenital and Adult Glioblastoma.

PURPOSE: Anaplastic lymphoma kinase (ALK) aberrations have been identified in pediatric-type infant gliomas, but their occurrence across age groups, functional effects, and treatment response has not been broadly established. EXPERIMENTAL DESIGN: We performed a comprehensive analysis of ALK expression and genomic aberrations in both newly generated and retrospective data from 371 glioblastomas (156 adult, 205 infant/pediatric, and 10 congenital) with in vitro and in vivo validation of aberrations. RESULTS: ALK aberrations at the protein or genomic level were detected in 12% of gliomas (45/371) in a wide age range (0-80 years). Recurrent as well as novel ALK fusions (LRRFIP1-ALK, DCTN1-ALK, PRKD3-ALK) were present in 50% (5/10) of congenital/infant, 1.4% (3/205) of pediatric, and 1.9% (3/156) of adult GBMs. ALK fusions were present as the only candidate driver in congenital/infant GBMs and were sometimes focally amplified. In contrast, adult ALK fusions co-occurred with other oncogenic drivers. No activating ALK mutations were identified in any age group. Novel and recurrent ALK rearrangements promoted STAT3 and ERK1/2 pathways and transformation in vitro and in vivo. ALK-fused GBM cellular and mouse models were responsive to ALK inhibitors, including in patient cells derived from a congenital GBM. Relevant to the treatment of infant gliomas, we showed that ALK protein appears minimally expressed in the forebrain at perinatal stages, and no gross effects on perinatal brain development were seen in pregnant mice treated with the ALK inhibitor ceritinib. CONCLUSIONS: These findings support use of brain-penetrant ALK inhibitors in clinical trials across infant, pediatric, and adult GBMs. See related commentary by Mack and Bertrand, p. 2567.

Disruption of polyhomeotic polymerization decreases nucleosome occupancy and alters genome accessibility.

Chromatin attains its three-dimensional (3D) conformation by establishing contacts between different noncontiguous regions. Sterile Alpha Motif (SAM)-mediated polymerization of the polyhomeotic (PH) protein regulates subnuclear clustering of Polycomb Repressive Complex 1 (PRC1) and chromatin topology. The mutations that perturb the ability of the PH to polymerize, disrupt long-range chromatin contacts, alter Hox gene expression, and lead to developmental defects. To understand the underlying mechanism, we combined the experiments and theory to investigate the effect of this SAM domain mutation on nucleosome occupancy and accessibility on a genome wide scale. Our data show that disruption of PH polymerization because of SAM domain mutation decreases nucleosome occupancy and alters accessibility. Polymer simulations investigating the interplay between distant chromatin contacts and nucleosome occupancy, both of which are regulated by PH polymerization, suggest that nucleosome density increases when contacts between different regions of chromatin are established. Taken together, it appears that SAM domain-mediated PH polymerization biomechanically regulates the organization of chromatin at multiple scales from nucleosomes to chromosomes and we suggest that higher order organization can have a top-down causation effect on nucleosome occupancy.

Genomic and biological study of fusion genes as resistance mechanisms to EGFR inhibitors.

The clinical significance of gene fusions detected by DNA-based next generation sequencing remains unclear as resistance mechanisms to EGFR tyrosine kinase inhibitors in EGFR mutant non-small cell lung cancer. By studying EGFR inhibitor-resistant patients treated with a combination of an EGFR inhibitor and a drug targeting the putative resistance-causing fusion oncogene, we identify patients who benefit and those who do not from this treatment approach. Through evaluation including RNA-seq of potential drug resistance-imparting fusion oncogenes in 504 patients with EGFR mutant lung cancer, we identify only a minority of them as functional, potentially capable of imparting EGFR inhibitor resistance. We further functionally validate fusion oncogenes in vitro using CRISPR-based editing of EGFR mutant cell lines and use these models to identify known and unknown drug resistance mechanisms to combination therapies. Collectively, our results partially reveal the complex nature of fusion oncogenes as potential drug resistance mechanisms and highlight approaches that can be undertaken to determine their functional significance.

Regulation of neuroendocrine plasticity by the RNA-binding protein ZFP36L1.

Some small cell lung cancers (SCLCs) are highly sensitive to inhibitors of the histone demethylase LSD1. LSD1 inhibitors are thought to induce their anti-proliferative effects by blocking neuroendocrine differentiation, but the mechanisms by which LSD1 controls the SCLC neuroendocrine phenotype are not well understood. To identify genes required for LSD1 inhibitor sensitivity in SCLC, we performed a positive selection genome-wide CRISPR/Cas9 loss of function screen and found that ZFP36L1, an mRNA-binding protein that destabilizes mRNAs, is required for LSD1 inhibitor sensitivity. LSD1 binds and represses ZFP36L1 and upon LSD1 inhibition, ZFP36L1 expression is restored, which is sufficient to block the SCLC neuroendocrine differentiation phenotype and induce a non-neuroendocrine "inflammatory" phenotype. Mechanistically, ZFP36L1 binds and destabilizes SOX2 and INSM1 mRNAs, two transcription factors that are required for SCLC neuroendocrine differentiation. This work identifies ZFP36L1 as an LSD1 target gene that controls the SCLC neuroendocrine phenotype and demonstrates that modulating mRNA stability of lineage transcription factors controls neuroendocrine to non-neuroendocrine plasticity.

Association of High Tumor Mutation Burden in Non-Small Cell Lung Cancers With Increased Immune Infiltration and Improved Clinical Outcomes of PD-L1 Blockade Across PD-L1 Expression Levels.

Importance: Although tumor mutation burden (TMB) has been explored as a potential biomarker of immunotherapy efficacy in solid tumors, there still is a lack of consensus about the optimal TMB threshold that best discriminates improved outcomes of immune checkpoint inhibitor therapy among patients with non-small cell lung cancer (NSCLC). Objectives: To determine the association between increasing TMB levels and immunotherapy efficacy across clinically relevant programmed death ligand-1 (PD-L1) levels in patients with NSCLC. Design, Setting, and Participants: This multicenter cohort study included patients with advanced NSCLC treated with immunotherapy who received programmed cell death-1 (PD-1) or PD-L1 inhibition in the Dana-Farber Cancer Institute (DFCI), Memorial Sloan Kettering Cancer Center (MSKCC), and in the Stand Up To Cancer (SU2C)/Mark Foundation data sets. Clinicopathological and genomic data were collected from patients between September 2013 and September 2020. Data analysis was performed from November 2021 to February 2022. Exposures: Treatment with PD-1/PD-L1 inhibition without chemotherapy. Main Outcomes and Measures: Association of TMB levels with objective response rate (ORR), progression-free survival (PFS), and overall survival (OS). Results: In the entire cohort of 1552 patients with advanced NSCLC who received PD-1/PD-L1 blockade, the median (range) age was 66 (22-92) years, 830 (53.5%) were women, and 1347 (86.8%) had cancer with nonsquamous histologic profile. A regression tree modeling ORR as a function of TMB identified 2 TMB groupings in the discovery cohort (MSKCC), defined as low TMB (≤19.0 mutations per megabase) and high TMB (>19.0 mutations per megabase), which were associated with increasing improvements in ORR, PFS, and OS in the discovery cohort and in 2 independent cohorts (DFCI and SU2C/Mark Foundation). These TMB levels also were associated with significant improvements in outcomes of immunotherapy in each PD-L1 tumor proportion score subgroup of less than 1%, 1% to 49%, and 50% or higher. The ORR to PD-1/PD-L1 inhibition was as high as 57% in patients with high TMB and PD-L1 expression 50% or higher and as low as 8.7% in patients with low TMB and PD-L1 expression less than 1%. Multiplexed immunofluorescence and transcriptomic profiling revealed that high TMB levels were associated with increased CD8-positive, PD-L1-positive T-cell infiltration, increased PD-L1 expression on tumor and immune cells, and upregulation of innate and adaptive immune response signatures. Conclusions and Relevance: These findings suggest that increasing TMB levels are associated with immune cell infiltration and an inflammatory T-cell-mediated response, resulting in increased sensitivity to PD-1/PD-L1 blockade in NSCLC across PD-L1 expression subgroups.

Activation of Tumor-Cell STING Primes NK-Cell Therapy.

Activation of the stimulator of interferon genes (STING) pathway promotes antitumor immunity but STING agonists have yet to achieve clinical success. Increased understanding of the mechanism of action of STING agonists in human tumors is key to developing therapeutic combinations that activate effective innate antitumor immunity. Here, we report that malignant pleural mesothelioma cells robustly express STING and are responsive to STING agonist treatment ex vivo. Using dynamic single-cell RNA sequencing of explants treated with a STING agonist, we observed CXCR3 chemokine activation primarily in tumor cells and cancer-associated fibroblasts, as well as T-cell cytotoxicity. In contrast, primary natural killer (NK) cells resisted STING agonist-induced cytotoxicity. STING agonists enhanced migration and killing of NK cells and mesothelin-targeted chimeric antigen receptor (CAR)-NK cells, improving therapeutic activity in patient-derived organotypic tumor spheroids. These studies reveal the fundamental importance of using human tumor samples to assess innate and cellular immune therapies. By functionally profiling mesothelioma tumor explants with elevated STING expression in tumor cells, we uncovered distinct consequences of STING agonist treatment in humans that support testing combining STING agonists with NK and CAR-NK cell therapies.

Astrocyte-neuron crosstalk through Hedgehog signaling mediates cortical synapse development.

Neuron-glia interactions play a critical role in the regulation of synapse formation and circuit assembly. Here we demonstrate that canonical Sonic hedgehog (Shh) pathway signaling in cortical astrocytes acts to coordinate layer-specific synaptic connectivity. We show that the Shh receptor Ptch1 is expressed by cortical astrocytes during development and that Shh signaling is necessary and sufficient to promote the expression of genes involved in regulating synaptic development and layer-enriched astrocyte molecular identity. Loss of Shh in layer V neurons reduces astrocyte complexity and coverage by astrocytic processes in tripartite synapses; conversely, cell-autonomous activation of Shh signaling in astrocytes promotes cortical excitatory synapse formation. Furthermore, Shh-dependent genes Lrig1 and Sparc distinctively contribute to astrocyte morphology and synapse formation. Together, these results suggest that Shh secreted from deep-layer cortical neurons acts to specialize the molecular and functional features of astrocytes during development to shape circuit assembly and function.

Standards for the classification of pathogenicity of somatic variants in cancer (oncogenicity): Joint recommendations of Clinical Genome Resource (ClinGen), Cancer Genomics Consortium (CGC), and Variant Interpretation for Cancer Consortium (VICC).

PURPOSE: Several professional societies have published guidelines for the clinical interpretation of somatic variants, which specifically address diagnostic, prognostic, and therapeutic implications. Although these guidelines for the clinical interpretation of variants include data types that may be used to determine the oncogenicity of a variant (eg, population frequency, functional, and in silico data or somatic frequency), they do not provide a direct, systematic, and comprehensive set of standards and rules to classify the oncogenicity of a somatic variant. This insufficient guidance leads to inconsistent classification of rare somatic variants in cancer, generates variability in their clinical interpretation, and, importantly, affects patient care. Therefore, it is essential to address this unmet need. METHODS: Clinical Genome Resource (ClinGen) Somatic Cancer Clinical Domain Working Group and ClinGen Germline/Somatic Variant Subcommittee, the Cancer Genomics Consortium, and the Variant Interpretation for Cancer Consortium used a consensus approach to develop a standard operating procedure (SOP) for the classification of oncogenicity of somatic variants. RESULTS: This comprehensive SOP has been developed to improve consistency in somatic variant classification and has been validated on 94 somatic variants in 10 common cancer-related genes. CONCLUSION: The comprehensive SOP is now available for classification of oncogenicity of somatic variants.