Machine learning identifies activation of RUNX/AP-1 as drivers of mesenchymal and fibrotic regulatory programs in gastric cancer.

Gastric cancer (GC) is the fifth most common cancer worldwide and is a heterogeneous disease. Among GC subtypes, the mesenchymal phenotype (Mes-like) is more invasive than the epithelial phenotype (Epi-like). Although gene expression of the epithelial-to-mesenchymal transition (EMT) has been studied, the regulatory landscape shaping this process is not fully understood. Here we use ATAC-seq and RNA-seq data from a compendium of GC cell lines and primary tumors to detect drivers of regulatory state changes and their transcriptional responses. Using the ATAC-seq data, we developed a machine learning approach to determine the transcription factors (TFs) regulating the subtypes of GC. We identified TFs driving the mesenchymal (RUNX2, ZEB1, SNAI2, AP-1 dimer) and the epithelial (GATA4, GATA6, KLF5, HNF4A, FOXA2, GRHL2) states in GC. We identified DNA copy number alterations associated with dysregulation of these TFs, specifically deletion of GATA4 and amplification of MAPK9 Comparisons with bulk and single-cell RNA-seq data sets identified activation toward fibroblast-like epigenomic and expression signatures in Mes-like GC. The activation of this mesenchymal fibrotic program is associated with differentially accessible DNA cis-regulatory elements flanking upregulated mesenchymal genes. These findings establish a map of TF activity in GC and highlight the role of copy number driven alterations in shaping epigenomic regulatory programs as potential drivers of GC heterogeneity and progression.

Detection of circulating tumor DNA with ultradeep sequencing of plasma cell-free DNA for monitoring minimal residual disease and early detection of recurrence in early-stage lung cancer.

BACKGROUND: In early-stage non-small cell lung cancer (NSCLC), recurrence is frequently observed. Circulating tumor DNA (ctDNA) has emerged as a noninvasive tool to risk stratify patients for recurrence after curative intent therapy. This study aimed to risk stratify patients with early-stage NSCLC via a personalized, tumor-informed multiplex polymerase chain reaction (mPCR) next-generation sequencing assay. METHODS: This retrospective cohort study included patients with stage I-III NSCLC. Recruited patients received standard-of-care management (surgical resection with or without adjuvant chemotherapy, followed by surveillance). Whole-exome sequencing of NSCLC resected tissue and matched germline DNA was used to design patient-specific mPCR assays (Signatera, Natera, Inc) to track up to 16 single-nucleotide variants in plasma samples. RESULTS: The overall cohort with analyzed plasma samples consisted of 57 patients. Stage distribution was 68% for stage I and 16% each for stages II and III. Presurgery (i.e., at baseline), ctDNA was detected in 15 of 57 patients (26%). ctDNA detection presurgery was significantly associated with shorter recurrence-free survival (RFS; hazard ratio [HR], 3.54; 95% confidence interval [CI], 1.00-12.62; p = .009). In the postsurgery setting, ctDNA was detected in seven patients, of whom 100% experienced radiological recurrence. ctDNA positivity preceded radiological findings by a median lead time of 2.8 months (range, 0-12.9 months). Longitudinally, ctDNA detection at any time point was associated with shorter RFS (HR, 16.1; 95% CI, 1.63-158.9; p < .0001). CONCLUSIONS: ctDNA detection before surgical resection was strongly associated with a high risk of relapse in early-stage NSCLC in a large unique Asian cohort. Prospective studies are needed to assess the clinical utility of ctDNA status in this setting.

A multimodal atlas of hepatocellular carcinoma reveals convergent evolutionary paths and 'bad apple' effect on clinical trajectory.

