Abstract 6689: Whole genome cell-free tumor DNA mutational signatures from blood for early detection of recurrence of low stage lung adenocarcinoma

Introduction: Lung cancer remains the leading cause of cancer-related deaths. Surgery is the best option for early lung cancer, and the role of adjuvant therapy remains controversial. Liquid biopsy offers a noninvasive approach to monitor cancer burden. Targeted sequencing of circulating cell free tumor DNA (ctDNA) in blood has shown success for diagnosis; however, low tumor burden and dynamic evolution of low stage disease is challenging for targeted panels. Thus, we hypothesized that a whole genome sequencing (WGS)-derived patient specific mutational signature from a matched tumor-normal WGS can provide sensitive and specific approach to detect mutations and copy numbers in ctDNA for monitoring of lung adenocarcinoma patients.

Methods: We successfully profiled 50 Stage 1 or 2 lung adenocarcinomas. ctDNA was extracted from 1-2 mL of plasma, tumor DNA was extracted from pathology tissue and normal germline DNA from the white blood cells. WGS using was performed on matched tumor and normal DNA, and ctDNA extracted from plasma. WGS coverage was 40x for matched tumor-normal and 20x for ctDNA. We derived a personalized mutational pattern for each tumor and used an AI-based error suppression model for quantification and ultra-sensitive detection of ctDNA in plasma samples. A patient-specific personalized genome-wide compendium of somatic mutations and copy numbers was established and ctDNA tested at 3 to 18 available time points during the therapy or follow up. A personalized mutational signature for detection ctDNA from WGS was quantified and the ctDNA Tumor Fraction (TF) was compared to the clinical status and time to recurrence.

Results: Tumor specific signatures were derived from matched tumor-normal samples with >5% tumor purity and <30% duplications rate. Out of all patients, 33 patients showed no recurrence and 12 recurred. Tumor-specific signatures detected the presence of the tumor signature in plasma with TF as low as 10−5. Based on positive minimal residual disease in plasma, the recurrence prediction sensitivity was 0.75 and specificity 0.82, with positive predictive value of 0.6 and negative predictive value 0.9. WGS ctDNA predicted recurrence with a median lead time of 508 days before clinical/imaging recurrence. In one case we were able to identify the second primary by deconvoluting known and novel ctDNA mutations. ctDNA mutational profiles enabled identification of smoking mutational signature matching clinical history, and APOBEC and ageing signatures as well as tumor mutational burden.

Conclusions: Patient-specific WGS tumor signature from plasma derived ctDNA enables specific and ultrasensitive tracking of minimal residual disease in low stage lung adenocarcinoma patients. Molecularly positive status can be used to predict recurrence and identify patients with clinical low stage disease that may benefit from adjuvant therapy.