BACKGROUND & AIMS: Hepatocellular Carcinoma (HCC) is a highly fatal cancer characterized by high intra-tumor heterogeneity (ITH). A panoramic understanding of its tumor evolution, in relation to its clinical trajectory, may provide novel prognostic and treatment strategies. METHODS: Through the Asia-Pacific Hepatocellular Carcinoma (AHCC) trials group (NCT03267641), we recruited one of the largest prospective cohorts of HCC with over 600 whole genome and transcriptome samples from 123 treatment-naïve patients. RESULTS: Using a multi-region sampling approach, we revealed seven convergent genetic evolutionary paths governed by the early driver mutations, late copy number variations and viral integrations, which stratify patient clinical trajectories after surgical resection. Furthermore, such evolutionary paths shaped the molecular profiles, leading to distinct transcriptomic subtypes. Most significantly, although we found the coexistence of multiple transcriptomic subtypes within certain tumors, patient prognosis was best predicted by the most aggressive cell fraction of the tumor, rather than by overall degree of transcriptomic ITH level - a phenomenon we termed the 'bad apple' effect. Finally, we found that characteristics throughout early and late tumor evolution provide significant and complementary prognostic power in predicting patient survival. CONCLUSIONS: Taken together, our study generated a comprehensive landscape of evolutionary history for HCC and provided a rich multi-omics resource for understanding tumor heterogeneity and clinical trajectories. CLINICAL TRIAL NUMBER: NCT03267641 (Observational cohort) IMPACT AND IMPLICATIONS: This prospective study, utilizing comprehensive multi-sector, multi-omics sequencing and clinical data from surgically resected HCC, reveals critical insights into the role of tumor evolution and intra-tumor heterogeneity (ITH) in determining the prognosis of Hepatocellular Carcinoma (HCC). These findings are invaluable for oncology researchers and clinicians, as they underscore the influence of distinct evolutionary paths and the 'bad apple' effect, where the most aggressive tumor fraction dictates disease progression. These insights not only enhance prognostic accuracy post-surgical resection but also pave the way for developing personalized therapies tailored to specific tumor evolutionary and transcriptomic profiles. The co-existence of multiple sub-types within the same tumor prompts a re-appraisal of the utilities of depending on single samples to represent the entire tumor and suggests the need for clinical molecular imaging. This research thus marks a significant step forward in the clinical understanding and management of HCC, underscoring the importance of integrating tumor evolutionary dynamics and multi-omics biomarkers into therapeutic decision-making.

Landscape of germline pathogenic variants in patients with dual primary breast and lung cancer.

BACKGROUND: Cancer predisposition is most often studied in the context of single cancers. However, inherited cancer predispositions can also give rise to multiple primary cancers. Yet, there is a paucity of studies on genetic predisposition in multiple primary cancers, especially those outside of well-defined cancer predisposition syndromes. This study aimed to identify germline variants associated with dual primary cancers of the breast and lung. METHODS: Exome sequencing was performed on germline DNA from 55 Singapore patients (52 [95%] never-smokers) with dual primaries in the breast and lung, confirmed by histopathology. Using two large control cohorts: the local SG10K_Health (n = 9770) and gnomAD non-cancer East Asians (n = 9626); and two additional local case cohorts of early-onset or familial breast cancer (n = 290), and lung cancer (n = 209), variants were assessed for pathogenicity in accordance with ACMG/AMP guidelines. In particular, comparisons were made with known pathogenic or likely pathogenic variants in the ClinVar database, pathogenicity predictions were obtained from in silico prediction software, and case-control association analyses were performed. RESULTS: Altogether, we identified 19 pathogenic or likely pathogenic variants from 16 genes, detected in 17 of 55 (31%) patients. Six of the 19 variants were identified using ClinVar, while 13 variants were classified pathogenic or likely pathogenic using ACMG/AMP guidelines. The 16 genes include well-known cancer predisposition genes such as BRCA2, TP53, and RAD51D; but also lesser known cancer genes EXT2, WWOX, GATA2, and GPC3. Most of these genes are involved in DNA damage repair, reaffirming the role of impaired DNA repair mechanisms in the development of multiple malignancies. These variants warrant further investigations in additional populations. CONCLUSIONS: We have identified both known and novel variants significantly enriched in patients with primary breast and lung malignancies, expanding the body of known cancer predisposition variants for both breast and lung cancer. These variants are mostly from genes involved in DNA repair, affirming the role of impaired DNA repair in the predisposition and development of multiple cancers.

Regulatory enhancer profiling of mesenchymal-type gastric cancer reveals subtype-specific epigenomic landscapes and targetable vulnerabilities.

OBJECTIVE: Gastric cancer (GC) comprises multiple molecular subtypes. Recent studies have highlighted mesenchymal-subtype GC (Mes-GC) as a clinically aggressive subtype with few treatment options. Combining multiple studies, we derived and applied a consensus Mes-GC classifier to define the Mes-GC enhancer landscape revealing disease vulnerabilities. DESIGN: Transcriptomic profiles of ~1000 primary GCs and cell lines were analysed to derive a consensus Mes-GC classifier. Clinical and genomic associations were performed across >1200 patients with GC. Genome-wide epigenomic profiles (H3K27ac, H3K4me1 and assay for transposase-accessible chromatin with sequencing (ATAC-seq)) of 49 primary GCs and GC cell lines were generated to identify Mes-GC-specific enhancer landscapes. Upstream regulators and downstream targets of Mes-GC enhancers were interrogated using chromatin immunoprecipitation followed by sequencing (ChIP-seq), RNA sequencing, CRISPR/Cas9 editing, functional assays and pharmacological inhibition. RESULTS: We identified and validated a 993-gene cancer-cell intrinsic Mes-GC classifier applicable to retrospective cohorts or prospective single samples. Multicohort analysis of Mes-GCs confirmed associations with poor patient survival, therapy resistance and few targetable genomic alterations. Analysis of enhancer profiles revealed a distinctive Mes-GC epigenomic landscape, with TEAD1 as a master regulator of Mes-GC enhancers and Mes-GCs exhibiting preferential sensitivity to TEAD1 pharmacological inhibition. Analysis of Mes-GC super-enhancers also highlighted NUAK1 kinase as a downstream target, with synergistic effects observed between NUAK1 inhibition and cisplatin treatment. CONCLUSION: Our results establish a consensus Mes-GC classifier applicable to multiple transcriptomic scenarios. Mes-GCs exhibit a distinct epigenomic landscape, and TEAD1 inhibition and combinatorial NUAK1 inhibition/cisplatin may represent potential targetable options.

Comprehensive molecular phenotyping of ARID1A-deficient gastric cancer reveals pervasive epigenomic reprogramming and therapeutic opportunities.

OBJECTIVE: Gastric cancer (GC) is a leading cause of cancer mortality, with ARID1A being the second most frequently mutated driver gene in GC. We sought to decipher ARID1A-specific GC regulatory networks and examine therapeutic vulnerabilities arising from ARID1A loss. DESIGN: Genomic profiling of GC patients including a Singapore cohort (>200 patients) was performed to derive mutational signatures of ARID1A inactivation across molecular subtypes. Single-cell transcriptomic profiles of ARID1A-mutated GCs were analysed to examine tumour microenvironmental changes arising from ARID1A loss. Genome-wide ARID1A binding and chromatin profiles (H3K27ac, H3K4me3, H3K4me1, ATAC-seq) were generated to identify gastric-specific epigenetic landscapes regulated by ARID1A. Distinct cancer hallmarks of ARID1A-mutated GCs were converged at the genomic, single-cell and epigenomic level, and targeted by pharmacological inhibition. RESULTS: We observed prevalent ARID1A inactivation across GC molecular subtypes, with distinct mutational signatures and linked to a NFKB-driven proinflammatory tumour microenvironment. ARID1A-depletion caused loss of H3K27ac activation signals at ARID1A-occupied distal enhancers, but unexpectedly gain of H3K27ac at ARID1A-occupied promoters in genes such as NFKB1 and NFKB2. Promoter activation in ARID1A-mutated GCs was associated with enhanced gene expression, increased BRD4 binding, and reduced HDAC1 and CTCF occupancy. Combined targeting of promoter activation and tumour inflammation via bromodomain and NFKB inhibitors confirmed therapeutic synergy specific to ARID1A-genomic status. CONCLUSION: Our results suggest a therapeutic strategy for ARID1A-mutated GCs targeting both tumour-intrinsic (BRD4-assocatiated promoter activation) and extrinsic (NFKB immunomodulation) cancer phenotypes.