Citation Format: Ivy Tran, Alejandro Vargas, Reid Wilkins, Isabella Pizzillo, Kenneth Tokoro, Danielle Afterman, Tomer Lauterman, Maja Kuzman, Santiago Gonzalez, Dunja Glavas, James Smadbeck, Dillon Maloney, Jurica Levatic, Samuel Phillips, Sunil Deochand, Michael Yahalom, Ryan Ptashkin, Iman Tavassoly, Zohar Donenhirsh, Eric White, Ravi Kandasamy, Ury Alon, Paz Polak, Boris Oklander, Asaf Zviran, Matija Snuderl, Harvey I. Pass. Whole genome cell-free tumor DNA mutational signatures from blood for early detection of recurrence of low stage lung adenocarcinoma. [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 6689.

Abstract 3401: Whole genome cell-free tumor DNA mutational signatures for noninvasive monitoring of pediatric brain cancers

Introduction: Liquid biopsy offers a noninvasive approach to monitor cancer burden during therapy and surveillance period. However, in pediatric brain cancers, liquid biopsy methods from the blood have been unsuccessful due to a low tumor burden and low number of mutations in coding regions. We hypothesized that a whole genome sequencing (WGS)-derived patient specific mutational signature from a matched tumor-normal WGS can provide a sensitive and specific approach to detect mutations in circulating cell free tumor DNA (ctDNA) and provide blood-based monitoring in pediatric patients with brain tumor.

Methods: All tumors were analyzed and molecularly subclassified using whole genome DNA methylation profiling and machine learning classifier. Tumor DNA was extracted from pathology tissue and normal germline DNA from the white blood cells, while ctDNA was extracted from 1-2 mL of post-surgery or follow-up plasma samples, WGS was applied to sequence DNA from matched tumor-normal and plasma samples. WGS coverage was 40x for matched tumor-normal DNA and 20x for ctDNA. Using the C2i assay, we derived a personalized mutational pattern for each tumor and used an AI-based error suppression model for quantification and ultra-sensitive detection of ctDNA in plasma samples. A patient-specific personalized genome-wide compendium of somatic mutations was established and ctDNA tested at 1 to 3 available time points during the therapy or surveillance period. An AI-based error suppression model was implemented to filter out the noise in the cell free DNA (cfDNA) while the personalized mutational signature was used to detect the ctDNA in the cfDNA and to amplify the somatic signal contained in it. The ctDNA Tumor Fraction (TF) was compared to the clinical status and MR-based imaging.

Results: We profiled 7 pediatric brain tumors, including 2 medulloblastomas (one Group 3, one Group 4), 3 pediatric glioblastomas IDH wild-type, 1 ependymoma PFA subtype and one low grade ganglioglioma. Tumor specific signatures were identified and detected in the plasma of 5 patients with clinical disease with a TF range 0.02-0.0005 but not in 2 patients with no tumor at the time of blood collection. In two children with a medulloblastoma and glioblastoma, the decrease of tumor fraction in ctDNA over 2 (TF: 0.002 to 0.0009) and 3 time points (TF: 0.0005 to undetectable), respectively, correlated with response to therapy based on imaging.

Conclusions: Patient-specific WGS tumor signature in ctDNA from blood can be used for sensitive monitoring of children with brain tumors.

Citation Format: Ivy Tran, Kristyn Galbraith, Guisheng Zhao, Robyn Borsuk, Joyce Varkey, Sharon Gardner, Jeffrey Allen, David Harter, Jeffrey Wisoff, Eveline T. Hidalgo, Sunil Deochand, Dillon Maloney, Danielle Afterman, Tomer Lauterman, Noah Friedman, Imane Bourzgui, Nidhi Ramaraj, Zohar Donenhirsh, Ronel Veksler, Jonathan Rosenfeld, Ravi Kandasamy, Iman Tavassoly, Boris Oklander, G. Praveen Raju, Theodore Nicolaides, Asaf Zviran, Matija Snuderl. Whole genome cell-free tumor DNA mutational signatures for noninvasive monitoring of pediatric brain cancers [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 3401.

Abstract 1959: Sensitive detection of circulating tumor DNA by whole genome sequencing: Validation of MRDetect using serial blood samples from stage III colorectal cancer patients

Background: While detection of circulating tumor DNA (ctDNA) is associated with poor cancer prognosis, the clinical utility for guiding treatment decisions is unresolved. Patients with minimal residual disease (MRD) often have less than one genome equivalent of ctDNA per 10 mL blood. Consequently, it is stochastic whether a 10 mL sample contains ctDNA from a particular genomic locus. Consequently, the sensitivity of ctDNA detection methods targeting a limited number of tumor loci is heavily affected by sampling bias. To overcome this challenge, we developed MRDetect; a whole genome sequencing (WGS) approach, which detects ctDNA using the patient-specific cumulative signal from tens of thousands of mutations throughout the genome. Recently, we showed how MRDetect found ctDNA fractions down to 10-4. Here, we performed a validation study to confirm the prognostic impact of MRDetect.

Aim: Validation of MRDetect for sensitive ctDNA detection to monitor residual disease in stage III colorectal cancer (CRC) patients treated with curative intent.

Methods: From a large, uniform cohort of stage III CRC patients n = 146), we had plasma samples collected every third month (n = 938, median = 9 per patient) and a median follow-up of 34 months. For each patient, a genome-wide mutational signature was established by WGS of tumor and matched normal DNA. Enhanced by an AI-based error suppression model, this signature was used to detect ctDNA in 1-2 mL plasma samples using WGS (20x coverage). We used de-novo point mutation and copy number variation analysis to investigate cancer evolution after treatment. To evaluate the reproducibility of MRDetect, aliquot samples (n = 2x190 samples) from 5 recurrence and 10 non-recurrence patients were processed and sequenced at two independent laboratories. Outcome measures: ctDNA status, tumor fraction, false positive rate, Time To ctDNA Recurrence (TTcR), and Time To radiological Recurrence (TTrR).

Results: Analysis of paired samples showed great reproducibility with high agreement between both ctDNA status calls (Cohens Kappa = 0.81) and the estimated tumor fractions (r2 = 0.99). MRDetect revealed post-operative ctDNA in all recurrence patients (5/5) with detected tumor fractions down to 2 x 10-4. Median TTcR was 0.9 month (range 0.5 - 7.3 months) while median TTrR was 12.8 months (range 11.3 - 31.1 months). The false positive rate was 1% (1/100), assessed in longitudinal samples from the 10 non-relapsing patients. Tumor evolution dynamics in plasma samples revealed novel amplification and deletions, which were absent in the primary tissue but confirmed in metachronous metastases. We will present results from the full cohort at AACR 2022.

Conclusion: MRDetect detects ctDNA with high sensitivity and specificity and enables effective postoperative assessment of MRD, cancer evolution dynamics and early relapse detection.

Citation Format: Amanda Frydendahl, Thomas Reinert, Jesper Nors, Sunil Deochand, Dillon Maloney, Noah Friedman, Tomer Lauterman, Danielle Afterman, Imane Bourzgui, Nidhi Ramaraj, Zohar Donenhirsh, Ronel Veksler, Ravi Kandasamy, Iman Tavassoly, Jonathan Rosenfeld, Anders Husted Andersen, Uffe S. Løve, Per V. Andersen, Ole Thorlacius-Ussing, Lene Hjerrild Iversen, Kåre Andersson Gotschalck, Boris Oklander, Asaf Zviran, Claus Lindbjerg Andersen. Sensitive detection of circulating tumor DNA by whole genome sequencing: Validation of MRDetect using serial blood samples from stage III colorectal cancer patients [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 1959.

Abstract 540: Genome-wide circulating tumor DNA for monitoring treatment response and metastatic relapse in bladder cancer

Background: Neoadjuvant chemotherapy (NAC) followed by radical cystectomy (CX), is gold standard treatment in localized muscle-invasive bladder cancer (MIBC). About 45% of patients with MIBC develop metastatic relapse within 2 years after CX. The response rate to chemotherapy and immune checkpoint inhibitors (ICI) is relatively low, and biomarker tests for monitoring response are needed. Furthermore, biomarkers for early detection of minimal residual disease (MRD) after CX is needed to enable earlier treatment initiation. Tumor-informed detection of mutations in cell-free DNA (cfDNA) from peripheral blood has shown promising results in its ability to monitor MRD. However, the low tumor fraction after surgery and limited input material obtained from a typical plasma sample limits the probability of detecting low metastatic burden scenarios. Here we implemented and applied locally a whole-genome sequencing (WGS) approach to circulating tumor DNA (ctDNA) monitoring for improving ctDNA detection.

Methods: A total of 140 MIBC patients undergoing NAC and CX were enrolled, including a test cohort (n=19) and a validation cohort (n=120). cfDNA was extracted from ~1mL plasma (n=1100) and procured from longitudinal plasma sampling during NAC (response measure), pre-cystectomy (response measure), post-surgery (relapse monitoring) and during immunotherapy (ICI treatment). WGS was applied to tumor/germline pairs (coverage >30x/20x) and plasma cfDNA (>20x) facilitating detection of genome wide genomic alterations and quantification of ctDNA using the MRDetect method.

Results: We developed a personalized tumor-informed WGS model by integrating genome-wide mutation and copy number variation data coupled with advanced signal processing and AI-based error suppression. Patient-specific somatic variant patterns were then used for detecting and measuring the ctDNA levels in low-input blood samples by WGS. The assay sensitivity allowed for detection of tumor fractions down to 8*10-5. Furthermore, in our test cohort of 19 patients, we detected ctDNA after CX in 7 of 8 patients with clinical relapse (88% sensitivity) and detected no ctDNA in 11 of 11 patients with no clinical relapse (100% specificity). We observed a positive lead-time for MRD-based recurrence detection compared to CT-based reccurence detection (9 months on average). The full dataset is currently being processed and will be presented at the AACR 2022 meeting.

Conclusions: For precision oncology, we need to develop quantitative and non-invasive methodologies to help tailor the treatments to individual patients and monitor them for further clinical decision-making. The results indicate the clinical potential of personalized genome-wide mutation integration as an ultra-sensitive, non-invasive method for MRD detection and treatment response monitoring which could aid in clinical management of patients with bladder cancer.

Citation Format: Iver Nordentoft, Karin Birkenkamp-Demtröder, Emil Christensen, Sunil Deochand, Dillon Maloney, Danielle Afterman, Tomer Lauterman, Noah Friedman, Imane Bourzgui, Nidhi Ramaraj, Zohar Donenhirsh, Ronel Veksler, Sia Viborg, Mads Agerbæk, Jørgen Bjerggaard Jensen, Jonathan Rosenfeld, Ravi Kandasamy, Iman Tavassoly, Boris Oklander, Asaf Zviran, Lars Dyrskjøt. Genome-wide circulating tumor DNA for monitoring treatment response and metastatic relapse in bladder cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 540.

Abstract 5114: Ultra-sensitive detection of minimal residual disease (MRD) through whole genome sequencing (WGS) using an AI-based error suppression model in resected early-stage non-small cell lung cancer (NSCLC)

Background: Early detection of recurrence and monitoring of MRD post-surgery is critical for clinical decision-making to tailor adjuvant therapy. In early-stage NSCLC, circulating tumor DNA (ctDNA) detection is especially challenging, requiring highly sensitive and specific assays. Therefore, we used a WGS approach (MRDetect) for ultra-sensitive ctDNA detection in NSCLC patients (pts) undergoing curative surgery.

Methods: We conducted a pilot study to evaluate the MRDetect approach in serial plasma samples (including pre-surgery, post-surgery and follow-up [f/u] timepoints) from resected stage IB-IIIA NSCLC pts. Pts underwent routine surveillance by computed tomography scans. ctDNA was extracted from ~1mL plasma. MRDetect uses WGS by a tumor-informed approach (sequencing coverage 40x for tumor, 20x for plasma DNA) combined with AI-based error suppression models (trained and calibrated with a non-cancer cohort, n=17) to increase the signal to noise ratio for precise ctDNA detection, and improve the accuracy of readouts especially for low tumor burden scenarios. The assay reports the detection and quantification of ctDNA burden in blood with a prognostic value for risk of recurrence. The ability of the assay to predict recurrence from a single sample, taken at the clinical landmark point (median 1.6 mths post-surgery, range 0.1-6.5) was evaluated.

Results: Overall, 52 NSCLC pts were enrolled (n=88 plasma samples) with median clinical f/u of 32.6 mths (range 3.1-98.6). There were 43 pts with post-surgery landmark samples, with median age 62 years, 70% were male, 79% were adenocarcinoma and 49% were EGFR mutated. 26% were stage IB and 37% each were stage II and III. There were 15/18 (sensitivity 83%) pts with confirmed radiological recurrence in which MRDetect was positive, including 6/7 (86%) EGFR mutated pts. The median RFS in MRDetect positive pts was 15.2 mths (range 3.7-33.4). Among 25 pts with no recurrence (median f/u 25.6 mths), MRDetect reported 4 pts to be MRD positive (specificity 84%). These results were consistent between EGFR mutated (sensitivity 86%, specificity 86%) and wildtype pts (sensitivity 82%, specificity 82%). For longitudinal samples (n=17 pts), negative ctDNA was associated with absence of recurrence in 14/15 pts (specificity 93%). At the AACR meeting, results from a planned larger validation study will be presented.

Conclusion: Using a robust WGS implemented AI-based computational platform (MRDetect), we demonstrate high sensitivity and specificity detection of MRD in both EGFR mutated and wildtype NSCLC. With an increasing number of therapeutic options in the adjuvant setting for NSCLC, an ultra-sensitive MRD assay has the potential to facilitate personalized clinical decision-making for tailoring both the need and choice of adjuvant therapies.

Citation Format: Aaron C. Tan, Stephanie P. Saw, Gillianne G. Lai, Kevin L. Chua, Angela Takano, Boon-Hean Ong, Tina P. Koh, Amit Jain, Wan Ling Tan, Quan Sing Ng, Ravindran Kanesvaran, Tanujaa Rajasekaran, Sunil Deochand, Dillon Maloney, Danielle Afterman, Tomer Lauterman, Noah Friedman, Imane Bourzgui, Nidhi Ramaraj, Zohar Donenhirsh, Ronel Veksler, Jonathan Rosenfeld, Ravi Kandasamy, Iman Tavassoly, Boris Oklander, Asaf Zviran, Wan-Teck Lim, Eng-Huat Tan, Anders J. Skanderup, Mei-Kim Ang, Daniel S. Tan. Ultra-sensitive detection of minimal residual disease (MRD) through whole genome sequencing (WGS) using an AI-based error suppression model in resected early-stage non-small cell lung cancer (NSCLC) [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 5114.

Lab validation of an ultrasensitive ctDNA pan-cancer MRD assay using whole-genome sequencing

e13582

Background: Minimal residual disease (MRD) monitoring using liquid biopsy for solid tumors requires a highly sensitive and specific assay that can overcome the limitation of low abundance cfDNA in a standard blood draw. We developed a whole-genome sequencing (WGS)-based assay to detect the presence of circulating tumor DNA (ctDNA) in plasma. The C2i assay is a tumor-informed assay that uses personalized tumor signature, advanced noise models, and artificial intelligence (AI) modalities to interrogate plasma for the presence of ctDNA longitudinally. Methods: The C2i test was developed in accordance with CAP/CLIA and New York state validation principles. We used contrived samples to establish analytical validation of the assay performance, which was then validated with a large clinical cohort of early-stage patients across various cancer types. Briefly, aggregated tumor signatures derived from cancer cell lines were fragmented and spiked into a contrived healthy plasma pool; the mixed samples were used to assess the presence of tumor DNA signature down to tumor fractions of 10e-4. Positive samples are identified by tumor-derived variants detected above the noise levels. Noise modeling was established using a panel of normal (PON) approach. We assessed the analytical sensitivity, specificity, and accuracy using 348 contrived samples derived from five different cancer cell lines. Reproducibility and precision were assessed with multiple replicates, and statistical concordance was reported. This validation was complemented by a cohort of 200 patients and ̃1000 plasma samples across a variety of cancer types including, NSCLC, MIBC, CRC, GBM, Breast Cancer, and a mixture of other cancer types. Results: Cancer cell lines, representing the five most prevalent disease indications, used for determining analytical sensitivity are as follows: CRC HT-29, Breast SK-BR3, Bladder HT-1376, Lung HCI-H526, and Prostate LNCaP. The cell line DNA was enzymatically fragmented and size-selected to mimic ctDNA. This ctDNA was spiked into cfDNA extracted from healthy volunteers at various dilution levels, varying from 10e0 to 10e-4. The 95% probability of detecting ctDNA was established at 10e-4. The reproducibility of tumor signature between replicates was assessed to be greater than 90%. The assay was performed using both normal and maximum input amounts. These performance estimates were then validated on a cohort of plasma collected from early-stage (stage I-III) patients across various cancer types. Conclusions: C2i MRD test is an ultrasensitive pan-cancer MRD monitoring assay used in several clinical trials across the world. We present an extensive analytical and clinical validation of the assay supporting its high performance.

BIOM-23. A PILOT STUDY OF WHOLE GENOME SEQUENCING (WGS) OF PLASMA CELL-FREE DNA (cfDNA) FOR ULTRASENSITIVE DETECTION OF TUMOR DNA IN PATIENTS WITH GLIOBLASTOMA (GBM)

BACKGROUND

Plasma circulating tumor DNA (ctDNA) is rarely detectable by traditional methods in patients with GBM. As a result, unlike in lung and other cancers, serial next generation sequencing of ctDNA for monitoring GBM tumor burden has been challenging. In light of the low tumor fraction (TF) of DNA fragments in GBM patient plasma and the urgent need to improve upon MRI for tracking GBM tumor burden, we conducted a pilot study in patients with newly diagnosed GBM using the C2 intelligence platform (C2i Genomics), which leverages genome-wide mutational integration for highly sensitive ctDNA detection.

METHODS

Plasma was collected pre- and post-operatively in patients with newly diagnosed GBM undergoing surgical resection/biopsy. cfDNA was extracted, quantified, and analyzed for fragment size. Genomic DNA (gDNA) was extracted from matched tumor tissue. Whole genome sequencing (WGS) was performed on both gDNA and cfDNA. A specific copy number alteration (CNA) compendium was created for each patient to generate a readout of TF (Zviran, Nat Medicine 2020). We assessed the association between TF at post-operative day 1 (a surrogate for residual disease) and OS, adjusting for other prognostic factors using Cox regression.

RESULTS

37 patients were enrolled. For samples with high tumor fraction (n=5), a statistically significant (p< 1e-4) correlation between CNA profiles of tumor tissue and plasma samples was observed. Post-operative TF above the median value was associated with inferior OS (median 7.7 vs. 19.3 months, p=0.019). This association persisted after adjusting for age, O6-methylguanine-DNA methyltransferase methylation status, extent of resection, and performance status (adjusted HR 2.5, 95% CI 1.1-5.6, p=0.03).

CONCLUSION

Genome-wide mutational integration enables ultra-sensitive detection of ctDNA in GBM patient plasma. Post-operative TF measured by the C2i test is independently associated with OS in newly diagnosed GBM, providing the foundation to evaluate this technology for personalized prognostication and disease monitoring.

Genome-wide circulating tumor DNA monitoring for bladder cancer treatment management and organ preservation

e16527

Background: Bladder cancer (BC) is the 9th most commonly diagnosed cancer worldwide and each year responsible for 165,000 deaths. Neoadjuvant combination chemotherapy, followed by radical cystectomy, is used for the management of localized muscle-invasive bladder cancer. One of the critical challenges in this therapeutic regimen is monitoring the tumor load to assess therapeutic efficacy – this is typically performed by assessing pathological downstaging in the cystectomy specimen. A high frequency of patients presents with T0N0 at cystectomy (no indication of residual disease), and consequently, it is vital to investigate organ preservation approaches to identify those patients who may qualify for bladder preservation. For precision oncology, we need to develop quantitative and non-invasive diagnostic methodologies to help the oncologist tailor the treatments to individual patients and monitor them for further clinical decision-making. Methods: Cell-free DNA (cfDNA) mutation detection has shown significant promise in its ability to monitor minimal residual disease and disease relapse by detection of cancer mutations in the peripheral blood. However, the combination of low tumor fraction and limited input material obtained from a typical plasma sample restricts the probability of detecting low metastatic burden in cfDNA through current deep targeted sequencing methods. Results: Here we present results from applying whole-genome sequencing (WGS) of cfDNA. We integrate a genome-wide mutation and copy number monitoring approach coupled with advanced signal processing and Artificial Intelligence (AI) for measuring the tumor load from low-input blood samples (̃1mL of plasma) with ultra-sensitive detection. The increased sensitivity allowed clinical detection of tumor fraction down to 8*10-5 and recurrence detection sensitivity achieving > 65% at the first two months post-surgery. The WGS cfDNA approach is being evaluated on a patient cohort of more than 50 bladder cancer patients with longitudinal plasma sampling during neoadjuvant chemotherapy (response measure), pre-cystectomy (complete response measure), and post-surgery (relapse monitoring). Conclusions: The results indicate the clinical potential of genome-wide mutation integration as an ultra-sensitive, non-invasive diagnostic method for bladder cancer clinical management and bladder organ preservation.

Acute Myeloid Leukemia iPSCs Reveal a Role for RUNX1 in the Maintenance of Human Leukemia Stem Cells

Leukemia stem cells (LSCs) are believed to have more distinct vulnerabilities than the bulk acute myeloid leukemia (AML) cells, but their rarity and the lack of universal markers for their prospective isolation hamper their study. We report that genetically clonal induced pluripotent stem cells (iPSCs) derived from an AML patient and characterized by exceptionally high engraftment potential give rise, upon hematopoietic differentiation, to a phenotypic hierarchy. Through fate-tracking experiments, xenotransplantation, and single-cell transcriptomics, we identify a cell fraction (iLSC) that can be isolated prospectively by means of adherent in vitro growth that resides on the apex of this hierarchy and fulfills the hallmark features of LSCs. Through integrative genomic studies of the iLSC transcriptome and chromatin landscape, we derive an LSC gene signature that predicts patient survival and uncovers a dependency of LSCs, across AML genotypes, on the RUNX1 transcription factor. These findings can empower efforts to therapeutically target AML LSCs.

Genome-wide cell-free DNA mutational integration enables ultra-sensitive cancer monitoring

In many areas of oncology, we lack sensitive tools to track low-burden disease. Although cell-free DNA (cfDNA) shows promise in detecting cancer mutations, we found that the combination of low tumor fraction (TF) and limited number of DNA fragments restricts low-disease-burden monitoring through the prevailing deep targeted sequencing paradigm. We reasoned that breadth may supplant depth of sequencing to overcome the barrier of cfDNA abundance. Whole-genome sequencing (WGS) of cfDNA allowed ultra-sensitive detection, capitalizing on the cumulative signal of thousands of somatic mutations observed in solid malignancies, with TF detection sensitivity as low as 10-5. The WGS approach enabled dynamic tumor burden tracking and postoperative residual disease detection, associated with adverse outcome. Thus, we present an orthogonal framework for cfDNA cancer monitoring via genome-wide mutational integration, enabling ultra-sensitive detection, overcoming the limitation of cfDNA abundance and empowering treatment optimization in low-disease-burden oncology care.

Genotype-Fitness Maps of EGFR-Mutant Lung Adenocarcinoma Chart the Evolutionary Landscape of Resistance for Combination Therapy Optimization

Cancer evolution poses a central obstacle to cure, as resistant clones expand under therapeutic selection pressures. Genome sequencing of relapsed disease can nominate genomic alterations conferring resistance but sample collection lags behind, limiting therapeutic innovation. Genome-wide screens offer a complementary approach to chart the compendium of escape genotypes, anticipating clinical resistance. We report genome-wide open reading frame (ORF) resistance screens for first- and third-generation epidermal growth factor receptor (EGFR) inhibitors and a MEK inhibitor. Using serial sampling, dose gradients, and mathematical modeling, we generate genotype-fitness maps across therapeutic contexts and identify alterations that escape therapy. Our data expose varying dose-fitness relationship across genotypes, ranging from complete dose invariance to paradoxical dose dependency where fitness increases in higher doses. We predict fitness with combination therapy and compare these estimates to genome-wide fitness maps of drug combinations, identifying genotypes where combination therapy results in unexpected inferior effectiveness. These data are applied to nominate combination optimization strategies to forestall resistant disease.

Neutralizing Gatad2a-Chd4-Mbd3/NuRD Complex Facilitates Deterministic Induction of Naive Pluripotency

Mbd3, a member of nucleosome remodeling and deacetylase (NuRD) co-repressor complex, was previously identified as an inhibitor for deterministic induced pluripotent stem cell (iPSC) reprogramming, where up to 100% of donor cells successfully complete the process. NuRD can assume multiple mutually exclusive conformations, and it remains unclear whether this deterministic phenotype can be attributed to a specific Mbd3/NuRD subcomplex. Moreover, since complete ablation of Mbd3 blocks somatic cell proliferation, we aimed to explore functionally relevant alternative ways to neutralize Mbd3-dependent NuRD activity. We identify Gatad2a, a NuRD-specific subunit, whose complete deletion specifically disrupts Mbd3/NuRD repressive activity on the pluripotency circuitry during iPSC differentiation and reprogramming without ablating somatic cell proliferation. Inhibition of Gatad2a facilitates deterministic murine iPSC reprogramming within 8 days. We validate a distinct molecular axis, Gatad2a-Chd4-Mbd3, within Mbd3/NuRD as being critical for blocking reestablishment of naive pluripotency and further highlight signaling-dependent and post-translational modifications of Mbd3/NuRD that influence its interactions and assembly.

Abstract 3247: Genome-wide cell-free DNA mutation integration for sensitive cancer detection

Solid malignancies are often diagnosed at a late stage with dismal prognosis. Even after cancer is diagnosed, we lack sensitive tools to guide difficult therapeutic decisions such as adjuvant therapy. Sensitive cancer detection by blood biopsy can therefore transform care by enabling early detection and residual disease monitoring.

Cell free DNA mutation detection has shown significant promise in its ability to survey the somatic genome and enable detection of cancer mutations in the peripheral blood. However, the combination of low tumor fraction and limiting number of DNA fragments in a typical plasma sample, restrict the probability of detecting early stage cancer in cfDNA through current deep targeted sequencing methods.

Focusing on non-small cell lung cancer (NSCLC), we reasoned that we would need to supplant depth of sequencing with breadth of sequencing to overcome the fundamental limitation of low input of cfDNA. To do so, we apply whole genome sequencing (WGS) that allows us to base sensitive detection on the cumulative signal provided by 10,000-30,000 somatic mutations observed in a substantial proportion of NSCLC. We developed an analytic method that integrates genome-wide mutation signal to obtain a tumor fraction (TF) estimate, and thus allow sensitive detection of residual disease and quantitative dynamic monitoring of disease burden. Benchmarking on artificial plasma showed TF detection sensitivity as low as 1:100,000, two orders of magnitude more sensitive than currently available methods.

To test this method, we performed WGS on resected NSCLC and matched germline samples of 8 NSCLC patients, as well as on matched pre- and post-surgery cfDNA. Patient-specific somatic mutations were identified in the tumor/normal pairs and used for the estimation of TF in the matched plasma samples. We detect pre-surgery circulating tumor DNA (ctDNA) in all of the early-stage pre-operative samples and in ~40% post-operative patients, correlated with post-operative disease progression.

In early cancer detection, tumor DNA is not available, requiring de-novo mutation detection in cfDNA. To do so, we first trained a convolutional neuronal network to distinguish between cancer altered sequencing reads and reads affected by sequencing errors. This was followed by genome-wide pattern matching to a specific genomic signature that mark lung cancer mutations (Tobacco signature) indicating the presence of ctDNA in the patient plasma. Applying this method to the pre-operative early stage lung cancer samples and plasma samples from 5 patients with benign nodules (CT-detected) showed an accurate discrimination between malignant and benign nodules, suggesting a potential role in improving the positive predictive value of lung cancer screening in at-risk populations.

These results show that genome-wide mutation integration is a promising novel approach for ultra-sensitive early detection and residual disease monitoring.

Citation Format: Asaf Zviran, Steven T. Hill, Rafael Schulman, Minita Shah, Sunil Deochand, Gavin Ha, Sarah Reed, Denisse Rotem, Greg Gydush, Justin Rhoades, Kevin Huang, Will Liao, Dillon Maloney, Nathan Omans, Murtaza Malbari, Cathy F. Spinelli, Selena Kazancioglu, Nicolas Robine, Viktor Adalsteinsson, Brian Houck-Loomis, Nasser Altorki, Dan A. Landau. Genome-wide cell-free DNA mutation integration for sensitive cancer detection [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 3247.

Abstract 1178: Genotype-fitness maps guide targeted therapy combination in lung adenocarcinoma

Targeted EGFR inhibition in lung cancer leads to dramatic responses. Nonetheless, disease evolution to resistance is the rule. As this evolutionary process is fueled by intra-tumoral genetic diversity, a comprehensive mapping of clonal fitness is required to inform strategies to overcome resistance.

To model intra-tumoral diversity in vitro, we performed a genome-wide, over-expression genetic perturbation assay in EGFR-driven non-small cell lung cancer (NSCLC) PC9 cells. Our screen covered 17,255 ORFs (open reading frame constructs) representing 12,728 wildtype and mutated genes, and examined the effects of a first-generation EGFR inhibitor (erlotinib), a third-generation EGFR inhibitor (osimertinib) and a MEK inhibitor (binimetinib) on evolutionary selection, alone or in combination. Specifically, to obtain genotype-to-fitness maps for each drug, we measured clonal abundance up to 4 times during the screen and mathematically resolved their growth behavior. Finally, for each drug, we applied multiple doses to obtain dose-fitness relationships, which may impact tumor evolution in patients due to high inter- and intra-tumoral variability in drug delivery. In total, we have performed 78 genome-wide screens to map the evolutionary landscape of PC9 resistance to targeted therapy.

Erlotinib and osimertinib both result in an overwhelming reduction of fitness, as expected from their clinical benefit. Known clinical resistance mechanisms involving ERBB2, PIK3CA, AXL, and BRAF confer a pronounced fitness advantage in the presence of both drugs, while EGFR (T790M) improves fitness only with erlotinib. Novel resistance mechanisms include alternative tyrosine kinases (NTRK, PDGFRB, CSF1R, KIT), KRAS, G-coupled protein receptors, transcription factors (SOX15, FOXA1), and cellular transporters (ABCG2). Notably, we observe significant divergence in dose-fitness relationships between different resistance mechanisms. For example, PIK3CA confers a modest but persistent advantage across the dose range, in contrast to BRAF which results in a fitness advantage only when sufficiently high drug levels are added.

The fitness landscape for binimetinib resistance appeared to be largely orthogonal to that of the EGFR inhibitors, with decreased fitness noted for many tyrosine kinases as well as KRAS mutations. MEK inhibitor resistance results from RAF family member overexpression (ARAF, BRAF, RAF1) and MAPK activation. The orthogonal resistance landscapes of EGFR inhibition and MEK inhibition translate into highly synergistic effects in drug combination, with effective prediction of combinatorial fitness from the single-agent fitness landscapes using generalized linear models. These results validate this approach as a systematic method to address the combinatorial problem of optimizing drug combinations and doses to directly anticipate and address cancer evolution.

Citation Format: Asaf Zviran, Patrick Bolan, Lisa Brenan, Amy Goodale, Denisse Rotem, Viktor Adalsteinsson, Federica Piccioni, Cory Johannessen, Dan Landau. Genotype-fitness maps guide targeted therapy combination in lung adenocarcinoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 1178.

The evolutionary landscape of chronic lymphocytic leukemia treated with ibrutinib targeted therapy

Treatment of chronic lymphocytic leukemia (CLL) has shifted from chemo-immunotherapy to targeted agents. To define the evolutionary dynamics induced by targeted therapy in CLL, we perform serial exome and transcriptome sequencing for 61 ibrutinib-treated CLLs. Here, we report clonal shifts (change >0.1 in clonal cancer cell fraction, Q < 0.1) in 31% of patients during the first year of therapy, associated with adverse outcome. We also observe transcriptional downregulation of pathways mediating energy metabolism, cell cycle, and B cell receptor signaling. Known and previously undescribed mutations in BTK and PLCG2, or uncommonly, other candidate alterations are present in seventeen subjects at the time of progression. Thus, the frequently observed clonal shifts during the early treatment period and its potential association with adverse outcome may reflect greater evolutionary capacity, heralding the emergence of drug-resistant clones.

In a subset of patients with chronic lymphocytic leukemia (CLL) treated with targeted agents, such as ibrutinib, drug resistant subclones emerge. Here, the authors report on transcriptional changes in CLL patients treated with ibrutinib and identify early clonal shifts associated with evolution of resistant clones.

Transient acquisition of pluripotency during somatic cell transdifferentiation with iPSC reprogramming factors

Somatic cells can be transdifferentiated to other cell types without passing through a pluripotent state by ectopic expression of appropriate transcription factors^1,2. Recent reports have proposed an alternative transdifferentiation method in which fibroblasts are directly converted to various mature somatic cell types by brief expression of the induced pluripotent stem cell (iPSC) reprogramming factors Oct4, Sox2, Klf4 and c-Myc (OSKM) followed by cell expansion in media that promote lineage differentiation^3–6. Here we test this method using genetic lineage tracing for expression of endogenous Nanog and Oct4 and for X chromosome reactivation, as these events mark acquisition of pluripotency. We show that the vast majority of reprogrammed cardiomyocytes or neural stem cells obtained from mouse fibroblasts by OSKM-induced transdifferentiation pass through a transient pluripotent state, and that their derivation is molecularly coupled to iPSC formation mechanisms. Our findings underscore the importance of defining trajectories during cell reprogramming by different methods.

Lucky iPSCs

The probabilistic behavior of direct induction of pluripotency has been a subject of intense research interest. Here we discuss recently published reports on this topic.

Derivation of novel human ground state naive pluripotent stem cells

Mouse embryonic stem (ES) cells are isolated from the inner cell mass of blastocysts, and can be preserved in vitro in a naive inner-cell-mass-like configuration by providing exogenous stimulation with leukaemia inhibitory factor (LIF) and small molecule inhibition of ERK1/ERK2 and GSK3β signalling (termed 2i/LIF conditions). Hallmarks of naive pluripotency include driving Oct4 (also known as Pou5f1) transcription by its distal enhancer, retaining a pre-inactivation X chromosome state, and global reduction in DNA methylation and in H3K27me3 repressive chromatin mark deposition on developmental regulatory gene promoters. Upon withdrawal of 2i/LIF, naive mouse ES cells can drift towards a primed pluripotent state resembling that of the post-implantation epiblast. Although human ES cells share several molecular features with naive mouse ES cells, they also share a variety of epigenetic properties with primed murine epiblast stem cells (EpiSCs). These include predominant use of the proximal enhancer element to maintain OCT4 expression, pronounced tendency for X chromosome inactivation in most female human ES cells, increase in DNA methylation and prominent deposition of H3K27me3 and bivalent domain acquisition on lineage regulatory genes. The feasibility of establishing human ground state naive pluripotency in vitro with equivalent molecular and functional features to those characterized in mouse ES cells remains to be defined. Here we establish defined conditions that facilitate the derivation of genetically unmodified human naive pluripotent stem cells from already established primed human ES cells, from somatic cells through induced pluripotent stem (iPS) cell reprogramming or directly from blastocysts. The novel naive pluripotent cells validated herein retain molecular characteristics and functional properties that are highly similar to mouse naive ES cells, and distinct from conventional primed human pluripotent cells. This includes competence in the generation of cross-species chimaeric mouse embryos that underwent organogenesis following microinjection of human naive iPS cells into mouse morulas. Collectively, our findings establish new avenues for regenerative medicine, patient-specific iPS cell disease modelling and the study of early human development in vitro and in vivo.

Deterministic direct reprogramming of somatic cells to pluripotency

Somatic cells can be inefficiently and stochastically reprogrammed into induced pluripotent stem (iPS) cells by exogenous expression of Oct4 (also called Pou5f1), Sox2, Klf4 and Myc (hereafter referred to as OSKM). The nature of the predominant rate-limiting barrier(s) preventing the majority of cells to successfully and synchronously reprogram remains to be defined. Here we show that depleting Mbd3, a core member of the Mbd3/NuRD (nucleosome remodelling and deacetylation) repressor complex, together with OSKM transduction and reprogramming in naive pluripotency promoting conditions, result in deterministic and synchronized iPS cell reprogramming (near 100% efficiency within seven days from mouse and human cells). Our findings uncover a dichotomous molecular function for the reprogramming factors, serving to reactivate endogenous pluripotency networks while simultaneously directly recruiting the Mbd3/NuRD repressor complex that potently restrains the reactivation of OSKM downstream target genes. Subsequently, the latter interactions, which are largely depleted during early pre-implantation development in vivo, lead to a stochastic and protracted reprogramming trajectory towards pluripotency in vitro. The deterministic reprogramming approach devised here offers a novel platform for the dissection of molecular dynamics leading to establishing pluripotency at unprecedented flexibility and resolution.