Epigenetic promoter alterations in GI tumour immune-editing and resistance to immune checkpoint inhibition.

OBJECTIVES: Epigenomic alterations in cancer interact with the immune microenvironment to dictate tumour evolution and therapeutic response. We aimed to study the regulation of the tumour immune microenvironment through epigenetic alternate promoter use in gastric cancer and to expand our findings to other gastrointestinal tumours. DESIGN: Alternate promoter burden (APB) was quantified using a novel bioinformatic algorithm (proActiv) to infer promoter activity from short-read RNA sequencing and samples categorised into APBhigh, APBint and APBlow. Single-cell RNA sequencing was performed to analyse the intratumour immune microenvironment. A humanised mouse cancer in vivo model was used to explore dynamic temporal interactions between tumour kinetics, alternate promoter usage and the human immune system. Multiple cohorts of gastrointestinal tumours treated with immunotherapy were assessed for correlation between APB and treatment outcomes. RESULTS: APBhigh gastric cancer tumours expressed decreased levels of T-cell cytolytic activity and exhibited signatures of immune depletion. Single-cell RNAsequencing analysis confirmed distinct immunological populations and lower T-cell proportions in APBhigh tumours. Functional in vivo studies using 'humanised mice' harbouring an active human immune system revealed distinct temporal relationships between APB and tumour growth, with APBhigh tumours having almost no human T-cell infiltration. Analysis of immunotherapy-treated patients with GI cancer confirmed resistance of APBhigh tumours to immune checkpoint inhibition. APBhigh gastric cancer exhibited significantly poorer progression-free survival compared with APBlow (median 55 days vs 121 days, HR 0.40, 95% CI 0.18 to 0.93, p=0.032). CONCLUSION: These findings demonstrate an association between alternate promoter use and the tumour microenvironment, leading to immune evasion and immunotherapy resistance.

A phase 1b study of OXIRI in pancreatic adenocarcinoma patients and its immunomodulatory effects.

Pancreatic ductal adenocarcinoma (PDAC) is a highly lethal and debilitating disease with limited therapeutic options. The aim of this clinical study was to evaluate the safety, efficacy and pharmacokinetics of a novel regimen comprised of metronomic oxaliplatin (O), chronomodulated capecitabine (X) and UGT1A1 genotype-guided dosing of irinotecan (IRI) [OXIRI] as well as its immunomodulatory effects. Thirty-six patients were enrolled into either dose-escalation or expansion cohorts. In the dose escalation phase, capecitabine doses (2000, 2650, 3500 and 4500 mg/day) were administered at midnight on days 1 to 14 while oxaliplatin and irinotecan were intravenously infused at fixed doses of 50 and 75 mg/m2 respectively on days 1, 8 in a 21-day cycle. The maximum tolerated dose of capecitabine was 2650 mg/day and the most common grade 3 adverse events were neutropenia (30.6%) and diarrhea (13.9%). No grade 4 toxicity was observed. UGT1A1-genotype directed dosing resulted in similar exposure levels of irinotecan, SN-38 and SN-38G in all patients. Objective response rate was 22.2%. Median overall survival and progression-free survival were 8.1 and 5.2 months, respectively. Exploratory immunoprofiling by flow cytometry and quantitative spatial localization analysis of infiltrated immune cells performed on biopsy and plasma samples revealed significant declines in CCL22, CCL2 and TNFα levels at end of first cycle and an active host immune response. Our study showed that OXIRI was well-tolerated and exhibited good efficacy, with immunomodulatory effects. It may be considered as an alternative to FOLFIRINOX in patients intolerant to the latter.

Integrated paired-end enhancer profiling and whole-genome sequencing reveals recurrent CCNE1 and IGF2 enhancer hijacking in primary gastric adenocarcinoma.

OBJECTIVE: Genomic structural variations (SVs) causing rewiring of cis-regulatory elements remain largely unexplored in gastric cancer (GC). To identify SVs affecting enhancer elements in GC (enhancer-based SVs), we integrated epigenomic enhancer profiles revealed by paired-end H3K27ac ChIP-sequencing from primary GCs with tumour whole-genome sequencing (WGS) data (PeNChIP-seq/WGS). DESIGN: We applied PeNChIP-seq to 11 primary GCs and matched normal tissues combined with WGS profiles of >200 GCs. Epigenome profiles were analysed alongside matched RNA-seq data to identify tumour-associated enhancer-based SVs with altered cancer transcription. Functional validation of candidate enhancer-based SVs was performed using CRISPR/Cas9 genome editing, chromosome conformation capture assays (4C-seq, Capture-C) and Hi-C analysis of primary GCs. RESULTS: PeNChIP-seq/WGS revealed ~150 enhancer-based SVs in GC. The majority (63%) of SVs linked to target gene deregulation were associated with increased tumour expression. Enhancer-based SVs targeting CCNE1, a key driver of therapy resistance, occurred in 8% of patients frequently juxtaposing diverse distal enhancers to CCNE1 proximal regions. CCNE1-rearranged GCs were associated with high CCNE1 expression, disrupted CCNE1 topologically associating domain (TAD) boundaries, and novel TAD interactions in CCNE1-rearranged primary tumours. We also observed IGF2 enhancer-based SVs, previously noted in colorectal cancer, highlighting a common non-coding genetic driver alteration in gastric and colorectal malignancies. CONCLUSION: Integrated paired-end NanoChIP-seq and WGS of gastric tumours reveals tumour-associated regulatory SV in regions associated with both simple and complex genomic rearrangements. Genomic rearrangements may thus exploit enhancer-hijacking as a common mechanism to drive oncogene expression in GC.

SMuRF: portable and accurate ensemble prediction of somatic mutations.

SUMMARY: Somatic Mutation calling method using a Random Forest (SMuRF) integrates predictions and auxiliary features from multiple somatic mutation callers using a supervised machine learning approach. SMuRF is trained on community-curated matched tumor and normal whole genome sequencing data. SMuRF predicts both SNVs and indels with high accuracy in genome or exome-level sequencing data. Furthermore, the method is robust across multiple tested cancer types and predicts low allele frequency variants with high accuracy. In contrast to existing ensemble-based somatic mutation calling approaches, SMuRF works out-of-the-box and is orders of magnitudes faster. AVAILABILITY AND IMPLEMENTATION: The method is implemented in R and available at https://github.com/skandlab/SMuRF. SMuRF operates as an add-on to the community-developed bcbio-nextgen somatic variant calling pipeline. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Obstructive sleep apnea and Fuhrman grade in patients with clear cell renal cell carcinoma treated surgically.

PURPOSE: To assess the association between obstructive sleep apnea (OSA) and Fuhrman grade in patients with clear cell renal cell carcinoma (ccRCC). As secondary endpoints, we studied its association with tumor size, metastasis-free survival (MFS) and cancer-specific survival (CSS). METHODS: We reviewed the databases of two tertiary care centers, identifying 2579 patients who underwent partial or radical nephrectomy for ccRCC between 1991 and 2014. Descriptive statistics were used to compare pathologic variables between patients with and without OSA. Linear and logistic regression models were used to assess the association of OSA with Fuhrman grade and tumor size. A Cox proportional hazards model was used to determine OSA association with MFS and CSS. A pathway analysis was performed on a cohort with available gene expression data. RESULTS: In total, 172 patients (7 %) had self-reported OSA at diagnosis. More patients with OSA had high Fuhrman grade compared to those without OSA [51 vs. 38 %; 13 % risk difference; 95 % confidence interval (CI), 5-20 %; p = 0.003]. On multivariable analysis, the association remained significant (OR 1.41; 95 % CI 1.00-1.99; p = 0.048). OSA was not associated with tumor size (p > 0.5), MFS (p = 0.5) or CSS (p = 0.4). A trend toward vascular endothelial growth factor pathway enrichment was seen in OSA patients (p = 0.08). CONCLUSIONS: OSA is associated with high Fuhrman grade in patients undergoing surgery for ccRCC. Pending validation of this novel finding in further prospective studies, it could help shape future research to better understand etiological mechanisms associated.

Precision Endpoints for Contemporary Precision Oncology Trials.

SUMMARY: Traditional endpoints such as progression-free survival and overall survival do not fully capture the pharmacologic and pharmacodynamic effects of a therapeutic intervention. Incorporating mechanism-driven biomarkers and validated surrogate proximal endpoints can provide orthogonal readouts of anti-tumor activity and delineate the relative contribution of treatment components on an individual level, highlighting the limitation of solely relying on aggregated readouts from clinical trials to facilitate go/no-go decisions for precision therapies.

Transcriptome Deconvolution Reveals Absence of Cancer Cell Expression Signature in Immune Checkpoint Blockade Response.

Immune checkpoint therapy (ICB) has conferred significant and durable clinical benefit to some patients with cancer. However, most patients do not respond to ICB, and reliable biomarkers of ICB response are needed to improve patient stratification. Here, we performed a transcriptome-wide meta-analysis across 1,486 tumors from ICB-treated patients and tumors with expected ICB outcomes based on microsatellite status. Using a robust transcriptome deconvolution approach, we inferred cancer- and stroma-specific gene expression differences and identified cell-type specific features of ICB response across cancer types. Consistent with current knowledge, stromal expression of CXCL9, CXCL13, and IFNG were the top determinants of favorable ICB response. In addition, we identified a group of potential immune-suppressive genes, including FCER1A, associated with poor response to ICB. Strikingly, PD-L1 expression in stromal cells, but not cancer cells, is correlated with ICB response across cancer types. Furthermore, the unbiased transcriptome-wide analysis failed to identify cancer-cell intrinsic expression signatures of ICB response conserved across tumor types, suggesting that cancer cells lack tissue-agnostic transcriptomic features of ICB response. SIGNIFICANCE: Our results challenge the prevailing dogma that cancer cells present tissue-agnostic molecular markers that modulate immune activity and ICB response, which has implications on the development of improved ICB diagnostics and treatments.

Pushing the Frontiers of Cancer Research: Highlights from the Frontiers in Cancer Science Conference 2023.

The 15th annual Frontiers in Cancer Science (FCS) conference gathered scientific experts who shared the latest research converging upon several themes of cancer biology. These themes included the dysregulation of metabolism, cell death, and other signaling processes in cancer cells; using patient "omics" datasets and single-cell and spatial approaches to investigate heterogeneity, understand therapy resistance, and identify targets; innovative strategies for inhibiting tumors, including rational drug combinations and improved drug delivery mechanisms; and advances in models that can facilitate screening for cancer vulnerabilities and drug testing. We hope the insights from this meeting will stimulate further progress in the field.

Oncogene-Driven Non-Small Cell Lung Cancers in Patients with a History of Smoking Lack Smoking-Induced Mutations.

Non-small cell lung cancers (NSCLC) in nonsmokers are mostly driven by mutations in the oncogenes EGFR, ERBB2, and MET and fusions involving ALK and RET. In addition to occurring in nonsmokers, alterations in these "nonsmoking-related oncogenes" (NSRO) also occur in smokers. To better understand the clonal architecture and genomic landscape of NSRO-driven tumors in smokers compared with typical-smoking NSCLCs, we investigated genomic and transcriptomic alterations in 173 tumor sectors from 48 NSCLC patients. NSRO-driven NSCLCs in smokers and nonsmokers had similar genomic landscapes. Surprisingly, even in patients with prominent smoking histories, the mutational signature caused by tobacco smoking was essentially absent in NSRO-driven NSCLCs, which was confirmed in two large NSCLC data sets from other geographic regions. However, NSRO-driven NSCLCs in smokers had higher transcriptomic activities related to the regulation of the cell cycle. These findings suggest that, whereas the genomic landscape is similar between NSRO-driven NSCLC in smokers and nonsmokers, smoking still affects the tumor phenotype independently of genomic alterations. SIGNIFICANCE: Non-small cell lung cancers driven by nonsmoking-related oncogenes do not harbor genomic scars caused by smoking regardless of smoking history, indicating that the impact of smoking on these tumors is mainly nongenomic.

Liquid Biopsies for Cancer Genotyping: Coming of Age?

Precision oncology is predicated on optimal molecular profiling that is "fit for purpose" to identify therapeutic vulnerabilities. Liquid biopsies may compensate for inadequate genotyping, but remain less sensitive and specific compared with tissue biopsies. The liquid biopsy toolbox is poised to expand through novel assays and insights from longitudinal profiling. See related article by Sugimoto et al., p. 1506.

PUREE: accurate pan-cancer tumor purity estimation from gene expression data.

Tumors are complex masses composed of malignant and non-malignant cells. Variation in tumor purity (proportion of cancer cells in a sample) can both confound integrative analysis and enable studies of tumor heterogeneity. Here we developed PUREE, which uses a weakly supervised learning approach to infer tumor purity from a tumor gene expression profile. PUREE was trained on gene expression data and genomic consensus purity estimates from 7864 solid tumor samples. PUREE predicted purity with high accuracy across distinct solid tumor types and generalized to tumor samples from unseen tumor types and cohorts. Gene features of PUREE were further validated using single-cell RNA-seq data from distinct tumor types. In a comprehensive benchmark, PUREE outperformed existing transcriptome-based purity estimation approaches. Overall, PUREE is a highly accurate and versatile method for estimating tumor purity and interrogating tumor heterogeneity from bulk tumor gene expression data, which can complement genomics-based approaches or be used in settings where genomic data is unavailable.

Brief Report: Droplet Digital Polymerase Chain Reaction Versus Plasma Next-Generation Sequencing in Detecting Clearance of Plasma EGFR Mutations and Carcinoembryonic Antigen Levels as a Surrogate Measure.

Introduction: To compare the performance of droplet digital polymerase chain reaction (ddPCR) and plasma next-generation sequencing (NGS) in detecting clearance of plasma EGFR (pEGFR) mutations. Methods: Patients with treatment-naive advanced EGFR-mutated lung cancer treated with first-line tyrosine kinase inhibitors (TKIs) were included. pEGFR were measured at baseline and first response assessment using ddPCR and NGS. Clearance of pEGFR was defined as undetectable levels after a positive baseline result. Results were correlated with time-to-treatment failure (TTF). In exploratory analysis, corresponding change in carcinoembryonic antigen (CEA) levels was evaluated. Results: Between January 1, 2020, and December 31, 2021, 27 patients were recruited. Ex19del comprised 74% (20 of 27) and L858R 26% (seven of 27). Osimertinib was used in 59% (16 of 27), dacomitinib 4% (one of 27), and gefitinib/erlotinib 37% (10 of 27). Sensitivity of ddPCR and NGS in detecting pEGFR mutation at baseline was 70% (19 of 27) and 78% (21 of 27), respectively (p = 0.16). All patients with detectable pEGFR by ddPCR were detected by NGS.At a median of 8 (range 3-24) weeks post-TKI initiation, clearance of pEGFR was achieved in 68% (13 of 19) and 71% (15 of 21) using ddPCR and NGS, respectively. Concordance between ddPCR and NGS was 79% (kappa = 0.513, p = 0.013). Clearance of pEGFR was associated with longer median TTF (not reached versus 6 months, p = 0.03) and median decrease in CEA levels by 70% from baseline.In another cohort of 124 patients, decrease in CEA levels by greater than 70% within 90 days of TKI initiation was associated with doubling of both TTF and overall survival. Conclusions: Plasma NGS trended toward higher sensitivity than ddPCR in detecting pEGFR, although both had similar concordance in detecting pEGFR clearance. Our results support using NGS at diagnosis and interchangeability of NGS and ddPCR for monitoring, whereas CEA could be explored as a surrogate for pEGFR clearance.